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In  this  dissertation,  generation-recombination  (G-R)  noise  is 

used  as  a tool  for  investigating  the  properties  of  deep  defects.  The  first 

experiment  we  describe  deals  with  low-field  noise  measurements, 

which  were  done  on  p+-p"-p+  silicon  devices  doped  with  a boron 

concentration  of  10  cm"  and  an  unknown  gold  concentration.  The 

resistance  and  the  G-R  noise  spectral  density  of  the  device  were 

measured  as  a function  of  temperature  and  used  to  extract  the  trap 

energy,  hole  capture  cross  section,  degeneracy  factor,  and  trap 

density.  The  resistance  and  the  noise  data  show  the  same  activation 

energy  of  0.336  eV,  which  stems  from  the  generally  accepted  gold 

donor  level,  hereafter  referred  to  as  the  gold  first-donor  level.  The 

values  of  trap  parameters  so  obtained  are  in  good  agreement  with  the 
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published  data  in  the  literature,  indicating  that  this  method  is 
accurate  for  investigating  deep-defect  properties. 

The  second  experiment  described  in  this  thesis  was  done  on 
p+-p"-p+  silicon  devices  doped  with  a boron  concentration  of  2.7 
xlO  cm  . The  resistance  and  low  frequency  voltage  noise  of  the 
devices  were  measured  as  a function  of  temperature.  The  resistance 
data  indicate  that  gold  produces  an  activation  energy  of  0.267  eV,  and 
the  noise  data  indicate  that  gold  produces  an  activation  energy  of 
0.336  eV.  This  observation  is  explained  in  terms  of  double  donor 
centers  interacting  with  the  valence  band. 

Then  the  current-voltage  (I-V)  characteristics  of  the  devices 
used  in  the  first  experiment  were  measured.  They  showed  sublinear 
behavior  at  moderately  high  fields.  With  the  use  of  low  frequency 
voltage  noise  data  measured  as  a function  of  electric  field,  we 
determined  that  the  sublinear  I-V  behavior  was  caused  by  a decrease 
in  the  hole  concentration  with  increasing  electric  field.  This  decrease 
is  caused  by  a field-induced  decrease  in  the  hole  emission  coefficient, 
which  is  larger  than  the  associated  field-induced  decrease  in  the  hole 
capture  coefficient. 

Device  design  schemes  for  investigating  deep  defects  by  noise 
measurements  are  presented.  We  analyzed  the  noise  of  several  multi- 
level trap  centers  and  derived  the  conditions  for  which  the  noise  of  a 
multi-level  trap  center  reduces  to  the  noise  of  a single-level  trap 
center.  We  also  devised  a method  by  which  G-R  noise  measurements 
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can  be  used  to  discriminate  between  donor-like  and  acceptor-like  trap 
behavior. 

In  addition  to  the  above  described  macroscopic  G-R  noise 
models,  a microscopic  model  for  carrier  capture  in  deep  traps  is 
derived.  Finally,  theories  for  G-R  noise  generated  in  field-effect 
devices  and  polysilicon  resistors  are  presented. 
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CHAPTER  I 
INTRODUCTION 


All  solids  in  equilibrium  and  at  finite  temperatures  contain 
intrinsic  point  defects.  In  addition,  all  solids  inevitably  contain 
chemical  impurities.  The  presence  of  defects  can  affect  the  usefulness 
of  materials  and  performance  of  solid-state  devices  in  both  desirable 
and  undesirable  ways.  Defect  study  in  semiconductor  engineering 
started  originally  with  seeking  methods  to  suppress  generation  of 
harmful  defects  during  device  processing  in  order  to  achieve  a high 
yield  of  device  fabrication.  Then  a new  trend  developed  in  which 
crystal  defects  are  utilized  to  improve  the  device  performance  and 
reliability.  A typical  example  is  using  a deep  level  impurity  to  control 
the  lifetime  of  carriers.  The  latter  trend  accounts,  to  a considerable 
extent,  for  the  intensive  and  continuing  research  activity  in  the  defect 
solid  state  area. 

Defect  engineering  becomes  increasingly  important  for  two 
reasons:  (1)  Undesirable  defects  cause  nothing  but  trouble;  they  make 
devices  inoperable  or  accelerate  the  degradation  of  device 
performance.  The  trend  towards  VLSI  submicron  semiconductor 
circuitry  has  increased  the  vulnerability  of  devices  for  unwanted 
process-induced  defects.  (2)  Defects  regulate  many  optical  and 
electrical  properties  of  interest.  The  device  performance  can  be 
improved  by  controlling  doping  profiles  and  impurity  concentrations. 
Isovalent  impurities  [1]  can  be  used  to  make  radiative  recombination 


1 


2 


possible  in  an  indirect  band  gap  semiconductor,  which  thus  becomes  a 
material  suitable  for  light  emitting  diode  applications.  In  fast  switching 
applications,  efficient  recombination  mechanisms  are  necessary  to 
remove  mobile  carriers  quickly.  Deep  level  defects  with  large  hole  and 
electron  capture  coefficient  are  good  candidates  for  this  application. 

Gold  has  been  used  as  a desirable  impurity  in  silicon  device 
technology  to  reduce  carrier  lifetime.  Despite  the  fact  that  it  has 
received  much  attention,  its  behavior,  once  incorporated,  is  not  always 
predictable.  Since  different  techniques  measure  defect  properties 
under  different  conditions,  applying  various  techniques  to  characterize 
gold  in  silicon  will  provide  a complete  picture  of  this  defect. 

In  this  thesis,  the  term  noise  is  used  to  describe  fluctuations  in 
an  electronic  signal  measured  on  an  electronic  device.  Noise  arises 
from  the  statistical  nature  of  the  electron  and/or  phonon  processes 
which  affect  charge  transport.  Since  it  is  associated  with  the 
microscopic  motion  of  carriers  in  a crystal,  noise  measurements  can 
be  used  to  investigate  for  example  the  properties  of  defects  in  a 
semiconductor  host  as  we  will  show  in  this  thesis. 

1.1.  Basic  Concepts  of  Defects  and  Impurities  in  Solids 

A perfect  crystal  consists  of  a three-dimensional  array  of  atoms 
arranged  on  a periodic  lattice.  The  intentional  or  unintentional 
introduction  of  native  defects  and  impurities  in  such  a crystal  disrupts 
this  periodic  structure  and  may  alter  the  properties  of  the  material  in 
a significant  way.  The  term  "flaw"  or  "imperfection,"  hereinafter 
referred  to  as  defect,  refers  to  foreign  atoms,  vacancies,  interstitials  or 
combinations  of  any  of  these.  By  measuring  the  changes  in  the 
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properties  of  a bulk  semiconductor  after  defect  introduction,  we  are 
able  to  determine  the  characteristics  of  defects.  Defects  can  be 
classified  based  on  their  physical  location  in  the  lattice  or  on  their 
electronic  structure. 

There  are  three  main  categories  of  defects.  The  first  one 
comprises  lattice  defects,  i.e.,  intrinsic  defects,  which  can  be  point 
defects,  line  defects,  or  planar  defects  [1].  The  second  category  deals 
with  impurities,  i.e.,  extrinsic  defects.  And  the  third  category  consists 
of  any  combination  of  lattice  defects  and  impurities.  There  are  two 
possible  point  defects:  the  vacancy  (Fig.  1 .(b)),  i.e.,  a vacant  atomic 
site,  and  the  interstitial  (Fig.  1 .(c)),  i.e.,  an  extra  atom  occupying  an 
interstitial  site.  An  impurity  may  replace  one  of  the  host  atoms,  in 
which  case  it  is  known  as  a substitutional  impurity.  Alternatively,  it 
may  occupy  an  interstitial  site,  in  which  case  it  may  be  either  a simple 
interstitial,  or  what  is  known  as  an  interstitialcy  [1]. 

When  a defect  is  incorporated  into  a crystal,  the  defect  can  give 
rise  to  various  localized,  charge  states  in  the  band  gap  of  a 
semiconductor.  Positively  charged  states  of  a defect  are  defined  as 
donor  states,  and  negatively  charged  states  are  defined  as  acceptor 
states.  Neutral  states  bear  no  other  distinctive  name.  Some  defects  can 
behave  either  as  donors  or  acceptors  depending  on  the  position  of  the 
local  Fermi  level.  This  characteristic  is  referred  to  as  amphoteric. 

Gold  is  widely  used  in  silicon  device  technology  to  reduce  the 
minority  carrier  lifetime.  In  spite  of  being  a subject  of  interest  and 
investigation  for  many  years,  much  has  still  to  be  learned  concerning 
this  system.  The  values  quoted  in  the  literature  for  energy  levels  and 
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Figure  1.1.  Schematic  illustration  of  various  flaws:  (a) 
perfect  crystal,  (b)  ideal  vacancy,  (c)  self-interstitial, 
(d)  vacancy-interstial  pair. 
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capture  cross  sections  show  discrepancies.  For  example,  gold  donor 
energy  values  reported  in  the  literature  vary  from  Ey  +0.33  eV  to 
Ev+0.37  eV  and  reported  hole  capture  cross  sections  range  from 
1.1  x 1CT14  cm2  to  1.5  x 10'15  cm2. 

However,  the  consensus  is  that  the  main  energy  levels  attributed 
to  gold  in  silicon  are  an  acceptor  level  at  Ec  -0.54  eV  and  a donor  level 
at  about  Ey  +0.35  eV.  Recently,  L.  S.  Lu  et  al.  [2]  drew  the  conclusion 
that  the  donor  and  the  acceptor  levels  belong  to  the  same  center;  the 
magnitude  of  the  hole  capture  cross  section  and  its  temperature 
independent  data  support  a neutral  potential  for  hole  capture  by  the 
gold  donor.  As  for  the  atomic  configuration  of  the  gold  center,  the 
published  results  seem  to  favor  the  gold-vacancy  structure  (Fig.  1 . (d)) 
of  Watkins,  in  which  the  gold  is  predicted  to  behave  like  closed  d 
shells  inside  a vacancy. 

1.2.  Noise  as  a Tool  for  Measuring  Defect  Properties 

Electronic  devices  and  systems  exhibit  random  fluctuations  in 
voltage  (or  current)  at  their  terminals,  and  these  fluctuations  are 
usually  referred  to  as  noise.  As  a result,  noise  imposes  a practical  limit 
on  the  performance  of  a circuit;  noise  can  also  be  exploited  as  a means 
for  investigating  the  electrical  characteristics  of  a circuit.  Physically, 
noise  is  due  to  the  statistical  nature  of  the  charge  transport  process;  it 
may  arise  both  from  electronic  transitions  between  discrete  energy 
levels  and  from  scattering  during  transport  [3]. 

Generation-recombination  (G-R)  noise  is  due  to  carrier 
transitions  between  traps  and  the  associated  conduction  or  valence 
band  of  a semiconductor.  Generally,  the  noise  spectrum  from  a single 
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trap  has  a Lorentzian  form  as  shown  in  Fig.  1.2.,  which  can  be  well 
described  by  a plateau  part  and  a roll-off  time  constant.  In  the  past, 
low-frequency  (L.F.)  G-R  noise  spectra  have  been  used  to  determine 
dopant  parameters  [4,5,6]  in  samples  with  impurity  concentrations 
low  compared  to  the  intended  shallow  doping  in  uniformly  doped 
material.  The  trap  energy  is  then  extracted  from  the  roll-off  time 
constants  of  the  spectral  components.  Experimentally,  G-R  noise  data 
measured  on  uncompensated  devices  can  not  be  used  to  determine 
whether  a defect  is  a donor  or  an  acceptor  center,  and  this 
information  should  be  inferred  from  the  magnitude  of  the  extracted 
carrier  capture  coefficient.  Since  noise  may  exist  in  various  forms,  for 
example:  thermal  noise,  1/f  noise,  the  proper  evaluation  of  the  G-R 
noise  magnitude  from  an  overall  noise  spectrum  is  essential  for  trap 
parameter  extraction. 

The  G-R  noise  characteristics  from  junction  devices,  such  as 
junction  field-effect  transistors  (JFET),  metal-oxide-semiconductor 
(MOS)  field-effect  transistors,  have  also  been  investigated  for  decades. 
The  activation  energy  of  the  gold  acceptor  level  of  0.55  eV  was 
extracted  from  G-R  noise  characteristics  measured  on  a JFET  by  P.  O. 
Lauritzen  [7]  and  measured  on  a MOSFET  by  L.  D.  Yau  and  C.  T.  Sah  [8]. 
The  latter  also  measured  the  thermal  emission  rates  as  a function  of 
electric  field  normal  to  the  silicon-silicon  dioxide  surface.  Recently,  F. 
Scholz  et  al.  [9]  used  L.F.  noise  data  measured  as  a function  of 
temperature  as  a MOS  device  characterization  tool  in  a manner 
analogous  to  Deep  Level  Transient  Spectroscopy  (DLTS),  and  found 
that  it  is  a powerful  diagnostic  technique  for  investigating  deep 
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Frequency  (arbitraiy  units) 

Figure  1.2.  Generation-recombination  noise  spectrum 
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defects.  In  general,  it  is  difficult  to  extract  information  about  the  field 
dependence  of  G-R  processes  from  noise  measurements  on  device 
structures.  In  order  to  do  this,  G-R  noise  from  a simple  bulk 
semiconductor  resistor  should  be  analyzed. 

In  this  dissertation,  we  address  G-R  noise  measurements  on 
samples  with  overcompensated,  deep  impurity  concentrations.  As  we 
will  indicate,  the  low- field  noise  data  combined  with  the  resistance 
data  of  the  device  can  be  used  to  determine  trap  parameters  also 
under  overcompensated  conditions.  Under  these  conditions,  the  G-R 
noise  is  dominant  over  other  noise  components  in  the  low  frequency 
range  we  investigated  because  of  the  large  impurity  concentration 
present  in  these  devices.  As  we  will  show  later,  in  an 
overcompensated  p-type  device,  it  is  difficult  to  measure  the  G-R 
noise  of  a deep  acceptor,  while  with  suitable  doping  conditions  one 
can  measure  the  defect  properties  of  a deep  donor.  Consequently,  we 
are  able  to  determine  whether  a deep  center  forms  a donor  or  an 
acceptor  center. 

1.3.  Comparison  of  Noise  Measurements  with  Other  Trap 
Characterization  Techniques 

Defects  may  be  identified  using  parameters  related  to  the 
thermodynamics  of  defect  formation,  such  as  diffusion,  or  using 
electronic  macroscopic  parameters.  The  phenomenological  electronic 
parameters,  i.e.,  thermal-emission  rates,  capture  cross  sections, 
optical  cross  sections  and  activation  energies,  can  be  measured  by 
changing  the  carrier  population  in  defects  either  thermally,  optically, 
or  electrically.  In  order  to  obtain  electronic  characteristics  of  a defect, 
one  relies  on  the  sensitivity  of  experimental  techniques  to  changes  in 
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defect  characteristics  when  subjected  to  external  disturbances.  The 
experimental  methods  for  characterization  of  electronic  properties  of 
defects  can  be  divided  into  three  groups  namely  into  electrical 
techniques,  optical  techniques  and  techniques  employing  electron 
paramagnetic  resonance  [10]. 

The  most  widely  used  electrical  techniques  are  the  transient 
methods  [11]  among  which  are  DLTS,  junction  capacitance  transient 
methods  [12],  and  pulsed  junction  capacitance  methods  [13].  Most  of 
the  above  techniques  probe  the  thermal  emission  or  thermal  capture 
of  free  carriers  by  traps  located  in  the  depletion  region  of  a diode  by 
measuring  the  transient  response  of  the  centers.  The  field  in  the 
junction  is  very  large  and  nonhomogeneous.  In  the  case  that  field- 
dependent  emission  or  capture  occurs,  the  observed  activation 
energies  will  show  some  variation  with  the  applied  electric  field.  In 
addition,  it  is  difficult  to  use  these  techniques  to  characterize  low  field 
effects  in  defects.  However,  we  can  use  transient  methods  to  identify 
acceptor  or  donor  trap  behavior  by  applying  either  a minority  or 
majority  carrier  pulse  to  junction  devices  [11]. 

Generation-recombination  noise  is  caused  by  fluctuations  in 
electron  capture  and  emission  rates.  In  G-R  noise  measurements,  the 
magnitude  of  the  homogeneous  electric  field  is  determined  from  the 
externally  applied  bias  voltage  and  the  device  length  and  the  noise 
data  are  collected  under  steady  state  conditions.  Measured  activation 
energies  are  thermal  activation  energies. 

The  response  of  deep  centers  to  these  different  experimental 
situations  and  techniques  should  help  us  to  unravel  some  of  the 
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problems  associated  with  carrier  capture  and  emission  from  deep 
centers. 

1.4.  Field  Dependent  Noise  and  Generation-Recombination  Processes 
G-R  noise  is  proportional  in  most  cases  to  the  square  of  the  d.c. 

current  flowing  through  a device,  and  it  sets  therefore  a limit  to 
device  sensitivity  under  high  field  operation.  Consequently,  it  is 
important  to  investigate  high-field  noise  spectra.  Recently,  Gasquet  et 
al.  [5]  measured  the  field-dependent  noise  of  shallow  dopants  to 
extract  trap  parameters.  However,  high  field  noise  data  from  deep 
centers  are  unavailable  at  present.  Since  G-R  noise  can  be  measured 
simultaneously  with  current- voltage  (I-V)  characteristics,  it  allows  us 
to  determine  the  effects  of  G-R  processes  on  the  I-V  curves  of  devices. 

Gold  is  a deep  center  in  silicon  and  is  widely  investigated.  Its 
high  field  capture  and  emission  coefficient  data  obtained  by 
techniques  such  as  DLTS  can  be  found  in  the  literature.  Thus  a 
comparison  of  the  results  obtained  by  G-R  noise  and  DLTS 
experiments  should  give  us  insight  in  the  dynamics  of  carriers  in  the 
presence  of  this  deep  center  and  external  disturbances.  In  addition, 
since  data  for  medium  electric  field  strengths  have  not  been  reported 
in  the  literature,  it  is  also  one  of  our  goals  to  investigate  the  exact 
field-dependence  of  emission  and  capture  coefficients  going  from  near 
equilibrium  to  the  hot  carrier  regime. 

1.5.  Microscopic  Description  of  Generation-Recombination  Processes 
The  presence  of  nonlinear  effects  in  the  high  field  operation  of  a 

device  makes  a rigorous  analytical  approach  to  any  high-field  problem 
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very  difficult.  Alternatively,  the  Monte  Carlo  method  [14]  has  proved  to 
be  a very  powerful  technique  for  the  study  of  high  field  phenomena.  In 
order  to  use  this  technique,  a microscopic  probability  formula  for  the 
trapping  and  detrapping  of  carriers  in  traps  should  be  developed. 
Although  several  macroscopic  approaches  to  G-R  problems  can  be 
found  in  the  literature  [15],  microscopic  theories  are  very  scarce  and 
microscopic  probability  formulae  for  G-R  processes  related  to  gold 
centers  in  silicon  are  unavailable.  Therefore,  an  attempt  is  being  made 
to  derive  such  formulae  using  a multiphonon  carrier  capture  model. 

1.6.  Thesis  Outline 

The  organization  of  this  thesis  is  as  follows.  In  Chapter  II,  a 
macroscopic  single-donor  G-R  noise  model  for  trap  parameter 
extraction  is  developed.  A parameter  extraction  example  is  then 
illustrated.  The  resistance  and  the  low  frequency  voltage  noise  are 
measured  as  a function  of  temperature  on  an  overcompensated  device 
made  of  Si:Au:B,  and  used  to  extract  device  parameters.  The 
normalized,  low-frequency  noise  spectral  density  and  the  resistance  of 
the  device  versus  1000/T  show  an  activation  energy  of  0.336  eV, 
which  is  in  agreement  with  the  well-established  gold  donor  energy 
level.  The  other  extracted  trap  parameters  are  also  consistent  with 
published  data. 

Chapter  III  discusses  a double-donor  G-R  noise  and  resistance 
model.  With  this  model  we  are  able  to  explain  why  different  activation 
energies  were  observed  for  the  resistance  and  the  low  frequency 
voltage  noise  measured  as  a function  of  temperature  on  a p+-p"-p+ 
device  made  of  Si:Au:B  with  a boron  concentration  of  2.7  xlO  cm'  . 
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The  normalized,  low-frequency  noise  spectral  density  versus  1000/T 
shows  an  activation  energy  of  0.336  eV,  which  is  in  agreement  with 
the  well-established  gold  donor  energy  level.  The  plot  of  the 
resistance  versus  1000/T  reveals  an  activation  energy  of  0.267  eV, 
which  we  associate  with  a second  donor  level  of  gold.  The  double- 
donor model  reduces  to  the  single-donor  model  discussed  in  Chapter 
II  under  certain  conditions. 

In  Chapter  IV,  the  current-voltage  characteristics  of  a silicon 
p+-p“-p+  device  doped  with  boron  and  gold  were  measured  and 
showed  sublinear  behavior  at  moderately  high  fields.  With  the  use  of 
low  frequency  voltage  noise  data  measured  as  a function  of  electric 
field,  we  found  that  the  sublinear  current- voltage  characteristics  are 
caused  by  a decrease  in  the  hole  concentration  with  increasing 
electric  field.  This  decrease  is  caused  by  a field-induced  decrease  in 
the  hole  emission  coefficient,  which  is  larger  than  the  simultaneously 
observed,  field-induced  decrease  in  the  hole  capture  coefficient. 

Chapter  V discusses  how  to  design  a device  to  investigate  deep 
defects  by  G-R  noise  measurement  techniques.  Under  certain 
conditions,  a multi-level  G-R  noise  model  reduces  to  a single-level  G-R 
noise  model.  We  address  the  problem  of  using  G-R  noise 
measurements  to  discriminate  between  a donor  center  and  an 
acceptor  center.  In  Chapter  VI,  a microscopic  nonradiative  carrier 
capture  model  is  proposed.  Chapter  VII  discusses  how  to  use  G-R 
noise  data  from  field-effect  transistors  to  extract  trap  parameters  and 
compares  these  techniques  with  G-R  noise  measurements  on  bulk 
devices.  In  Chapter  VIII,  a G-R  noise  model  for  polysilicon  resistors  is 
proposed.  Finally,  suggestions  for  further  research  are  made. 


CHAPTER  II 

LOW-FIELD  INVESTIGATION  OF  THE  GOLD  DONOR  LEVEL  IN 

SILICON 

2.1.  Introduction 

Gold  is  one  of  the  widely  studied  deep-level  elements  in  silicon 
due  to  its  lifetime  killer  effect.  It  is  believed  that  gold  introduces  an 
acceptor  and  a donor  level  in  the  bandgap  of  silicon.  Recently,  Lu  et  al. 
[2]  drew  the  conclusion  that  the  donor  and  the  acceptor  levels  belong 
to  the  same  center;  the  magnitude  of  the  electron  capture  cross 
section  data  supports  a neutral  potential  model  for  electron  capture  by 
the  gold  acceptor;  the  magnitude  of  the  hole  capture  cross  section  and 
its  temperature  independent  data  support  a neutral  potential  for  hole 
capture  by  the  gold  donor.  As  for  the  atomic  configuration  of  the  gold 
center,  the  published  results  seem  to  favor  the  gold-vacancy  structure 
of  Watkins  [16],  in  which  the  gold  is  predicted  to  behave  like 
closed  d shells  inside  a vacancy. 

The  experimental  methods,  by  which  the  above  information  was 
obtained,  are  generally  based  on  transient  capacitance  related 
techniques  such  as  DLTS  and  probe  the  thermal  emission  and  capture 
of  carriers  from  traps  located  in  the  depletion  region  of  a diode.  The 
emission  and  capture  coefficient  are  measured  as  two  independent 
parameters.  The  trap  activation  energy  is  obtained  from  an  Arrhenius 
plot  of  emission  coefficient  versus  temperature  assuming  that  the 
internal  electric  field  has  no  influence  on  the  results.  However, 
depending  on  the  doping  level  and  bias  voltage,  this  electric  field  can 
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be  very  large  and  field  effects  on  the  capture  as  well  as  the  emission 
processes  should  be  expected.  One  must  be  cautioned  against  wrong 
information  being  extracted  from  the  emission  data  if  the  field  effect 
is  not  properly  accounted  for.  Since  the  electric  field  will  affect  the 
carrier  velocity  distribution  and  the  trap  potential  configuration,  and 
the  field  is  not  homogeneous  and  transient  conditions  prevail,  i.e., 
transient  rather  than  steady  state  carrier  velocity  distribution  are 
generated,  it  is  difficult  to  study  and  unravel  these  effects  from 
experimental  data  measured  with  the  DLTS  technique.  To  overcome 
these  problems,  the  G-R  noise  spectra  of  a homogeneous  doped 
semiconductor  resistor  are  measured.  The  electric  field  in  such  a 
resistor  operating  in  the  space-charge  free  domain  is  independent  of 
position  and  its  value  is  determined  by  the  applied  bias  voltage  and  the 
device  length.  In  addition,  G-R  noise  measurements  are  done  under 
steady  state  conditions  so  that  transient  carrier  distributions  are 
avoided. 

In  the  following  section,  we  derive  the  expressions  for  the  G-R 
noise  spectral  density  and  device  resistance  in  terms  of  the  shallow 
and  deep  impurity  concentrations.  Then  we  present  detailed 
descriptions  of  the  experiment  and  the  parameter  extraction 
procedures. 


2.2.  Theory 

In  this  section,  we  develop  and  use  a G-R  noise  theory  to  extract 

trap  parameters  from  a defect  with  a single-donor  level.  The  device 
+ - + 

structure  is  a p -p  -p  Si:Au:B  resistor  shown  in  figure  2.1.  Figure  2.2 
is  the  schematic  energy  band  diagram  of  this  device.  The  gold  donor 
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Figure  2.1.  Device  structure 
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Figure  2.2.  Schematic  energy  band  diagram  in  equilibrium 
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energy  level  is  0.336  eV  above  the  valence  band  edge,  the  gold 
acceptor  energy  level  is  0.54  eV  below  the  conduction  band  edge.  The 
shallow  acceptor  energy  level  is  0.045  eV  above  the  valence  band  edge. 

For  simplicity,  the  following  assumptions  are  made:  (1)  only  gold 
donors  interact  with  the  valence  band,  (2)  the  electron  density  is 
neglected,  (3)  space  charge  neutrality  and  space  uniformity  prevail. 
The  carrier  densities  can  be  written  as 

P = Po  + AP,  (2.1) 

Nau+=  Nauo+  + ANau+,  (2.2) 

where  P,  Nau  are  the  hole  concentration  in  the  valence  band  and  the 

charged  gold  donor  concentration,  respectively:  AP,  ANau+  are  the 

corresponding  fluctuations  in  P and  Nau  and  PG  and  Nauo+  are  the 

corresponding  densities  in  steady  state. 

The  condition  of  charge  neutrality  gives  the  following  equation 

for  the  fluctuations  in  P,  and  Nau+ 

ANau++AP  = 0.  (2.3) 

The  rate  equation  takes  the  form 

dP  + + 

■jjp  = ep  Nau  - cp  P (Nau  - Nau  ),  (2.4) 


where  cp  is  the  hole  capture  coefficient  and  ep  is  the  hole  emission 
coefficient. 

In  steady  state,  eq.  (2.4)  gives 

p 

ep  = cp  (Nau  -Nauo+  )( ° +),  (2.5) 

Nau 

Using  eqs.  (2.1)  through  (2.5),  we  get  after  linearization 
dAP 

= - [ep  + cp  P0  +cp  (Nau-  Nauo  )]AP. 


(2.6) 
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This  equation  can  be  further  simplified  as  we  will  show  in  the 
following  section. 


2.2.1.  Overcomnensated  Devices 

If  the  gold  concentration  is  larger  than  the  boron  concentration, 
i.e.,  the  device  is  overcompensated,  the  Fermi  level  will  be  located 
between  the  gold  donor  energy  level  and  the  valence  band  edge  at 
temperatures  between  150  and  250  K.  Then  P is  very  small,  more 
precisely  P « (Nau-  Nau+),  Nau+,  Nau. 

With  the  aid  of  eq.  (2.5),  we  can  assume  that 

[ep+cpP0]«  [cp  (Nau- Nauo+)].  (2.7) 

Then  eq.  (2.6)  simplifies  to 


dAP  , . AP 

“cj-j—  = - cp  (Nau-  Nauo  ) AP=  - — — , 


with 


t— = cp  (Nau  - Nauo+). 

Lp 


(2.8) 


(2.9) 


Assuming  Na-Na  -Nauo  , where  Na  is  the  concentration  of  boron,  eq. 
(2.9)  becomes 


1 


t— = cp  (Nau  - Na). 

tp 


(2.10) 


The  normalized,  spectral  density  of  G-R  noise  has  the  form  [17] 

2 


Sw((0)  SDD(C0) 


Vo 


'PP*W' 4 go  tp 

PoZ  L w d P02  (1+CO2  tp2) 


(2.11a) 


where  gG  is  the  free  carrier  generation  rate  given  by 

g0=  Nauo+  ep.  (2.11b) 

The  normalized  plateau  value  of  the  G-R  noise  can  then  be  written  as 
Sw«»  4 Nauo+  ep  tp2 


V' 


72 — = 


LwdP0^ 


(2.12) 
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The  density  Nauo+  can  be  expressed  in  terms  of  Nau  using 

Fermi-Dirac  statistics 

xt  + Nau 

Nauo  = v — . (2.13) 


E— E . 

i au 

l+gaue  KT 

where  Eau  is  the  gold  donor  energy  level,  Ef  is  the  Fermi  level,  K is  the 
Boltzmann  constant,  T is  the  temperature,  and  gau  is  the  degeneracy 
factor  of  the  gold  donor. 

Combining  eqs.  (2.5)  and  (2.13),  we  get 

e - E 
v au 

KT 

ep  = Cp  Ny  gau  e . (2.14) 

Using  the  charge  neutrality  equation 

Na  = P + Nau+  , (2.15a) 

and  assuming  that  the  shallow  acceptors  are  fully  ionized,  P can  be 
written  as 


Ev~  Eau 


P = 


gau  (Na  - P ) Nv  e 


KT 


(2.15b) 


Nau-  Na  + P 

Using  eqs.  (2.10),  (2.12),  (2.14)  and  (2.15b),  and  assuming  that 
the  hole  density  is  small,  we  get 


Sw(°) 

V 2 


4e 


E -E 
au  v 

KT 


LwdNvcpNagau  ’ 


(2.16) 


where  L,  w,  d are  the  length,  width,  and  thickness  of  the  device,  Ev  is 
the  valence  band  edge,  and  Nv  is  the  density  of  states  of  the  valence 
band. 

The  resistance  of  the  device  is  given  by 
L 


R = 


e w d 


(2.17) 
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Using  eqs.  (2.15b)  and  (2.17),  we  get 
L(N_ 


R = 


'au  ^a^  ^ 


E - E 
au  v 

'■  KT 


’au 


ewd  MpNaNv  g 

where  Jtp  is  the  hole  mobility,  and  -e  is  the  charge  of  an  electron. 


(2.18) 


Equations  (2.16)  and  (2.18)  show  that  the  activation  energies 
extracted  from  these  equations  should  have  the  same  value. 

With  the  use  of  a low  field  mobility  formula,  eqs.  (2.16)  and 
(2.18),  and  an  optimization  procedure  [18],  we  are  able  to  extract  trap 
parameters  from  the  previous  equations.  The  semi- empirical  mobility 
formula  that  we  use  is  [19] 


ftp-54. 3( 


-2.23 

T 0’57  1.36  xlO8  (T) 

300 ] + 17  t 2 4 T 0146  ’ 

l+[Nf/(2.35  x1017(3^q)  )]0.88  (^qq) 

(2.19) 


where  Nf  is  the  effective  charged  impurity  concentration,  and  all 
quantities  are  in  c.g.s.  units. 


2.2.2.  Lightly  Doped  Devices 

If  the  gold  concentration  is  smaller  than  the  boron 
concentration,  we  may  assume  that 

P » Nau+,  Nau.  (2.20) 

With  this  condition  eq.  (2.6)  simplifies  to 


dAP  r AP 

— " [Cp  + Cp  PQ  ]AP=  - , 

with 

^ = [Cp  + Cp  P0  ]. 

The  concentration  of  holes  is  given  by 


(2.21) 


(2.22) 
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E -E, 

p0  = nv  e~^~. 

Combining  eqs.  (2.14),  (2.22)  and  (2.23),  we  get 

Ef  ' Eau 

+ 1 ]• 


(2.23) 
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(2.24) 


When  the  Fermi  level  is  located  above  the  trap  level,  the  above 
equation  becomes 


E 


au 


tp  ” CP  Nv  ^ 


KT 


(2.25) 


The  trap  energy  can  be  extracted  from  this  equation  by  plotting 
Ln(tpT2)  versus  1000/T  where  it  is  assumed  that  the  temperature 
dependence  of  cp  is  weak. 

From  eq.  (2.11)  we  get 

Spp(O)  4 Nau+  en  tr 


Sw(°) 


tpv- 


T = 


tnP2 


L w 


(2.26) 


Combining  eqs.  (2.13),  (2.14),  (2.17),  (2.23)  and  (2.24),  we  get 
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(2.27) 


(2.28) 


Equation  (2.27)  shows  a peak  for  Y =0  where  T=T0.  From  this  peak 
value  the  gold  concentration  can  be  determined.  From  eqs.  (2.23)  and 
(2.28),  we  get 
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Wheny>0,  combining  eqs.  (2.23),  (2.27)  and  (2.28),  we  get 
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When  y<  0,  combining  eqs.  (2.23),  (2.27)  and  (2.28),  we  get 
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(2.30) 
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(2.31) 


The  trap  parameters  can  be  extracted  by  using  the  following 
procedures  as  outlined  in  [20] 

(1)  The  trap  energy  is  extracted  from  eq.  (2.30)  or  (2.31). 

(2)  The  degeneracy  factor  is  extracted  from  eq.  (2.29). 

(3)  The  gold  concentration  is  extracted  from  eq.  (2.27). 

(4)  The  capture  coefficient  is  extracted  from  eqs.  (2.23)  and  (2.24). 

(5)  The  emission  coefficient  is  extracted  from  eq.  (2.14). 


2.2.3.  Arbitrary  Doping  Levels 

If  the  impurity  concentration  does  not  satisfy  eq.  (2.7)  or  eq. 
(2.20),  the  full  expression  for  eq.  (2.6)  has  to  be  considered;  i.e.. 


dAP  . AP 

= - [ep  + cp  P0  +cp  (Nau-  Nauo  )]AP=  - — — , 


with 


— = [ep+  cp  P0  +cp  (Nau-  Nauo+)]. 


Substituting  for  ep,  PQ  and  Nauo+,  we  get 


(2.32) 


(2.33) 
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Equation  (2.12)  then  becomes 
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(2.35) 
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Combining  eqs.  (2.17),  (2.19),  (2.23),  (2.34)  and  (2.35),  and  using  an 
optimization  technique  [18],  we  can  get  the  degeneracy  factor,  gold 
concentration  and  trap  energy  level  from  eq.  (2.35),  and  then  from  eq. 
(2.34)  the  capture  coefficient.  The  emission  coefficient  is  extracted 
from  eq.  (2.14). 


2.3.  Experimental  Results  and  Discussion 

For  devices  with  small  deep  impurity  concentrations,  the 

method  of  section  2.2.2  applies,  as  is  experimentally  demonstrated  in 

[20].  The  experimental  data  obtained  from  our  devices  show  that  we 

deal  with  overcompensation.  Therefore  the  theory  of  section  2.2.1 

should  apply. 

+ - + 

A p -p  -p  resistor  (CM24B-82  #1)  shown  in  figure  2.1  is  used 
in  this  experiment.  An  epi-layer  of  thickness  20.1  um  is  grown  on  a 
highly  doped  substrate  with  a resistivity  of  0.02  Ohm-cm.  The  boron 
concentration  in  the  epi-layer  is  1.0  x 1015  cm'^.  The  drive-in 

O 

temperature  is  900  Celsius;  the  drive-in  time  is  10.75  hours.  The 
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junction  diameter  of  the  device  is  7.5  x 10~2  cm  and  the  calculated 
junction  depth  of  the  Ohmic  contact  is  1.93  um.  The  resistance  of  the 
device  at  room  temperature  is  64  Ohm.  A guard  ring  to  control  the 
surface  properties  of  the  p+-p  -p+  device  was  included  in  the  device 
design. 

When  the  gold  concentration  is  larger  than  the  boron 
concentration  and  the  shallow  impurities  are  completely  ionized,  the 
gold  impurities  will  be  partially  ionized  and  the  Fermi  level  will  be 
located  between  the  gold  donor  energy  level  and  the  valance  band 
edge.  As  it  will  be  shown  in  Chapter  V,  the  Fermi  level  is  close  to  and 
significantly  affected  by  the  gold  energy  level,  and  varies  linearly  with 
temperature. 

A block  scheme  of  the  experimental  arrangement  is  shown  in 
figure  2.3.  The  voltage  Vs  is  the  bias  voltage,  Rd  represents  the  device 
and  Rb  the  bias  resistance.  The  device  is  mounted  in  a temperature 
controlled  holder  which  allows  for  noise  measurements  between  77 
and  300  K.  The  bias  current  I,  supplied  from  a battery,  is  controlled  by 
a variable  series  metal  film  resistor  which  is  almost  20  times  larger 
than  the  resistance  of  the  device.  A Brookdeal  5004  ultra-low  noise 
preamplifier  and  a Hewlett  Packard  3561 A dynamic  signal  analyzer 
make  it  possible  to  measure  the  device  noise  at  frequencies  between  1 
and  105  Hz. 

A switch  is  used  to  connect  the  input  of  the  amplifier  to  a 
dummy  resistor,  by  which  we  are  able  to  estimate  the  magnitude  of 
the  undesirable  noise  from  the  amplifier  and  the  bias  circuit.  We  found 
that  the  noise  from  the  bias  circuit  is  negligible  when  compared  with 
the  generation-recombination  noise  from  the  device. 
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Figure.  2.3.  Block  scheme  of  the  experimental  set-up 
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The  measured  noise  spectra  show  frequency  independent 
excess  noise  up  to  frequencies  where  the  bandwith  of  the  measuring 
systems  limits  the  signal  transfer.  The  excess  noise  is  interpreted  in 
terms  of  G-R  noise.  Biasing  of  the  guard  ring  has  no  effect  on  the 
measured  spectral  density.  This  indicates  that  the  noise  is  generated 
in  the  bulk  of  the  resistor. 

Figure  2.4  shows  the  resistance  of  the  device  as  a function  of 
1000/T.  Figure  2.5  shows  the  measured  values  of  the  G-R  noise  plateau 
normalized  to  V2  as  a function  of  1000/T.  With  the  use  of  the  semi- 
empirical  formula  (2.19)  for  the  mobility  and  eqs.  (2.16)  and  (2.18), 
we  extract  a hole  capture  cross  section  of  <*pi  = 3.2  xlO’15  cm2,  a 
degeneracy  factor  of  16.25,  a gold  donor  energy  above  the  valence 
band  of  0.336  eV,  and  a gold  concentration  of  2.51  xlO15  cm’3. 

Note  that  the  extracted  value  of  the  gold  degeneracy  factor  is 
quite  close  to  16  as  found  from  Hall  coefficient  data  [21],  and  that  the 
value  of  the  hole  capture  cross  section  is  of  the  same  order  as 
published  values  [2,22],  This  indicates  that  the  above  described  noise 
method  is  an  accurate  technique  to  investigate  deep  level  properties. 


Resistance  (Ohm) 
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Figure  2.4.  The  resistance  versus  lOOO/T  of  device  CM24B-82  #1. 
The  squares  are  the  experimental  data  and  the  solid  line 
represents  the  best  fit  to  the  data. 


(Ol/lv^s) 
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Figure  2.5.  The  normalized,  low-frequency  G-R  noise  plateau  values 
versus  1000/T.  The  circles  represent  the  experimental  data  and  the 
solid  line  represents  the  best  fit  to  the  data. 


CHAPTER  III 

EXPERIMENTAL  EVIDENCE  FOR  A SECOND-DONOR  LEVEL  OF  GOLD 

IN  SILICON 

3.1.  Introduction 

The  complete  signature  of  an  impurity  in  a given  semiconductor 
host  should  include  information  about  its  atomic  configuration  and 
wavefunction  and  the  energy  levels  of  bound  electrons  and/or  holes. 
From  these  quantities  it  is  possible  to  deduce  all  its  electronic 
characteristics.  The  phenomenological  electronic  parameters,  i.e., 
thermal-emission  rates,  capture  cross  sections  and  optical  cross 
sections,  can  be  measured  by  changing  the  carrier  population  in 
defects  either  thermally,  optically,  or  electrically.  In  order  to  obtain 
the  defect's  electronic  characteristics  one  relies  on  the  sensitivity  of 
the  experimental  techniques  to  the  changes  in  defect  characteristics 
when  subjected  to  external  disturbances.  The  experimental  methods 
for  characterization  of  electronic  properties  of  defects  can  be  divided 
into  three  groups  namely,  electrical  techniques,  optical  techniques 
and  techniques  involving  electron  paramagnetic  resonance  [10].  The 
most  widely  used  electrical  techniques  are  the  transient  methods  [11]. 

These  transient  measurement  techniques  are  rather  involved 
and  mainly  based  on  the  assumption  that  there  is  only  one  dominant 
emission  or  capture  process  of  either  electrons  or  holes  for  each  deep 
level  during  the  transient.  The  effect  from  a small  or  slow  response  is 
not  easily  observed,  simply  because  it  has  a minor  contribution  to  the 
overall  signal.  For  example,  if  two  energy  levels,  say  Ei,  E2,  are  very 
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close  and  the  emission  rate  from  level  Ei  is  smaller  than  that  of  level 
E2,  then  we  are  unable  to  detect  energy  level  Ei  by  measuring 
emission  rates  versus  1000/T. 

The  activation  energy  for  emission  is  obtained  in  these 
experiments  by  relating  the  emission  rate  to  the  energy  level  depth  Ea 
by  detailed  balance.  If  the  emission  coefficient  depends  on  electric 
field,  then  the  activation  energy  extracted  from  these  methods  will 
show  a variation  with  electric  field  strength. 

In  chapter  II,  we  used  resistivity  and  noise  data  to  explore  the 
defect  properties  of  gold  in  silicon  and  found  that  the  results  are 
compatible  with  those  obtained  by  other  techniques.  The  variation  of 
resistance  versus  temperature  provides  under  certain  conditions  the 
activation  energy,  and  the  noise  data  give  us  the  capture  cross  section 
and  emission  coefficient.  In  this  noise  experiment  the  magnitude  of 
the  electric  field  was  controlled  precisely  by  the  externally  applied 
voltage. 

Despite  the  fact  that  gold  in  silicon  has  received  so  much 
attention,  its  behavior  is  not  always  predictable.  The  consensus  seems 
to  be  that  a donor  and  an  acceptor  level  at  0.35  and  0.54  eV  above  the 
valence  and  below  the  conduction  band  respectively  are  introduced  by 
gold  in  the  silicon  energy  gap.  This  is  the  so-called  two-level  model. 
There  are  some  indications,  however,  that  gold  can  produce  more 
than  two  levels.  Newman  [23]  measured  the  photoconductive  spectral 
response  of  p-type  gold-doped  silicon  and  pointed  out  that  it  is 
difficult  to  explain  his  data  in  terms  of  a two-level  model. 

A search  for  another  acceptor  level  closer  to  the  conduction 
band  was  conducted  [24].  This  was  done  by  changing  the 
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concentration  ratio  of  gold  (Nau)  to  residual  donors  (Nd-Na)  in  a single 
crystal  by  several  additions  of  phosphorus  from  the  case  where 
Nau  > (Nd-Na)  to  the  case  where  2 Nau  > (Nd-Na)  > Nau.  No  new  level 
was  discovered.  A similar  experiment  in  p-type  silicon  was  carried  out 
by  Taft  and  Horn  [25].  They  doped  p type  silicon  devices  with  different 
gold  concentrations.  No  exact  relationship  between  the  gold  and  the 
shallow  impurity  concentration  in  the  p-type  silicon  was  specified,  and 
no  new  energy  level  was  found. 

In  this  chapter  we  report  on  resistance  and  low  frequency 
voltage  noise  measurements  carried  out  as  a function  of  temperature 
on  a p+-p--p+  device  made  of  Si:Au  with  a boron  concentration 
of  2.7  x 10  cm’  . Two  different  activation  energies  are  observed.  The 
normalized,  low-frequency  noise  spectral  density  versus  1000/T 
shows  an  activation  energy  of  0.336  eV,  which  is  in  agreement  with 
the  well-established  gold  donor  energy  level.  The  plot  of  the 
resistance  versus  1000/T  reveals  an  activation  energy  of  0.267  eV.  It  is 
very  unlikely  that  this  energy  level  is  caused  by  an  unintended 
impurity  since  the  devices  were  fabricated  under  very  well  controlled 
conditions.  Therefore  a double-donor  model  is  proposed  to  interpret 
this  behavior. 


3.2.  Theory 

In  this  section,  we  discuss  a method  of  using  low-frequency 
generation-recombination  noise  involving  a double  donor  as  a tool  for 
trap  parameter  extraction.  The  device  is  a p+-p  -p+  Si:Au:B  resistor 
shown  in  figure  3.1.  Figure  3.2  is  the  schematic  energy  band  diagram 
of  the  device.  As  we  will  show  in  this  chapter,  gold  in  silicon  creates 
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Figure  3.1.  Device  structure 
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Figure  3.2.  Schematic  energy  band  diagram  in  equilibrium 
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two  donor  levels:  the  first  donor  level  lies  0.336  eV  above  the  valence 
band  edge  and  the  second  donor  level  0.267  eV  above  the  valence 
band  edge.  The  gold  acceptor  energy  level  is  0.54  eV  below  the 
conduction  band  edge.  If  the  Fermi  level  is  far  away  from  the  acceptor 
energy,  the  effects  from  the  acceptor  level  can  be  neglected  (see 
Chapter  V).  The  boron  activation  energy  is  0.045  eV. 

In  the  previous  chapter,  we  discussed  a single-donor  model  in 
which  the  Fermi  level  is  located  above  the  second  donor  level. 
Consequently,  the  gold  donor  centers  in  silicon  are  almost  singly 
ionized,  and  the  effects  of  the  second  donor  level  on  the  electronic 
properties  can  be  neglected.  In  this  chapter,  we  consider  a double- 
donor model,  that  under  certain  conditions  may  reduce  to  a single- 
donor model. 

The  calculations  of  the  resistance  and  the  voltage  noise  resulting 
from  double  donor  centers  interacting  with  the  valence  band  are  made 
based  on  the  following  assumptions:  1.  space  uniformity  prevails:  2. 
only  donor  levels  interact  with  the  valence  band,  i.e.,  we  neglect  the 
interaction  between  the  gold  acceptor  level  at  midgap  and  the  valence 
band;  3.  P > N;  4.  Na  = Na. 

The  pertinent  hole  transitions  are  described  in  Table  3.1.  The 
parameter  P is  the  hole  density,  N is  the  electron  concentration,  Na  is 
the  boron  concentration,  and  Na  is  the  ionized  boron  concentration. 
The  total  gold  concentration  is  Nau,  and  Nau+,  Nau++,  Nau°  are  the 
singly  ionized,  doubly  ionized,  and  neutral  gold  concentration 
respectively.  The  symbol  e+  denotes  a mobile  hole,  gi,  ri  are  the  hole 
generation  rate  and  the  hole  recombination  rate  respectively,  and  epi, 
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TABLE  3.1.  Generation  and  recombination  rates  in  terms  of 
ionized  gold  and  mobile  hole  concentrations. 


Rate 

Hole  Transitions  involved 

Description  of  the  rates 

gi 

+ o , 

Nau  ->  Nau  + e 

epi  Nau+ 

° + + 
Nau  + e ->  Nau 

cpi  Nau°  P 

& 

++  + , 
Nau  ->  Nau  + e 

eP2  Nau++ 

r2 

+ _L  + + 

Nau  + e ->  Nau 

cp2  Nau+  P 
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cpi,  and  eP2,  cP2  are  the  emission  and  capture  coefficient  for  holes  in 
the  first-donor  level  and  the  second-donor  level  respectively. 

According  to  W.  Shockley  and  J.  T.  Last  [26],  the  ionized  impurity 
densities  can  be  expressed  as 
1 


Elcf  Ef 


Nau+  = 


Nau++  = 


Saul 


o cx  K T 

Nau  e 


Sau2 


Nau0  e 


Eld+  E2d'  2 Ef 
ITT 


(3.1) 


(3.2) 


where  Eld  is  the  first  gold  donor  energy  level,  E2d  is  the  second  gold 
donor  energy  level,  Ef  is  the  Fermi  energy  level,  and  gaUl  and  gau2  are 
the  degeneracy  factors  for  the  first  and  second  donor  level 
respectively. 

The  total  gold  concentration  is  given  by 


Nau++  + Nau+  + Nau°  = Nau. 


(3.3) 


From  eqs.  (3.1),  (3.2)  and  (3.3),  the  ionized  gold  concentrations  are 
rewritten  as 

Nau+  = 


Nau 


E2d  Ef 


- . ’ 
^aul  „ KT  KT 

1+~ e +Saule 


(3.4) 


Sau2 


Nau++  = 


Nau 


,Eid+E2d'2  Ef*  ""  - (E2d-  Ef^ 

+ «au2e 


(3.5) 


1 + ^au2  ^ 


KT 


KT 


^aul 


The  value  of  P can  be  derived  from  eqs.  (3.3)  to  (3.5)  and  eq.  (2.23) 
with  the  use  of  the  charge  neutrality  equation 

2Nau++  + Nau+  + P = Na.  (3.6) 

Our  noise  calculation  applies  to  the  low  field  Ohmic  regime  of 
device  operation.  Then 
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I = ^=r^P  pv-  (3-7) 

where  |4p  is  the  hole  mobility,  V is  the  applied  voltage,  and  R is  the 
resistance. 

The  rate  equations  for  the  hole  and  ionized  gold  densities  read: 
dNau+ 

at  - • gi  + ri  + g2  - r2,  (3.8) 


dNau++ 

dt  - ' g2  + r2, 
dP 

dt  = gi  - ri  + g2  - r2. 


(3.9) 

(3.10) 


where  gi,  ri,  g2,  r2  are  described  in  Table  3.1. 

In  order  to  calculate  the  noise  spectral  density,  we  express  the  carrier 
densities  as  sums  of  a steady-state  term  and  a fluctuation  term, 

P = P0  + AP,  (3.11) 


Nau+  = Nauo+  + Nau+, 


Nau++  = Nauo++  + Nau++, 


++ 


.++ 


(3.12) 


(3.13) 


+ + 


where  P,  Nau  , Nau  are  the  hole  concentration  in  the  valence  band, 
the  singly  ionized  gold  concentration  and  the  doubly  ionized  gold 
concentration  respectively.  The  terms  AP,  ANau+,  ANau++  are  the 
corresponding  fluctuations  in  the  densities  of  hole  and  charged  traps 
and  P0,  Nauo  , Nauo  are  the  corresponding  densities  in  steady-state. 
From  eqs.  (3.3)  and  (3.6),  we  obtain 


ANau++  + ANau+  + ANau°  = 0, 

and 


(3.14) 


2ANau++  + ANau+  + AP  = 0. 


(3.15) 
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We  expand  the  carrier  densities  around  steady-state  in  first  order  and 
obtain  for  eq.  (3.10) 


dt  - gio  _ no  + §20  ■ r20  + epi  ANau+  - cpi  Nau°  AP 

- cpi  ANau°  P0  + eP2  ANau++-cP2  ANau+  P0-cP2  Nau+  AP,  (3.16) 
where  gio,  ri0,  g2o.  f2o  are  the  steady-state  values  of  gi,  ri,  g2,  r2, 
respectively. 

Using  eqs.  (3.14),  (3.15),  (3.16)  and  the  steady-state  solution  of  eqs. 
(3.9)  and  (3.10),  we  get 


dAP 

dt 


-(epi-cP2  Po+Cpi  Nauo°+cP2  Nauo++cpi  P0)AP- 


(_6p2"i"Cpi  Po+2epi_2cp2Po)  ANau  , (3.17) 

dANau++ 

— ^ = - (eP2  + 2cp2  P0)  ANau  + (cp2  Nauo  -cp2  P0)AP.  (3.18) 


Using  the  Langevin  approach,  the  rate  equations  including  the 
noise  can  then  be  written  as 


dAP 

dt 


= -MnAP-Mi2  ANau++  + Fi  + F2, 


(3.19) 


dANau++ 

— ^ = -M2i  AP  - M22  ANau  - F2, 

with  relaxation  matrix  elements 


(3.20) 


Mn  = epi-cp2  Po+Cpi  Nauo°+cP2  Nauo++cpi  Po, 

M22  = eP2  + 2cP2  P0, 

M12  = _ep2+CpiP0+2epi-2cp2Po, 

M2i  =-  cp2  Nauo+  + cp2  P0.  (3.21) 

The  Langevin  noise  source  term  Fi  models  the  noise  in  the  hole 
transition  rates  between  the  valence  band  and  the  first-donor  level, 
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whereas  F2  accounts  for  the  noise  in  the  hole  transition  rates  between 
the  valence  band  and  the  second-donor  level. 

Performing  a Fourier  transformation  on  eqs.  (3.19)  and  (3.20), 
we  obtain 


TjCO+Mn 

M12 

AP 

F1+F2 

L M21 

j(0+M22  _ 

_A~Nau++_ 

- -F2  _ 

Inverting  this  expression  results  in 


~ _ (](0+M22  ) (Fi  +F 2)  +M 1 2F2 
Det 


(3.23) 


with 

Det  = -(02+jC0(Mn+M22)"Mi2M2i  +M22M11,  (3.24) 

where  AP,  Fi,  F2  are  the  Fourier  components  of  AP,  Fi  and  F2, 
respectively. 

The  eigenvalues  A-i , X2  of  the  phenomenological  relaxation  matrix  M 
[3]  are  the  roll-off  frequencies  of  the  noise  spectra.  The  eigenvalue 
equation  of  M reads: 

X2  — A,(Mii+M22)"(Mi2M2l-  MuM22)=0. 

The  eigenvaules  are  therefore  given  by 

^ (Mh+M22)+V  (Mii+M22)2"4(MhM22_Mi2M2i)  , 

Ai  = 2 * (3.25) 

'i  (Mh+M22)W  (Mi1+M22)2-4(M1iM22-Mi2M2i) 

A2  = 2 • (3.26) 


Assuming  that  there  is  no  correlation  between  the  Langevin  noise 
sources  and  that  the  transition  rates  show  full  shot  noise,  we  obtain 


(«)= 


(C02+M222  )(Sf1  +Sf2)  + (Mi22+2Mi2M22)  Sf2 

I Det  I 2 


with 


Sfi  = 


4ri 

Lwd  ’ 


(3.27) 
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and 


Sp2= 


4r2 


Lwd  ’ 

where  L,  w,  d are  the  length,  width  and  thickness  of  the  device. 
Thus 

G H 

SPP(<0,=  (C02+V)  + (®W)  • 


(3.28) 


with 


G = 4' 


^l2(rl+r2)~M222  (rl+r2)-(Mi22+2Mi2M22)  r2 

Lwd  (X,i2-^22) 


(3.29) 


and 


tt  „ _^22(ri+r2)+M222  (ri+r2)+(Mi22+2Mi2M22)  r2 

H = 4 7 — 777TT1^ • (3-30) 

Lwd  (ki  -hi ) 

Note  that  eq.  (3.28)  predicts  two  Lorentzian  components  for  the 
overall  noise  spectrum.  Equations  (3.28)  through  (3.30)  can  be  further 
simplified. 


3.2.1.  Deep  Double-Donor  Defects  in  Nearly  Compensated  Material 

In  this  section,  a device  doped  with  a gold  concentration 
satisfying  Na  > Nau  > 0.5Na  is  studied.  From  the  steady-state  solution 
of  eqs.  (3.8)  and  (3.9),  we  find 


p 

ePi  = cpi  Nauo0  ( -'T_°  _+), 


and 


eP2  = cp2  Nauo+  ( 


Nauo 


Po 


a. 


(3.31) 


(3.32) 


+ + 


Nauo++ 

At  low  temperatures  we  may  assume  P0«  Nauo",  Nauo  ",  Nauo".  This 
assumption  is  substantiated  by  exact  numerical  calculations  presented 
in  section  3.3.  We  make  this  assumption  at  this  point  because  it  allows 
us  to  obtain  simple  analytical  expression  for  the  resistance  and  voltage 
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noise  of  the  double  donor  model.  With  eqs.  (3.31)  and  (3.32)  and  the 
above  assumption  eqs.  (3.17)  and  (3.18)  can  be  simplified  to 


dAP 

dt 


=-(cpi  Nau°+cp2 


Nau+)AP-(-eP2+CpiP+2epi-2Cp2P)ANau++ , 


dANau++ 

dt 


= (cp2  Nau+)AP  - (ep2  + 2cp2  P)  ANau++  . 


(3.33) 

(3.34) 


The  relaxation  matrix  elements  in  eq.  (3.21)  become 


Mn=  cpi  Nau°+cp2  Nau+,  M12  = - eP2+CpiP+2epi-2cp2P, 


M2i=  - cp2  Nau+,  M22  = ep2  + 2cp2  P.  (3.35) 

Note  that  eq.  (3.28)  predicts  two  Lorentzian  components  in  the  overall 
noise  spectrum.  However  in  the  following  we  argue  that  the  values  of 
^*1  and  X,2  are  such  that  in  our  experiments  where  we  measured  the 
noise  between  100  Hz  and  1000  Hz,  only  SPp(C0)  for  'k2«  CO  « A,i  was 
observed. 

Since  the  coefficient  cP2  is  associated  with  a repulsive  capture 


process,  one  may  expect  its  value  to  be  rather  small  [27],  of  the  order 
of  10" 11  to  10" 14  cm3/s.  The  coefficient  cPi  on  the  other  hand  i 


is 


related  to  a neutral  capture  process,  and  will  be  therefore  of  the  order 
of  10"8  cm3/s.  Consequently  we  may  assume 


cpi  Nau+ 
» 


Eld"  Ef 
KT 


Cp2  Nau  1 

Using  eqs.  (3.25),  (3.26)  and  (3.31),  (3.32),  (3.36),  we  obtain 


(3.36) 


Mu  » M21,  M22.  M12, 


>>  X2 , 

^1  = Ml  1 =Cpi  Nau°  = T— , 

Lp 


(3.37) 
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~ (M11+M22)  < M22  _ep2  + 2Cp2  P' 

Our  exact  numerical  calculations  indicate  an  upper  limit  of  1 Hz  for 
A,2/2x  reached  at  a temperature  of  221  K.  This  is  consistent  with  the 
fact  that  we  did  not  observe  any  Lorentzian  spectrum  in  our  low 
frequency  noise  window  of  100  Hz  to  1000  Hz. 

Using  eqs.  (3.31),  (3.32)  and  (3.37),  we  get 
4(ri+r2) 


^pp^  L w d (CO2  +^i2)  ’ 


(3.38) 


ri+  r2=  (cpi  Nau  +cp2  Nau  )P  = MnP, 


Sw(fO) 

V2 


Spp(CO)  4 cpi  Nau°  P tp2 
P2  'LwdP2  (1+co2  tD2) 


(3.39) 


Equation  (3.39)  can  be  approximated  by 


4 . 

“ 772 — = : TT; — o for  A2«  0)  « Ai.  (3.40) 

V L w d P cpi  Nau 

This  equation  can  be  simplified  to  eq.  (2.16)  if  we  use  the  hole 
concentration  given  by  eq.  (2.15).  This  shows  that  when  the  Fermi 
level  is  above  the  first  donor  energy  level,  the  double-donor  model 
reduces  to  a single-donor  model.  Equations  (3.6),  (3.7)  and  (3.40)  can 
be  used  to  extract  trap  parameters  from  the  resistance  data  and  noise 
data  if  we  use  the  mobility  formula  (2.19). 

The  theoretical  predictions  for  the  noise  and  resistance  can  be 
further  simplified  by  assuming  (Eld  - Ef)  > KT  log(gaui)  which  is 
certainly  valid  for  the  low  temperature  region  of  our  experiment.  Then 
eq.  (3.5)  can  be  approximated  by 


Nau++  = 


Nau 


1 + 


&au2 

^aul 


~7e 


2d~N 

KT 


(3.41) 
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Using  eqs.  (3.3),  (3.6),  neglecting  P in  eq.  (3.6)  and  assuming 
Nau++,  Nau+  > Nau°,  the  hole  concentration  at  low  temperature  is  given 


E - E0  , 
v 2d 


P = 


&au2  ,»T  »T  , „T  ^ kt 
— (Na  - Nau  ) N e 

&aul 

2 Nau-N„ 


Using  the  following  approximations. 


(3.42) 


Nau++  = Na-Nau,  (3.43) 

Nau+  = 2Nau-Na,  (3.44) 

and  substituting  eqs.  (3.1)  and  (3.42)  into  eq.  (3.40),  the  expression 
for  the  normalized  G-R  noise  becomes 


Sv>) 

7I~ 


4 e 


Eld-  Ev 
KT 


for  A*2«  CO  « A,i.  (3.45) 


Vz  L w d cpi  gau  l (2Nau-Na  ) Ny 

Equations  (3.42)  and  (3.45)  illustrate  that  the  noise  contains  the 
activation  energy  of  the  first  donor  level,  whereas  the  carrier 
concentration  or  resistance  shows  the  activation  energy  associated 
with  the  second  donor  level. 


3.2.2.  Simplified  Double-Donor  Noise  Theory 

Section  3.2.1  dealt  with  the  double-donor  model  rigorously. 
However,  the  same  results  can  be  obtained  using  a simple  idea. 
Suppose  that  carrier  fluctuations  can  be  caused  by  two  processes:  (1) 
carrier  transitions  between  the  first-donor  level  and  the  valence  band 
such  that  AP  = -ANau+,  ANau++=  0,  (2)  carrier  transitions  between  the 
second-donor  level  and  the  valence  band  such  that  APsANau+=-ANau++. 
We  therefore  approximate  the  rate  equation  as 
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For  case  1, 


dAP  . 

= epi  ANau  - cpi  Nau  AP  - cP2  Nau  AP-  cP2  ANau  P, 

= - (epi  + cpi  Nau0  + cP2  Nau+  - cP2  P)AP.  (3.46) 

Using  eqs.  (3.31),  (3.32)  and  the  conditions  P«  Nau++,  Nau+ at  low 
temperature,  we  get 
dAP 

= - ( cpi  Nau  + cp2  Nau  )AP.  (3.47) 

For  case  2 we  obtain 
dAP  + 

-^j-=  epi  ANau  -cpi  Nau  AP-  cP2  Nau  AP 

- cp2  ANau+  P+eP2  ANau++, 

= - (-ePi  + cpi  Nau0  + cp2  Nau+  + cp2  P - ep2)AP 

-AP 

= - ( cpi  Nau°  + cp2  Nau+  )AP  = -r — . (3.48) 

tp 

Equations  (3.47)  and  (3.48)  show  the  same  life  time  constant  in 
both  cases.  Since  the  carrier  fluctuations  in  the  valence  band  arise 
from  the  processes  listed  in  Table  3.1,  the  following  spectral  density 
expression  for  G-R  noise  is  valid  [17] 

2 


Sw<“) 


V' 


2 - 


4gV 


with 


L w d P2  (1+co2  tp2) 


(3.49) 


g =gi  +g2. 


where  g represents  the  total  hole  generation  process.  With  the 
condition  of  eq.  (3.36),  eq.  (3.39)  readily  follows. 


3.3.  Device  Description.  Experimental  Results  and  Discussion 
In  this  section  we  outline  how  we  used  the  above  procedures  to 
extract  trap  parameters.  When  the  electronic  properties  of  the  device 
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do  not  fit  the  deep  defect  density  approximation  made  in  the  previous 

section,  an  exact  numerical  analysis  should  be  employed. 

+ - + 

A p -p  -p  resistor  (M102-22  #1)  shown  in  figure  3.1  is  used  in 
this  experiment.  An  epi-layer  of  thickness  10.4  um  is  grown  on  a 
highly  doped  substrate  with  a resistivity  of  0.02  Ohm-cm.  The  boron 
concentration  in  the  epi-layer  is  2.7  x 1015  cm‘3. 

If  the  gold  concentration  is  Na  > Nau  > 0.5Na  and  the  shallow 
impurities  are  completely  ionized,  the  gold  impurities  are  either 
doubly  ionized  or  singly  ionized,  and  the  Fermi  level  will  be  very  close 
to  the  second-donor  energy  level.  The  device  performance  will  be 
affected  by  both  the  first  and  the  second  donor  level. 

In  chapter  II,  we  investigated  a p+-p’-p+  structure  (device 
CM24B-82  #1)  made  of  Si:Au:B  with  a boron  concentration  of  1.0  x 
10  cm  . In  that  particular  case  the  resistance  and  noise  data  versus 
1000/T  showed  identical  activation  energies  of  0.34  eV.  From  the 
experimental  data  a gold  concentration  2.51  x 1015  cm'3  was 
extracted.  In  the  present  work  we  use  again  a p+-p"-p+  structure 
(device  M10B2-22  #1)  but  this  time  with  a shorter  p"  region  and 
higher  boron  doping.  An  epi-layer  of  thickness  10.40  um  was  grown  on 
a highly  doped  substrate  with  a resistivity  of  0.02  Ohm-cm.  The  boron 
concentration  in  the  epi-layer  is  2.7  x 1015  cm'3.  The  thickness  of  the 
wafer  is  7.34  mil.  The  drive-in  temperature  is  900°  Celsius;  the  drive- 
in  time  is  10.75  hours.  The  junction  diameter  of  the  device  is  7.5  x 
10  2 cm  and  the  calculated  junction  depth  of  the  Ohmic  contact  is 
1.93  um.  The  resistance  of  the  device  at  temperature  T=300.8  K is 
41.76  Ohm.  The  plot  of  the  resistance  versus  1000/T  shows  an 
activation  energy  of  0.267  eV  as  is  indicated  in  figure  3.3.  The 
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normalized,  low-frequency  noise  spectral  density  plotted  versus 
1000/T  in  figure  3.4  on  the  other  hand  shows  an  activation  energy  of 
0.336  eV.  A combination  of  figure  3.3  and  figure  3.4  is  shown  in  figure 
3.5.  As  we  can  clearly  see  from  this  figure,  the  values  of  the  resistance 
show  a smaller  activation  energy.  The  externally  applied  voltage  during 
the  noise  measurement  was  about  0.12V  which  ensured  that  the 
device  operated  in  the  Ohmic  regime.  The  voltage  noise  spectral 
density  was  measured  between  100  and  1000  Hz. 

Since  the  gold  diffusion  conditions  for  the  devices  discussed  in 
Chapter  II  and  this  chapter  were  the  same,  but  the  thickness  of  the 
wafers  and  the  length  of  the  p~  region  of  the  devices  differed,  we 
expect  the  gold  concentration  in  M10B2-22  #1  to  be  larger  than  2.51 
x 10  cm  . We  have  no  reason,  however,  to  believe  that  it  is  larger 

ic  q 

than  2.7  x 10  cm'  , since  then  the  resistance  versus  1000/T  for  this 
device  is  expected  to  have  the  same  slope  as  in  Chapter  II,  which  is 
not  in  agreement  with  the  data  shown  in  figure  3.3.  On  the  other  hand 
if  the  concentration  of  gold  is  smaller  than  2.7  x 1015  cm'3  and  gold 
has  only  one  donor  level,  then  the  resistance  of  the  device  is  expected 
to  decrease  as  temperature  decreases  due  to  an  increase  of  the 
mobility.  Again  this  is  not  what  we  measured.  Therefore,  we  conclude 
that  gold  in  silicon  produces  a second-donor  level  and  that  in  our 
sample  the  condition  Na  > Nau  > 0.5  Na  is  met. 

Using  the  optimization  procedures  outlined  in  Chapter  II  and 
the  semi-empirical  mobility  formula  (2.19),  we  calculate  the  hole 
density  P from  the  charge  neutrality  equation  (3.6)  and  gold 
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1000/T  (K  ) 


Figure  3.3.  The  resistance  versus  lOOO/T  of  device  M10b2-22  #1. 
The  squares  represent  the  experimental  data.  The  dashed  line 
shows  the  best  fit  to  the  data. 
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Figure  3.4.  The  relative,  low-frequency  G-R  noise  plateau  values 
versus  1000/T.  The  circles  represent  the  experimental  data. 
The  dashed  line  shows  the  best  fit  to  the  experimental  data. 
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Figure  3.5.  The  relative,  low-frequency  G-R  noise 
plateau  values  and  the  resistance  versus  1000/T. 
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conservation  equation  (3.3)  exactly,  and  fit  the  noise  and  the 
resistance  data  to  eqs.  (3.7)  and  (3.40),  and  extract  the  trap 
parameters  listed  in  Table  3.2.  The  extracted  hole  capture  cross 
section  of  Opi  = 6.5  x 10  15  cm^  is  consistent  with  the  hole  capture 
cross  section  of  the  first-donor  level  published  in  the  literature.  The 
values  of  the  other  parameters  listed  in  Table  3.2  justify  the 
assumptions  made  in  the  previous  section.  The  accuracy  of  the  fitted 
results  is  limited  by  the  mobility  formula  we  used,  the  experimental 
accuracy  and  the  fact  that  we  assumed  the  activation  energies  of  the 
gold  donor  and  the  hole  capture  cross  section  of  the  first  donor  to  be 
temperature  independent.  The  above  data  may  therefore  not 
represent  the  exact  physical  quantities.  In  fact.  Eld  and  E2d  may  vary 
slightly  with  temperature  because  they  contain  a contribution  from  the 
entropy  and  because  Ev  varies  slightly  with  temperature.  These 
temperature-dependent  factors  have  been  accounted  for  in  the 
degeneracy  factors. 


3.4.  Summary 

At  low  temperatures,  when  Na  > Nau  > 0.5Na,  doubly  ionized 
gold  centers  become  electrically  active  and  the  Fermi  energy  level  is 
located  between  E2d  and  Ev.  The  average  number  of  carriers  in  the 
valence  band  is  then  controlled  by  the  second  gold  donor  level. 
However,  since  the  interaction  between  the  mobile  hole  and  the 
second-donor  level  is  repulsive,  the  hole  generation-recombination 
rates  involving  the  second  donor  level  will  be  smaller  than  the  rates 
associated  with  the  first  gold  donor  level.  The  second-donor's 
contribution  to  the  noise  will  therefore  be  small  and  the  observed 
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TABLE  3.2.  Trap  parameter  results  obtained  from  extraction 
procedure. 


Device  parameters 

Magnitude 

Unit 

^pi 

6.5  x io~15 

cm2 

Nau 

2.58  xio15 

cm'3 

^aul 

16 

^au2 

60.2 

Eld  - Ev 

0.336 

eV 

Eld  - Ev 

0.267 

eV 
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generation-recombination  noise  will  mainly  come  from  the  following 
processes:  (a)  hole  emission,  Nau+  ->  Nau°  + e+  and  (b)  hole  capture, 
Nau  + e ->  Nau  , This  is  why  the  resistance  and  noise  data  show 
different  activation  energies. 

From  the  above  discussion,  it  becomes  clear  that  the  activation 
energy  of  the  resistance  is  associated  with  the  second  gold  donor 
which  lies  0.267  eV  above  the  valence  band  edge,  and  that  this  level 
can  be  detected  only  within  a certain  window  of  doping  densities,  i.e., 
Na  > Nau  > o.5Na. 

3.5.  Comparison  of  Experimental  Results  for  the  Single-  and  Double- 

Donor  Model 

In  Chapter  II  and  in  sections  of  this  chapter,  we  discussed  trap 
parameter  extraction  using  a single-donor  and  double-donor  model. 
The  facts  are  that  the  noise  measured  on  devices  (Ml 02-22  #1) 
and  ( CM24B-82  #1)  are  produced  by  the  same  gold  centers,  and  that 
the  double  donor  model  applies  to  both  of  the  devices.  Using  the  same 
optimization  procedures  as  outlined  in  the  previous  chapter,  the 
double  donor  model  is  fitted  to  the  experimental  data  of  Chapter  II 
and  this  chapter,  and  the  final  results  are  shown  in  figure  3.6.  The 
agreement  is  good  in  terms  of  the  simplified  models,  except  that  the 
hole  capture  cross  sections  for  the  two  devices  are  different:  the  hole 
capture  cross-section  Opi  = 6.5  x io'15  cm2  for  device  (M102-22  #1) 
while  the  hole  capture  cross-section  dpi  = 3.217  x io15  cm2  for 
device  (CM24B-82  #1).  Other  parameters  are  listed  in  Table  3.2.  The 
difference  in  the  hole  capture  cross-sections  may  be  caused  by  the  fact 
that  the  device  wafers  are  supplied  by  different  manufacturers.  Device 
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(CM24B-82  #1)  uses  a wafer  supplied  by  Cincinetti  while  device 
(M 102-22  #1)  is  made  on  a wafer  supplied  by  Monsanto. 
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Figure  3.6.  The  normalized,  low-frequency  G-R  noise 
plateau  values  and  the  resistance  vs  1000/T. 
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CHAPTER  IV 

THE  EFFECT  OF  FIELD- DEPENDENT  EMISSION  ON  THE 

CURRENT-VOLTAGE  CHARACTERISTICS  OF  A P+-P'-P+  Si:Au:B 

DEVICE 

4.1.  Introduction 

Theoretical  and  experimental  investigations  of  the  current- 
voltage  characteristics  of  high-low-high  structures  made  of  semi- 
insulating  semiconductor  devices  with  deep  level  doping  have  been 
carried  out  for  decades.  An  excellent  review  of  carrier  injection  in 
these  materials  is  presented  in  the  book  by  Lampert  and  Mark  [28]. 
The  current-voltage  (I-V)  behavior  in  high-resistivity  materials  can  be 
described  by  (1)  Ohm's  law,  (2)  the  trap  free  law,  and  (3)  the  trap- 
filled  limited  law  [28].  The  current  in  these  devices  is  predominantly 
due  to  drift  and  might  be  affected  by  hot  carrier  phenomena. 
However,  these  laws  are  no  longer  sufficient  to  describe  all  aspects  of 
the  I-V  characteristics  in  the  domain  where  the  capture  and  emission 
of  carriers  are  field-dependent.  Both  superlinear  and  sublinear 
behavior  can  result  from  these  competing  mechanisms.  Our 
experiments  indicate  that  in  the  devices  we  investigated  in  Chapter  II 
sublinear  behavior  prevails.  We  therefore  develop  in  this  chapter  a 
model  for  sublinear  I-V  behavior  and  do  not  address  superlinear 
aspects. 

The  occurrence  of  sublinear  I-V  behavior  not  associated  with  hot 
carrier  mobility  degradation  has  been  reported  by  several  researchers 
[29-32].  The  possible  mechanisms  responsible  for  this  phenomenon 
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are  carrier  exclusion  [30,31]  and  field  enhanced  carrier  capture  into 
deep  trap  levels  [32].  In  steady  state,  the  carrier  concentration  in  a 
semi-insulating  device  is  determined  by  the  balance  of  emission  and 
capture  processes,  in  addition  to  junction  effects.  As  a consequence, 
either  a field-induced  decrease  in  the  emission  coefficient,  or  a field- 
enhanced  capture  coefficient,  or  both  can  be  reasons  for  observing 
sublinear  I-V  characteristics. 

Most  of  the  previous  investigations  on  nonlinear  I-V 
characteristics  were  based  on  I-V  curves  only.  In  Chapter  II  and  III,  we 
have  used  a combination  of  resistivity  and  noise  data  to  explore  the 
defect  properties  of  gold  in  silicon,  and  we  have  shown  that  low-field 
device  parameters  can  be  extracted  successfully  by  this  method.  In 
this  chapter,  we  extend  our  noise  measurement  to  a moderate  high- 
field  region  to  help  unravel  the  generation-recombination  effects  on 
the  I-V  characteristics. 

We  measured  the  low  temperature  I-V  characteristics  of  a 
+ - + 

p -p  -p  device  made  of  Si:Au:B  at  various  temperatures  between  150 
and  200K.  The  low-voltage  region  of  the  I-V  characteristics  shows 
Ohmic  behavior,  and  is  followed  by  a sublinear  region  at  higher  voltage. 
With  the  use  of  noise  data,  we  found  that  the  sublinear  I-V  behavior  in 
our  device  is  due  to  a significant  decrease  in  the  hole  emission 
coefficient  of  gold  donor  centers  with  increasing  bias. 

In  section  4.2,  we  describe  the  method  we  employed  to 
interpret  the  I-V  characteristics  and  noise  spectra  observed  in  our 
experiment.  In  section  4.3,  the  experimental  details  are  presented 
and  discussed. 
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4.2.  Theory 

In  this  section,  we  describe  the  theory  used  to  interpret  the 

sublinear  I-V  characteristics  and  the  noise  data.  The  device  we 

+ - + 

investigated  was  a p -p  -p  silicon  resistor,  with  boron  dopants  and  a 

large  gold  concentration.  We  assume  that  the  device  is  homogeneous 

and  that  the  electric  field  is  independent  of  position.  In  addition  we 

assume  that:  (1)  only  the  gold  donors  (0.336  eV  above  the  valence 

band  edge)  interact  with  the  valence  band,  (2)  the  electron  density 

can  be  neglected,  (3)  the  gold  and  boron  densities  are  uniform,  (4)  the 

gold  concentration  Nau  is  larger  than  the  boron  concentration  Na,  (5) 

diffusion  currents  can  be  neglected. 

The  I-V  characteristics  are  then  given  by  eq.  (3.7).  The  rate 

equation  takes  the  form 

dP  + + 

= ep  Nau  - cp  P (Nau  - Nau  ),  (4. 1) 

where  Nau  is  the  charged  gold  donor  concentration,  cp  is  the  hole 

capture  coefficient,  and  ep  is  the  hole  emission  coefficient. 

Next  we  expand  the  fluctuating  parameters  P and  Nau+  around 

their  steady  state  values.  Neglecting  second  and  higher  order  terms 

we  write 

P = P0  + AP 

and 

+ + + 

Nau  = Nauo  + ANau  , 

where  the  subscript  o indicates  steady  state  and  AP  and  ANau+  are  the 
time-dependent  fluctuations  in  P and  Nau+  respectively. 

Substituting  the  above  expressions  in  equation  (4.1)  we  derive 
after  separating  d.c.  and  a.c.  terms 
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eP  = cp  (Nau  - Nauo+  )( ■ 0 +),  (4.2) 

F Nauo 

dAP 

“£[£“=  - [ep  + cp  PQ  +cp  (Nau-  Nauo+)]AP.  (4.3) 

In  the  derivation  of  equation  (4.3)  we  used  that  AP  s - ANau+,  which 

can  be  justified  by  considering  the  condition  for  space  charge 

neutrality,  P + Nau  = Na , or  in  terms  of  density  fluctuations 
+ 

AP  + ANau  = ANa  . Our  low-frequency  noise  experiment  does  not 
probe  the  rapid  fluctuations  associated  with  the  trapping  and 
detrapping  of  holes  in  the  shallow  acceptor  level,  consequently 
we  may  assume  ANa  s 0 or  AP  = - ANau+. 

In  our  overcompensated  p+-p'-p+  device,  and  at  the 
temperatures  selected  for  our  experiments,  150<  T < 250K,  the 
equilibrium  Fermi  level  lies  below  the  deep  gold  donor  level  and  far 
above  the  shallow  acceptor  level,  resulting  in  equilibrium  values  for  P 
between  109  cm'3  and  1012  cm'3,  specifically  P=2.2  xlO11  cm'3  at 
T=200K.  The  equilibrium  value  of  Nau+  is  of  the  order  of  the  value  of 
Na  since  the  gold  donor  electrons  will  occupy  the  shallow  acceptor 
sites.  Due  to  the  field  dependence  of  ep  and  cp,  P and  Nau  will  deviate 
from  their  equilibrium  values  when  an  electric  field  is  applied.  As  we 
will  demonstrate  later,  the  changes  in  ep  and  cp  in  the  field  range  we 
investigated  are  relatively  minor  so  that  P and  Nau+  do  not  change 
much  from  their  equilibrium  values.  Consequently  we  may  assume, 
using  eq.  (4.2)  and  the  charge  neutrality  condition,  that 

[ep  + cp  PQ  ] « [cp  (Nau-  Nauo+)],  (4.4a) 

and 

Nauo  = Nao  = Na- 


(4.4b) 
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These  conditions  simplify  equation  (4.3)  to 


dAP 

dt 


AP 

= - cp  (Nau-  Na)  AP=  - -r — , 


(4.5) 


with 


T~=  cp  (Nau  - Na). 

Lp 


(4.6) 


The  normalized,  plateau  value  of  the  G-R  noise  has  the  form 
S^O)  4 Nauo+  ep  tp2 


V2-  LwdPp^ 

With  the  use  of  the  previous  assumptions  and  equations,  eq.  (4.7)  can 
be  written  as 


(4.7) 


sw(°)  _ 4 

V2  LwdepNau  ' (4.8) 

If  we  now  know  the  hole  mobility  as  a function  of  electric  field, 
then  the  hole  concentration  P can  be  calculated  from  eq.  (3.7)  and  the 
hole  emission  coefficient  from  eq.  (4.8).  In  addition  the  value  of  cp  can 
be  derived  from  eq.  (4.2). 

The  extraction  of  device  parameters  requires  a model  for  the 
hole  mobility  since  this  quantity  is  not  independently  measured  in  our 
experiment.  The  experimental  drift  velocity  versus  electric  field  curve 
of  lightly  doped  silicon  for  moderate  electric  field  values  can  be  found 
in  the  literature  [33].  The  measured  field  variation  of  velocity  can  be 
approximated  with  the  following  empirical  expression  proposed  by 
Caughey  and  Thomas  [34], 


V =U 


E_ 


E P 1/P 
(1+(|H  ) 

c 


(4.9) 
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where  Us  is  the  saturation  velocity,  (3  is  a fitting  parameter,  Ec  is  a 
critical  field  given  by  Ec  = t)s/p0,  and  p0  is  the  low  field  carrier 

mobility. 

Using  \>s=0.85lxl07cm/s  and  po=1280  cm2/(Vs)  [35],  we  fitted 

eq.  (4.9)  to  the  experimental  data  [33],  and  obtained  a best  fit  for  p 
= 1.252  at  T=200.0K.  The  agreement  between  the  experimental  data 
and  eq.  (4.9)  is  shown  in  Fig.  4.1.  The  average  error  is  5%. 

Equation  (4.9)  is  strictly  valid  only  for  lightly  doped  silicon.  In 
order  to  estimate  the  effect  of  ionized  impurities  on  the  hole  mobility, 
we  use  the  low-field,  semi-empirical,  mobility  formula  (2.19)  to 
estimate  the  effects  of  a charged  impurity  concentration.  We  found 
from  this  equation  that  the  effect  of  ionized  impurity  scattering  is 
negligible  for  concentrations  less  than  3.0  x 1015  cm'3.  Since  our 
devices  have  an  ionized  impurity  concentration  of  the  order  of 
2.0x10  cm  , mobility  degradation  due  to  impurities  can  be 
neglected. 


4.3.  Experimental  Details  and  Discussion 
A p+-p'-p+  structure  (CM24B-82  #2)  made  of  Si:Au:B  with  a 
boron  concentration  of  1.0  xlO15  cm'3  is  used  in  this  work.  The 
device  details  are  described  in  Chapter  II.  From  the  experimental 
data,  a gold  concentration  of  2.51  x 1015  cm  3 was  extracted.  The 
calculated  length  of  the  p'-device  region  is  18.1  um.  The  resistance  of 
the  device  at  T=298.9  K is  69.9  Ohm.  The  bias  current  I,  supplied 
from  a battery,  is  controlled  by  a variable  series  metal  film  resistor. 
The  voltage  noise  spectral  density  was  measured  between  100  and 
1000  Hz  and  was  independent  of  frequency.  The  device  was  placed 
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Figure  4.1.  Hole  drift  velocity  versus  electric  field  at  T=200.0K 
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inside  a cryostat  which  allowed  us  to  vary  the  device  temperature 
between  77  and  300  K.  Calculations  show  that  the  conditions  for 
which  the  simple  noise  description  given  in  section  4.2  are  valid  are 
met  at  temperatures  between  150  and  250  K.  We  therefore 
concentrate  on  this  temperature  regime.  The  I-V  characteristics  at 
T=200.0  K measured  as  a function  of  bias  polarity  are  shown  in  Fig. 
4.2.  Figure  4.3  shows  the  normalized  plateau  value  of  the  G-R  noise 
measured  as  a function  of  bias  voltages  at  T=200.0  K.  The 
measurements  performed  at  other  temperatures  between  150  and 
250  K show  qualitatively  similar  behavior.  We  focus  on  the  200.0  K 
data  however,  since  at  this  temperature  accurate  drift  velocity  data  are 
available.  In  the  moderate  field  regime  of  our  investigations 
(E<1. 5 kV/cm)  the  noise  spectral  density  was  independent  of  bias 
polarity. 

As  the  voltage  across  the  device  increases,  the  electric  field 
inside  the  device  increases  and  the  hole  mobility  will  decrease  at 
moderate  to  high  electric  field  strength.  However,  the  mobility 
reduction  is  very  small  in  the  field  strength  range  of  0 to  1.5  kV/cm 
where  the  I-V  characteristics  deviate  from  ohmic  behavior.  Space 
charge  injection  is  also  negligible  at  these  field  strengths  due  to  the 
large  deep  gold  concentration.  Hence,  the  sublinear  I-V 
characteristics  must  be  caused  by  other  mechanisms  among  which  are 
junction  effects  and  generation-recombination  processes. 

In  order  to  ensure  that  Joule  heating  was  negligible  and  not 
responsible  for  the  observed  sublinearity  in  the  I-V  characteristics,  we 
measured  the  noise  data  at  T=192.5  K where  the  resistance  of  the 
device  is  significantly  larger,  and  subsequently  Joule  heating  is  less 
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Figure  4.2.  The  current  voltage  characteristics  at  T=200.0K 
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Figure  4.3.  The  normalized,  low  frequency  G-R  noise  plateau 
versus  voltage  and  electric  field  at  T=200.0  K 
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Figure  4.4.  The  normalized,  low  frequency  G-R  noise  plateau 
versus  voltage  and  electric  field  at  T=192.5K 
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important.  The  noise  data  are  shown  in  Fig.  4.4.  The  fact  that  the  noise 
behavior  observed  at  both  temperatures  is  similar  proves  that  our 
measurements  were  not  affected  by  Joule  heating. 

Sublinear  I-V  characteristics  observed  in  n+  -n'  - n+  structures 
made  of  Si:Au  are  attributed  to  carrier  exclusion  [30,31].  Carrier 
exclusion  occurs  at  a reverse-biased  high-low  junction  in  a device  and 
reduces  the  majority  carrier  density.  A condition  for  exclusion  to 
occur  is  a strong  asymmetry  of  the  electron  and  hole  lifetimes. 
Exclusion  effects  can  therefore  not  be  used  to  account  for  the 
sublinear  I-V  characteristics  of  our  device,  since  at  low  temperatures 
the  lifetimes  of  the  hole  and  electron  carriers  are  of  the  same  order, 
and  the  minority  carrier  density  is  negligible. 

Hence,  we  opt  for  an  interpretation  of  the  sublinear  I-V 
characteristics  in  terms  of  field-dependent  generation-recombination 
processes.  At  low  temperatures,  the  hole  concentration  is  smaller  than 
the  gold  concentration,  and  therefore  a variation  of  the  gold  hole 
emission  and  capture  coefficient  with  electric  field  will  significantly 
affect  the  hole  density  in  the  device.  Using  the  equations  presented  in 
the  previous  section,  we  extract  from  the  experimental  data  the  hole 
emission  coefficient  and  the  capture  coefficient.  The  results  are 
presented  in  Fig.  4.5  and  show  that  the  emission  coefficient  is 
independent  of  electric  field  in  the  low  field  region  and  then  starts  to 
decrease  as  the  applied  voltage  approaches  1 V.  In  the  same  region, 
the  I-V  characteristics  show  sublinear  behavior.  Because  the  field- 
induced  decrease  in  the  hole  emission  coefficient  is  larger  than  the 
field-induced  decrease  in  the  hole  capture  coefficient,  the  net  effect  is 
that  the  hole  carrier  density  in  the  valence  band  is  reduced  as  the 
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Figure  4.5.  The  hole  emission  and  capture  coefficients 
versus  voltage  at  T=200K 
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electric  field  increases,  resulting  in  the  observed  sublinear  I-V 
characteristics. 


4.4.  High  Field  Effects 

In  the  preceding  sections,  we  have  shown  that  the  I-V 
characteristics  are  independent  of  device  polarity  at  medium  electric 
fields.  At  higher  electric  field  strengths  however,  the  device  shows  a 
bias  polarity  dependence.  Figure  4.6  depicts  the  I-V  characteristics  of 
device  (CM24B-82  #2).  Positive  polarity  refers  to  a positive  bias  to  the 
Ohmic  contact  at  the  diffusion  region.  Figures  4.7  and  4.8  represent 
the  I-V  characteristics  of  device  (Ml 02-22  #1).  The  asymmetry  of  the 
device  geometry  and  the  doping  profiles,  combined  with  the  unknown 
field  dependence  of  the  defect  parameters  will  affect  the  I-V 
characteristics  in  the  high  field  regime  and  make  it  therefore  difficult 
to  analyze  the  measured  polarity  dependence  of  the  I-V  characteristics 
and  the  noise  data  and  to  extract  trap  parameters  for  this  electric  field 
regime. 
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Figure  4.6.  The  I-V  characteristics  of  device  (CM24B-82#2) 
at  T=192.5K. 
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Figure  4.7.  The  I-V  Characteristics  of  device  (CM10B2-22#1) 
at  T=  156.0  K.  The  solid  line  represents  the  results  of  the 
measurements  with  negative  polarity,  the  crosses  refer  to 
positive  polarity  results. 
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Figure  4.8.  The  I-V  Characteristics  of  device  (CM10B2-22#1) 
at  T=  156.0  K.  The  solid  line  represents  the  results  of  the 
measurements  with  negative  polarity,  the  crosses  refer  to 
positive  polarity  results. 


CHAPTER  V 

LOW-FREQUENCY  GENERATION-RECOMBINATION  NOISE  AS  A TOOL 
FOR  INVESTIGATING  DEEP  DEFECTS,  DEVICE  DESIGN 

CONSIDERATIONS 

5.1.  Introduction 

Electronic  devices  and  systems  exhibit  random  fluctuations  in 
voltage  (or  current)  at  their  terminals.  These  fluctuations  are  referred 
to  as  electronic  noise.  Since  noise  reflects  the  microscopic  motion  of 
carriers,  it  has  been  used  as  a tool  for  investigating  the  electrical 
characteristics  of  a device.  Generation-recombination  (G-R)  noise  is 
due  to  carrier  transitions  between  traps  and  their  associated 
conduction  or  valence  band,  and  it  can  be  used  to  investigate  the 
properties  of  a defect  in  a semiconductor. 

Up  to  now,  most  of  the  investigations  of  defects  are  done  using 
transient  capacitance  methods  [11].  These  methods  investigate 
transient  capture  and  transient  emission  processes  as  two 
independent  events.  However,  most  problems  of  interest  in  device 
application  are  in  steady-state  or  quasi-steady  state  situations  in  which 
generation  and  recombination  processes  occur  simultaneously. 
Therefore  the  data  from  steady-state  G-R  noise  measurements  will 
provide  without  a doubt  useful  information  concerning  G-R 
mechanisms. 

In  the  past,  most  G-R  noise  studies  [4-6]  have  focused  on  traps 
producing  a single  energy  level.  Little  work  has  been  done  on  multi- 
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level  G-R  traps.  Since  defects  can  exist  in  several  charge  states  in  a 
device,  other  states  may  affect  the  measured  data  collected  for  a 
particular  energy  level.  It  is  essential  to  understand  the  effects  of 
other  charge  states  on  the  observed  energy  level  in  order  to  obtain  an 
accurate  interpretation  of  measured  data.  To  our  knowledge,  no  work 
has  been  done  on  discrimination  between  donor  and  acceptor  behavior 
of  a defect  energy  level  with  G-R  noise  measurements.  In  Chapters  II 
and  III,  we  investigated  G-R  noise  on  samples  with  nearly-  and  over- 
compensated, deep  impurity  concentrations,  and  discovered  a second- 
donor  level  of  gold  in  silicon  [Chapter  III].  These  investigations  have 
motivated  us  to  further  explore  this  technique  as  a deep  defect 
diagnostic  tool. 

Because  device  design  is  critical  for  the  successful  investigation 
of  defects  with  G-R  noise  measurements,  this  chapter  is  written  to 
provide  information  for  designing  a bulk  semiconductor  device  for 
investigating  the  properties  of  defects.  In  the  following  section,  we 
present  several  G-R  noise  models  for  defects  with  either  one,  or  two 
energy  levels  in  the  lower  half  of  the  band  gap  of  a p-type 
semiconductor,  and  study  several  special  cases  which  make  an 
analytical  solution  and  an  insightful  device  design  scheme  possible. 
These  cases  are:  (a)  the  energy  levels  of  a defect  are  separated  by  more 
than  a few  KT  from  each  other  and  (b)  the  energy  levels  of  a defect  are 
within  a few  KT  from  each  other.  We  layout  the  general  theoretical 
background  in  section  5.2  and  then  explore  various  models.  Finally, 
guidelines  for  proper  device  design  are  given. 
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5.2.  Theory 

In  this  chapter,  we  investigate  the  properties  of  defects  with 
energy  levels  located  in  the  lower  half  of  the  band  gap  of  a p-type 
semiconductor.  The  resistance  and  voltage  noise  resulting  from  deep 
defects  interacting  with  the  valence  band  are  calculated  for  the  low 
field  ohmic  regime  of  device  operation  and  are  based  on  the  following 
assumptions:  (1)  space  charge  neutrality  and  space  uniformity  prevail, 
(2)  the  electron  density  is  neglected,  (3)  shallow  impurities  are  fully 
ionized  (Na  = Na). 

According  to  W.  Shockley  and  J.  T.  Last  [26],  the  ionized  defect 
densities  for  a multi-level  defect  can  be  expressed  as 

Ef  -Ela 

Nf  = glaNi°  e KT  , (5.1.1) 


Ni+  = -7—  Ni°  e 


Eld"  Ef 

TT 


&ld 
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&2d 


Ni++  = -r-  Ni°  e 


Eld+  E2d  2 Ef 

ITT 


(5.1.2) 

(5.1.3) 


O + + + 

where  Ni,  Ni  , Ni  , Ni  . Ni  are  the  total,  singly  negatively  ionized, 
neutral,  singly  positively  ionized,  and  doubly  positively  ionized  defect 
concentration  respectively.  The  symbols  gld,  g2d*  Sia  are  the 

degeneracy  factors  for  the  first-donor,  second-donor,  and  the  acceptor 
level  respectively.  The  energies  Ela,  Eld,  E2d  are  the  defect  acceptor, 
first-donor,  second-donor  energy  level  respectively,  and  Ef  is  the 

Fermi  energy  level.  In  this  chapter,  we  restrict  our  discussion  to  a 
defect  with  energy  levels  satisfying  Ela  > Eld  > E2d,  such  as  for  example 

gold  in  silicon. 

In  order  to  calculate  the  noise  spectral  density,  we  express  the 
carrier  densities  as  sums  of  a steady-state  term  and  a fluctuation  term 
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Ni°  = Nio°  + ANi°, 

(5.1.4) 

Ni+=  Nio+  + ANi+, 

(5.1.5) 

Ni'=  Nio'  + ANi', 

(5.1.5b) 

Ni++=  Nio++  + ANi++, 

(5.1.6) 

P = P0  +AP, 

(5.1.6b) 

where  Nio  , Ni0°,  Nio+,  Nio++  are  the  singly  negatively  charged, 
neutral,  singly  positively  charged,  and  doubly  positively  charged  defect 
densities  in  steady  state  respectively,  and  ANio  , ANi0°,  ANi0+,  ANio++ 
are  the  corresponding  fluctuations. 

The  carrier  generation  and  recombination  rates  associated  with 
a multi-level  defect  in  a semiconductor  are  described  in  Table  5.1.  The 
parameter  P is  the  hole  density,  and  N is  the  electron  concentration. 
The  symbol  e+  denotes  a mobile  hole,  and  gj  and  rl5  where  i=  1,2,3,  are 
the  hole  generation  and  recombination  rates  respectively.  The  trap 
parameters  epo,  cpo,  epi,  cpi,  and  ep2,  cp2  are  the  emission  and  capture 
coefficients  for  holes  in  the  acceptor,  first-donor  and  second-donor 
level  respectively. 

In  equilibrium,  with  the  use  of  eqs.  (5.1.1)  through  (5.1.3)  and 
the  principle  of  detailed  balance,  we  get 
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(5.1.7) 

(5.1.8) 

(5.1.9) 


where  the  subscript  e denotes  the  equilibrium  value. 
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TABLE  5.1.  Generation  and  recombination  rates  in  terms  of 
ionized  defect  and  mobile  hole  concentrations. 


Rate 

Hole  Transitions  involved 

Description  of  the  rates 

& 

Ni°  ->  Ni’  + e+ 

epo  Ni° 

rQ 

Ni  + e+  ->  Ni° 

cpo  Ni  P 

gi 

Ni+  ->  Ni°  + e+ 

epi  Ni+ 

0 , + 

Ni  + e ->  Ni 

cpi  Ni°  P 

62 

++  + , 
Ni  ->  Ni  + e 

ep2  Ni++ 

12 

+ . ++ 
Ni  + e ->  Ni 

cp2  Ni+  P 

77 


5.2.1.  Single  Level  Noise  Model 

We  first  deal  with  the  G-R  noise  model  for  a defect  with  a single 
donor  level  in  a device  with  a shallow  background  doping.  The 
following  results  have  been  derived  in  Chapter  II.  With  the  carrier 
densities  given  by  eqs.  (5.1.6b)  and  (5.1.5),  the  condition  of  charge 
neutrality  gives  the  fluctuations  in  the  carrier  densities  at  low 
frequencies 

AP  + ANi+=0.  (5.2.1) 

The  rate  equation  takes  the  form 
dP  + + 

cfi"  = epi  Ni  - CP1  p (Ni  - Ni  )•  (5.2.2) 


Using  eqs.  (5.1.5),  (5.1.6b)  and  (5.1.8),  we  get  after  linearization 


dAP 

dt  = ‘ ^p1  + cp1  Po  ^p1 


(Ni-Nio+)]AP. 


(5.2.3) 


This  equation  can  be  further  simplified  as  we  will  show  later. 


5.2. 1.1.  Deep  Defects  in  Overcompensated  Devices 

If  the  deep  donor  defect  concentration  is  larger  than  the  shallow 
impurity  concentration,  the  Fermi  level  is  located  close  to  the  defect 
donor  energy  level.  At  low  temperatures,  P is  small,  we  may  assume 
P « (Ni-  Nio  ),  Nio  , Ni. 

With  the  use  of  eq.  (5.1.8),  we  can  assume  that 

[epi  + cpi  P0  ] « [Cpi  (Ni  -Nio+)].  (5.2.4) 

Then  eq.  (5.2.3)  simplifies  to 
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Assuming  Na  = Nio+,  eq.  (5.2.6)  becomes 

1 o 

— = Cpi  Nio  = cpi  (Ni  - Na).  (5.2.7) 

The  temperature  dependence  of  the  capture  coefficient  can  be 
derived  from  this  equation.  The  normalized  spectral  density  of  G-R 
noise  has  the  form  [17] 
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where  Nio+  epi  is  the  steady-state  free  carrier  generation  rate. 
The  normalized  plateau  value  of  the  G-R  noise  is 
Sw(Q)  _ 4 Nio+  ePi  tp2 


LwdP„ 


LwdP0  cpi  Ni 

At  low  temperatures  P can  be  written  as  [Chapter  II] 
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(5.2.10) 


Ni-  Na 

with  the  aid  of  the  charge  neutrality  equation 

Na’  = P + Ni+.  (5.2.11) 

Using  eqs.  (5.1.8),  (5.2.7),  (5.2.9)  and  (5.2.10),  and  assuming  the 
hole  density  is  small,  we  get 
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Using  eqs.  (3.7)  and  (5.2.10),  we  get 
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Equations  (5.2.12)  and  (5.2.13)  show  that  the  activation  energies 
extracted  from  these  equations  should  have  the  same  value. 


5.2. 1.2.  Deep  Defects  in  Uncompensated  Devices 

If  the  deep  impurity  concentration  is  smaller  than  the  shallow 
impurity  concentration,  we  may  assume  from  eq.  (5.2.6)  that 

P » Ni+,  Ni.  (5.2.14) 

Equations  (5.2.5)  and  (5.2.6)  simplify  to 
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Combining  eqs.  (2.23),  (5.1.8)  and  (5.2.15),  we  get 
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When  the  Fermi  level  is  located  above  the  trap  level,  the  above 
equation  becomes 
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(5.2.16b) 


The  trap  energy  can  be  extracted  from  this  equation  by  plotting 
Ln(tpT2)  versus  1000/T  where  it  is  assumed  that  the  temperature 
dependence  of  the  hole  capture  cross  section  is  weak. 

From  eq.  (5.2.8),  we  get 


sw(0)  SPP(°)  4 Nio+  ep  i tp 
tDV  2 " tDP2  ” LwdP2 


(5.2.17) 


The  trap  energy,  degeneracy  factor,  trap  concentration  can  be 
extracted  using  the  previous  equations  as  discussed  in  section  2.2.2. 
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5.2. 1.3.  Device  Design  Consideration 

The  magnitude  of  G-R  noise  increases  when  the  Fermi  level 
approaches  the  trap  energy  level.  This  condition  is  critical  for  a device 
to  be  used  for  trap  parameter  extraction.  For  an  overcompensated 
device,  P is  very  small  at  low  temperatures.  Using  eqs.  (2.23)  and 
(5.2.10),  the  Fermi  level  can  be  expressed  as 


Once  a device  is  overcompensated,  the  Fermi  level  will  be  close  to  the 
trap  energy  level.  As  the  temperature  is  lowered,  the  Fermi  level  gets 
closer  to  the  trap  energy  level,  and  the  noise  signal  will  therefore 
increase.  The  only  design  criterion  to  be  considered  is  to  have  a 
device  with  suitable  values  of  resistance  for  low-frequency  noise 
measurements. 

For  a device  that  is  lightly  doped  with  deep  defects,  and  in 
which  the  deep  defect  concentration  is  smaller  than  the  shallow 
acceptor  concentration,  the  Fermi  level  is  controlled  by  the  shallow 
impurity  concentration  and  temperature.  If  the  shallow  dopants  are 
almost  fully  ionized,  then  the  hole  density  P can  be  written  as 


(5.2.18) 


= Na. 


(5.2.19) 


Thus  the  Fermi  level  becomes 


Ef=  Ev  + KT  log 


K1 

L NaJ 


(5.2.20) 


From  this  equation,  we  find  that  the  Fermi  level  can  be  moved  across 
the  lower  half  of  energy  band  gap.  This  method  is  good  for 
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investigating  an  unknown  energy  level.  When  Na  « Nv  is  fulfilled,  it 
can  be  used  to  investigate  deep  level  defects. 

5.2. 1.4.  Single-Acceptor  Noise  Model 

When  the  device  is  lightly  doped  with  deep  single-acceptor 
defects,  we  obtain,  following  the  same  procedures  as  outlined  in 
section  5.2. 1.2, 


Equations  (5.2.15)  and  (5.2.22)  have  the  same  form,  and  as  a result  we 
are  unable  to  discriminate  between  acceptor-like  and  donor-like 
behavior  of  a defect  by  measuring  G-R  noise  in  a device  with  doping 
concentrations  satisfying  Ni«  Na.  However,  we  can  use  the  magnitude 
of  the  extracted  capture  coefficient  to  infer  the  donor  or  acceptor 
behavior  of  a defect  energy  level.  At  low  temperatures,  the  coefficient 
cp2  (see  Table  5.1)  is  associated  with  a repulsive  capture  process,  and 
one  may  expect  its  value  to  be  rather  small  [27],  of  the  order  of  10" 11 
to  10" 14  cm3/s;  the  coefficient  cpi  is  related  to  a neutral  capture 
process,  and  will  be  therefore  of  the  order  of  10"8  cm^/s;  the 
coefficient  cpo  on  the  other  hand  is  associated  with  an  attractive 
capture  process,  and  one  may  expect  its  value  to  be  rather  large,  of  the 
order  of  10"5  to  10"8  cm^/s. 

Now,  let's  consider  a device  that  is  doped  with  a deep  defect 
concentration  larger  than  a large  shallow  dopant  concentration.  When 
the  deep  defect  is  acceptor-like,  the  Fermi  level  will  behave,  at  lower 


with 


(5.2.21) 


(5.2.22) 
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temperatures,  according  to  eq.  (5.2.20)  provided  that  all  shallow 
acceptors  are  ionized.  With  Na  large,  the  temperature  variation  of  Ef 
will  be  small  making  it  difficult  for  Ef  to  reach  values  close  to  the  deep 
acceptor  level.  As  a consequence  the  hole  transition  rates  between  the 
deep  acceptor  level  and  the  valence  band  will  be  small,  providing 
negligible  noise  contributions  to  the  overall  noise  spectrum.  This 
makes  it  very  difficult  to  detect  deep  acceptor-like  traps  in  p-type 
material.  However,  when  the  deep  defect  is  donor-like,  the  Fermi 
level  is  above  the  shallow  acceptor  level  and  close  to  the  deep  donor 
level  (see  eq.  (5.2.18)),  and  a strong  noise  component  results. 
Consequently,  with  the  aid  of  a device  with  a deep  defect 
concentration  larger  than  the  shallow  dopant  concentration,  we  are 
able  to  distinguish  a donor  from  an  acceptor  by  noise  measurements. 


5.2.2.  Double-Donor  Noise  Model 

In  this  section,  we  deal  with  a double-donor  defect  of  which 
some  aspects  have  been  discussed  in  Chapter  III.  The  total  defect 
concentration  is  given  by 

Ni++  + Ni+  + Ni°  = Ni.  (5.3.1) 

From  eqs.  (5.1.2),  (5.1.3)  and  (5.3.1),  the  ionized  defect 


concentrations  can  be  written  as 

Ni 


Ni+  = 


i + — e 

®2d 


E2d"  Ef 
KT 


Sid  6 


IE” E7 
KT 


(5.3.2) 


and 
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Ni++  = 


Ni 

- (Eld+E2d-2  EjJ 


* + &2d  ^ 


KT 


+re 

&ld 


-te— E7  . 


(5.3.3) 


KT 


The  value  of  P can  be  derived  from  eqs.  (2.23)  and  (5.3.1)  to  (5.3.3) 
with  the  use  of  the  charge  neutrality  equation 


2Ni++  + Ni+  + P = Na. 

The  rate  equations  for  this  system  read: 
dNi° 

~d T~  = ' ri* 

dNi+ 

"dt  = - gi  + ri  + g2  - T2, 

dNi++ 

"dt  = " g2  + P2, 


dt*  = gl  ' rl  + g2  ‘ r2. 
where  gi,  ri,  g2,  r2  are  described  in  Table  5.1. 
From  eqs.  (5.3.1)  and  (5.3.4),  we  obtain 


(5.3.4) 


(5.3.5a) 

(5.3.5) 

(5.3.6) 

(5.3.7) 


ANi++  + ANi+  + ANi°  = 0,  (5.3.8) 

2ANi++  + ANi+  + AP  = 0.  (5.3.9) 

Expanding  the  carrier  densities  around  steady-state  in  first  order,  and 
using  the  Langevin  approach,  the  rate  equations  including  the  noise 
can  be  written  as 


dAP 

dt 

dANi++ 

dt 


= -Mu  AP  - M12  ANi++  +Fi  + F2, 


= -M21  AP  - M22  ANi++  - F2, 


with  relaxation  matrix  elements 


Mn=  epi-cp2  Po+Cpi  Nio°+Cp2  Nio++cpi  P0, 
M22=  6p2  + 2cp2  PQ. 


(5.3.10) 

(5.3.11) 
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Mi2=  - ep2+cpi  P0+2epi-2cP2  P0, 

M2i=  - cp2  Nio+  + cp2  P0.  (5.3. 12) 

The  Langevin  noise  source  terms  Fj  and  F2  model  the  noise  in  the 
hole  transition  rates  between  the  first-donor  level,  the  second-donor 
level  and  the  valence  band  respectively 

Using  the  same  procedures  as  outlined  in  Chapter  III,  and  assuming 

that  there  is  no  correlation  between  the  Langevin  noise  sources  and 

that  the  transition  rates  show  full  shot  noise,  we  obtain  from  eqs. 

(5.3.10)  through  (5.3.12) 

4ri  4r2 

^F1_  L w d ’ Sf2-  l w d ’ 

„ , , G H 

V®'-  (©2Ui2)  + (COW)  • 

_ /. j2(r  1 + I"2)-M222  (ri+r2)-(Mi22+2Mi2M22)  r2 

L w d (Ki'-'ki) 

TT  „ _^22(rl+r2)+M222  (ri+r2)  + (Mi22+2Mi2M22)  P2 

H = 4 , “ _ 

L w d (X,i2-^22) 

with 

^ (Mii+M22)+V  (Ml  i+M22)2'4(MiiM22"Mi2M2l) 

Ai  = 2 

and 

^ (Mii+M22)_V  (Ml  l+M22)2"4(MiiM22"Mi2M2l) 

A.2  = 2 

5.2. 2.1.  Deep  Defects  with  Closely  Spaced  Energy  Levels  in  Nearly- 
compensated  Material 

In  this  subsection,  a device  with  a double-donor  defect 
concentration  satisfying  Na  > Ni  > 0.5Na  is  studied.  At  low 

+ ++  O 

temperatures  we  may  assume  P « Ni  , Ni  , Ni  . This  assumption  is 
correct  as  long  as  the  defect  energy  levels  are  closely  spaced,  i.e.,  not 


(5.3.13) 

(5.3.14) 

(5.3.15a) 

(5.3.15b) 

(5.3.16) 

(5.3.17) 
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more  than  a few  KT  apart.  The  relaxation  matrix  elements  (5.3.12) 
become 

Mii=cpi  Ni°+cp2  Ni+,  Mi2=  ~eP2+CpiP+2epi-2cp2  P0. 

M2i=  -cP2  Ni+,  M22  =Cp2  + 2cP2  P0.  (5.3.18) 

Note  that  eq.  (5.3.14)  predicts  two  Lorentzian  components  in  the 
overall  noise  spectrum.  However,  if  A,i  is  too  small,  only  Spp(CO)  for 
CO  » A*2  will  be  observed. 

If  cp2  « cpi,  we  may  assume 


Cpl 


» 


Ni 


.+ 


1 


Eld~  Ef 

KT 


cp2  Ni°  gld  ~ 

This  assumption  is  valid  if  (Eid-E2d)  is  smaller  than  a few  KT.  Then  Ef 
is  close  to  Eid  and  E2d,  and  Ni+  and  Ni°  are  of  the  same  order. 

Using  eqs.  (5.1.8),  (5.1.9),  (5.3.16),  (5.3.17)  and  (5.3.19),  we  obtain 

Mu  » M21,  M22,Mi2, 

» X-2 , 


(5.3.19) 


= Mu  =cpi  Ni°  - 7 , 

Lp 


Xo  = 


(MnM22‘Mi2M2l) 


< M22  =Cp2  + 2cp2  P. 


(Mi  1+M22) 

Using  eqs.  (5.1.8),  (5.1.9),  and  (5.3.20),  we  get 

q 4(ri  +r2) 

pp  - L w d (CO2  +^i2)  ’ 
ri+  r2=  (cpi  Nio°  +cp2  Nio+)P0  = MnP0, 


(5.3.20) 


toj  _ Spp((Q)  4 Cpl  Mo0  pQ  tp2 
V0^-  “ Po2  " Lwd  P02  (1+CO2  tp2) 

Equation  (5.3.21)  can  be  approximated  by 

4 


for  A,2«  CO.  (5.3.21) 


SJCO) 

“Vo2- 


LwdP0  Cpi  Nio° 


for  A,2«  co  « Xi. 


(5.3.22) 
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Equations  (3.7),  (5.3.1)  to  (5.3.4)  and  (5.3.22)  can  be  used  to  extract 
trap  parameters  from  the  resistance  data  and  noise  data  if  we  use  the 
mobility  formula  (2.19). 

The  theoretical  predictions  for  the  noise  and  resistance  can  be 
further  simplified  by  assuming  (Eld  - Ef)  > KT  log(gid)  which  is 
certainly  valid  at  low  temperatures.  Equation  (5.3.3)  can  be 
approximated  by 


Ni++  = 


Ni 


1 + rle 

&ld 


KT 


Using  eqs.  (5.3.1)  and  (5.3.4),  neglecting  P in  eq.  (5.3.4)  and  assuming 
Ni++,  Ni+  > Ni°,  the  hole  concentration  at  low  temperatures  is  given  by 


e - En , 

v 2d 


P = 


^2d  (M  XT-  P KT 

— (Na  - Ni ) Nv  e 


2 Ni-Na 

Using  the  following  approximations, 


(5.3.23) 


Ni++  = Na-Ni, 

T-  + 


(5.3.23a) 

NiT  = 2Ni-Na,  (5.3.23b) 

and  substituting  eqs.  (5.1.2)  and  (5.3.23)  into  eq.  (5.3.22),  the 
expression  for  the  normalized  G-R  noise  becomes 


S (CO) 

72  = 


4 e 


E.  E 
Id  v 

i rr~ 


for  \q<<  (0  « A«i.  (5.3.24) 


v L w d cpi  gld(2Ni-Na  ) Nv 

Equations  (5.3.23)  and  (5.3.24)  illustrate  that  the  noise  contains  the 
activation  energy  of  the  first-donor  level,  whereas  the  carrier 
concentration  or  resistance  shows  the  activation  energy  associated 
with  the  second-donor  level. 
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Equation  (5.3.22)  is  the  same  as  eq.  (5.2.9).  This  shows  that 
when  the  doping  concentrations  satisfy  Na«Ni,  and  the  hole  density 
can  be  approximated  by  eq.  (5.2.10),  the  double-donor  model  reduces 
to  the  single-donor  model. 


5. 2.2.2.  Deep  Defects  with  Closely  Spaced  Energy  Levels  in 
Uncompensated  Material 

If  the  defect  concentration  is  smaller  than  the  boron 
concentration,  we  may  assume  that  P = Na,  P » Ni,  Ni°,  Ni+,  Ni++ 
hold.  In  the  following,  we  limit  our  discussion  to  the  conditions: 
cpi  » cp2  as  assumed  in  section  5.2.2. 1,  and  epi  » ep2,  from  which  we 
derive,  with  the  use  of  eqs.  (5.1.8)  and  (5.1.9), 


e 

CP2  gld 


Eld  " E2d 
KT 


The  relaxation  matrix  elements  (5.3.12)  become 


Mn=  epi  + cpi  Po,  M2i=  eP2  + 2cP2  Po, 

Mi2=  cpiPo+2epi,  M22  = cP2  Po.  (5.3.25) 

In  this  section,  we  discuss  the  situation  where  the  defect  energy  levels 
are  close  to  each  other.  With  the  use  of  eqs.  (5.1.2),  and  (5.1.3),  we 
obtain  that  Ni°,  Ni+,  Ni++  are  almost  of  the  same  order.  With  the  use  of 
eqs.  (5.1.8),  and  (5.1.9),  we  infer  that  epi,  cpi  P0,  and  eP2,  cP2  P0  are 
almost  of  the  same  order,  respectively.  Using  these  results  and  eqs. 
(5.3.16),  (5.3.17)  and  (5.3.25),  we  obtain 


Mu,  M12  » M22.  M21, 

A-i  = Mu  = epi  + cpi  P0  = — , 


X.2  = 


(Mi  iM22~Mi2M2l) 
(M11+M22) 
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6p2  Cpi  Po  +2  6pi  ep2  + Cpi  Cp2  Pp  +3epi  Cp2  Pp 

epl  + cpl  Po+  Cp2  P0 

Since  Mp  is  of  the  same  order  as  epi,  X2  is  of  the  same  order  of  -eP2, 
the  observed  noise  spectra  will  have  only  one  Lorentzian  component  in 
a certain  frequency  window  because  -X2  is  smaller  than  and  the 
same  reason  used  in  section  5.2.2. 1 applies  to  this  case.  The  roll-off 
time  constant  is  the  same  as  given  by  eq.  (5.2.15)  of  the  single-donor 
model. 

From  eq.  (5.3.19)  and  Table  5.1  we  get  ri»  r2.  Equation  (5.3.14) 
can  be  approximated  by 

SPP,m)  = L w for  -*-2<<co'  (5-3-261 

ri  = cpi  Nio°  P0, 


S^(CO)  Spp(CO)  4 c Nio”  P0  tpz  „ , ..  

tt  - LWdVu'^T  for'X2<<“'  ,5-3'26a) 

Using  ri  = gi,  eq.  (5.3.26a)  transforms  to 

Sw(°)  SPP(°)  4 Nio+  ePi  tp  r ^ ^ r-  n --1 

VV7  = t^o2" = LwdpJ  f°r  ‘ l2<<  “ <<:  Xl'  |5-3'27) 

Note  that  eqs.  (5.2.17)  and  (5.3.27)  are  of  the  same  form  with  the 
same  lifetime  constant.  In  this  case,  we  are  unable  to  measure  the 
second-donor  level  of  defects. 


r.  O 


VV 

V, 


5. 2.2. 3.  Deep  Defects  with  Energy  Levels  Far  Apart  from  Each  Other 
in  Overcompensated  Material 

For  an  overcompensated  device,  with  Ni  >Na,  but  of  the  same 
order,  and  with  the  energy  levels  spaced  by  more  than  a few  KT,  we 
may  assume  at  low  temperatures  that  Na  = Ni+,  and  Ni  - Na  = Ni° 
using  the  same  reasoning  as  described  in  5.2. 1.1.  The  Fermi  level  is 
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close  to  the  first-donor  level  and  Nio++,  Po  « Nio+,  Na  holds.  Using 
cpi  »cp2  as  assumed  in  section  5.2.2. 1,  and  assuming  that  epi«  ep2, 
which  is  possible  because  the  conditions  Eid-E2d  > a few  KT  and 
cpi  »cp2  hold,  eq.  (5.3.12)  becomes 

Mn=  cpi  Nio°,  M22=  ep2. 


Mi2=  -eP2,  M2i  = -cP2  Nio+. 

If  cpi  Nio°  » eP2  holds,  i.e.. 


(5.3.29) 


(Ni  -Na)  » ( 7*7}  N„  ( -^  e 

cpi  v gid 


&2d. 


Ev‘E2d 


KT 


we  obtain  with  the  use  of  eqs.  (5.3.14)  through  (5.3.17)  and  (5.3.29) 

Mu  » M21* 


Mn»  -Mi2=  M22. 

A,i  = Mu  = cpi  Nio°  = t — » X2, 

Ln 


Sw(“) 


Spp«0) 


v -2  P : 

vo 


2" 


4 cPi  Nio°  P0  tp2 
L w d P02  (1+CO2  tp2) 


for  X,2«  CO. 


(5.3.30) 


Equation  (5.3.30)  can  be  approximated  by 

4 


Sw^ 
7T~  = 


for  X,2«  (0  « X-i,  (5.3.31) 


VQZ  L w d P0  cpi  Nio° 
which  is  the  same  as  given  by  eq.  (5.2.9)  of  the  single-donor  model. 

If  the  hole  concentration  P is  approximated  by  eq.  (5.2.10),  the 
double-donor  reduces  to  the  single-donor  model  discussed  in  5.2. 1.1. 

For  a device  with  a defect  concentration  satisfying 
Na>  Ni  > 0.5Na,  we  may  assume  cpi/cp2  « Ni+/Ni°  because  the  Fermi 
level  is  close  to  the  second-donor  level  and  far  away  from  the  first- 
donor  level. 

We  start  with  the  rate  equations  including  the  noise 
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dAP 

dt 


= -Mu  AP  - M12  ANi°  + Fi  + F2, 


(5.3.32) 


dANi° 

dt 


= -M2i  AP  - M22  ANi°  + Fi, 


(5.3.33) 


with  relaxation  matrix  elements 

Mu  = -epi+cp2  P0+cp2  Nio++cpi  Nio°+ep2, 

Mi2  = 2epi-2cp2  PG  + Cpi  PD  - Cp2» 

M2i  = cpi  Nio°  - epi  , 

M22  = 2epi  + cpi  PG.  (5.3.34) 

Using  the  condition  PQ  « Nj0+,  eq.  (5.1.8)  and  similar  approximation 
procedures  as  outlined  in  section  5.2.2. 1,  we  get 


sw<<°>  = 4 

V02  L w d P0  cp2  Nio+ 


for  X2«  CO  « A,i, 


(5.3.35) 


which  is  in  the  same  form  as  given  by  eq.  (5.2.9)  of  the  single-donor 
model. 

Using  eq.  (5.3.23b),  eq.  (5.3.35)  becomes 

4 


Swl“) 
~2~  = 


V, 


o 


LwdP0  cp2  [2Ni-Na] 


for  X2«  CO  « X,i.  (5.3.36) 


Under  the  conditions  studied  here,  the  normalized  G-R  noise  is 
inversely  proportional  to  the  hole  concentration  P,  indicating  the 
noise  and  resistance  show  the  same  donor  activation  energy. 


5. 2.2. 4.  Deep  Defects  with  Energy  Levels  Far  Apart  from  Each  Other 
in  Uncompensated  Material 

When  the  Fermi  level  is  close  to  the  second-donor  level  and  far 
away  from  the  first-donor  level  in  a situation  where  the  levels  are 
spaced  by  more  than  a few  KT,  we  can  assume  Sf2  » Spi  which  can  be 
transformed  to  ep2  Ni++  » epi  Ni+  with  the  use  of  eq.  (5.3.13)  and 
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Table  5.1.  From  eqs.  (5.1.2),  (5.1.3),  (5.1.8),  (5.1.9)  and  (5.3.13),  this 
condition  translates  into 


Eld"  Ef 

T-e  KT 

Cp2  gld 


(5.3.37) 


The  condition  Sf2  » Sfi  implies  F2  » Fi.  The  carrier  fluctuations  will 
then  mainly  come  from  the  fluctuations  in  the  emission  and  capture 
rates  associated  with  the  second-donor  level,  and  therefore  eq.  (5.3.9) 
becomes 

-ANi+  = - AP  = ANi++.  (5.3.38) 

The  rate  equations  including  the  noise  can  then  be  written  as 
d AP 


= -MnAP  - M12  ANi++  + F2, 


dt 

d ANi++ 
dt 


= -M21AP  - M22  ANi++  - F2, 


(5.3.39) 

(5.3.40) 


with  relaxation  matrix  elements 

Mu  = epi-cp2  Po+Cpi  Nio°+cP2  Nio++cpi  PQ, 

M21  = eP2  + 2cP2  P0, 

M12  = -eP2+CpiPo+2epi-2cp2Po. 

M22  =-  cp2  Nio+  + cp2  P0.  (5.3.41) 

With  additional  conditions  epi/eP2  «1,  and  P » Ni,  Ni°,  Ni+,  Ni++,  and 
with  the  aid  of  eq.  (5.1.9),  one  can  show  that  eqs.  (5.3.39)  and  (5.3.40) 
have  the  same  time  constant  defined  by 

■j”  = ep2  + cp2P0.  (5.3.42) 

Under  the  above  stated  conditions,  the  double-donor  model  reduces  to 
the  single-donor  model  discussed  earlier.  However,  this  time  we  can 
only  measure  the  activation  energy  of  the  second-donor  level. 
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When  the  Fermi  level  is  close  to  the  first-donor  level  and  far 
away  from  the  second-donor  level,  we  can  assume  Sp2  « Spi,  which 


becomes,  with  the  use  of  eqs.  (5.1.2),  (5.1.3),  (5.1.8),  (5.1.9),  (5.3.13) 
and  Table  5.1, 


Eld'  Ef 


£ei»  J_e  KT 

Cp2  gld 


(5.3.43) 

The  condition  Sp2  « Spi  implies  F2  « Fi . The  carrier  fluctuations  will 


come  mainly  from  first-donor  level  transitions,  and  therefore  eq. 
(5.3.9)  becomes 

ANi+  = - AP  = -ANi°. 


The  rate  equations  including  the  noise  can  then  be  written  as 


d AP 
dt 


= -Mu  AP  - Mi2  ANi°  + Fi, 


d ANi° 
dt 


= -M21  AP  - M22  ANi°  + Fi, 


(5.3.44) 

(5.3.45) 


with  relaxation  matrix  elements  given  by  (5.3.34). 

With  additional  conditions  cP2/cpi  «1  and  P » Ni,  Ni°,  Ni+,  Ni++,  eqs. 
(5.3.44)  and  (5.3.45)  have  the  same  time  constant  defined  by 

= epi  + cpi  P0.  (5.3.46) 

Again,  the  double-donor  model  reduces  to  the  single-donor  model 
discussed  earlier.  In  this  case  however  we  can  only  measure  the  G-R 
noise  caused  by  the  first-donor  level. 


5.2.3.  Acceptor-Donor  Noise  Model 

In  this  section,  we  calculate  the  resistance  and  the  voltage  noise 
resulting  from  deep  defects  with  an  acceptor  and  a donor  level  that 
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interacts  with  the  valence  band.  The  total  deep  defect  concentration  is 
given  by 


Ni’  + Ni+  + Ni°  = Ni.  (5.4.1) 

From  eqs.  (5.1.1),  (5.1.2)  and  (5.4.1),  the  ionized  deep  defect 
concentrations  are  rewritten  as 

Mi 

(5.4.2) 


Ni+  = 


Ni 


1 + Sid  &la  ^ 


2E..-E,  -E.. 
f la  Id 

KT 


Sid  e 


KT 


and 


Ni'  = 


nq-Ej- 

1 -hr-e  KT 


Ni 


®ld+Ela'2  Ef* 


(5.4.3) 


KT 


Sla  Sia  Sid 

The  value  of  P can  be  derived  from  eqs.  (5.4.1)  to  (5.4.3)  and  eq.  (2.23) 
with  the  use  of  the  charge  neutrality  equation 

Ni+  + P = Na  + Ni'.  (5.4.4) 

The  rate  equations  for  the  carrier  densities  read: 


dNi+ 

“dt”  =-gi+ri’ 
dNi' 

dt  - So  - rG, 
dNi° 

~fa — = - go  + r0+ gi  - ri, 

dP 

(jj;  - Si  _ ri  + So  _ ro. 
where  gi,  ri,  go,  rD  are  described  in  Table  5.1. 
From  eqs.  (5.4.1)  and  (5.4.4),  we  obtain 


(5.4.5) 


(5.4.5a) 

(5.4.6) 

(5.4.7) 


ANi'  + ANi+  + ANi°  = 0, 

and 


(5.4.8) 


-ANi'  + ANi+  + AP  = 0. 


(5.4.9) 
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Expanding  the  densities  around  steady-state  and  using  the  Langevin 
approach,  the  rate  equations  including  the  noise  can  be  written  as 


dAP 

dt 


= -MnAP  - M12  ANi+  + F0  + Fi, 


dANi+ 

dt 


= -M2iAP  - M22  ANi+  - Fi, 


(5.4.10) 

(5.4.11) 


with  relaxation  matrix  elements 

Mu  = epo-Cpi  Po+Cpo  Nio  +cpi  Nio°+Cpo  P0, 
M21  = epi+  2cpi  P0, 


M12  — _6pi+Cpo  P0+2ePo-2cpi  Po, 

M22  = -Cpi  Nio°  + Cpi  Po,  (5.4.12) 

where  F0 , Fi  are  Langevin  noise  source  terms  modeling  the  noise  in  the 
hole  transition  rates  between  the  acceptor  level,  the  first-donor  level, 
and  the  valence  band  respectively. 

Assuming  that  there  is  no  correlation  between  the  Langevin  noise 
sources  and  that  the  transition  rates  show  full  shot  noise,  we  obtain 

4ri 


S 4r° 
&Fo_  L w d 


Sfi  = 


L w d 


(5.4.13) 


„ , , G H 

Spplt0'-  (co2+7i,2)  + (coW)  ' 

\i2(ri+r0l-M222  (r]+r0)-(Mi22+2M12M22)  r„ 

L w d (ki'-'ki) 

tt  „ -A,22(ri+r0)+M222  (ri+r0)+(Mi22+2Mi2M22)  rQ 

H = 4 . 0 ^ 

L w d (Xi2-X-22) 


(5.4.14) 

(5.4.15) 

(5.4.16) 


with 


^ (Mn+M22)+V (Mi  i+M22)2-4(Mi  1M22-M12M21)  fr.  A , ^ 

Ai  = 2 » (5.4.17) 


and 
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^ (Mh+M22)W  (Mii+M22)^-4(MiiM22-Mi2M2l) 

M= 2 • (5.4.18) 


5.2.3. 1.  Deep  Defects  with  Energy  Levels  Close  to  Each  Other  in 
Overcompensated  Devices 

+ - o 

At  low  temperatures  we  may  assume  P«  Ni  , Ni  , Ni  in  an 
overcompensated  device  with  Ni  > Na.  The  assumption  is  valid  as  long 
as  the  energy  levels  of  the  deep  defects  are  within  a few  KT  of  each 
other.  With  the  aid  of  eq.  (5.1.7),  the  relaxation  matrix  elements 
(5.4.12)  become 


Mn-Cpo  Nj  +Cpi  Ni  , M12  - -Cpi+Cpo  Po+2epo_2cpi  Po, 

M2i  =epi  + 2cpi  P0,  M22=  - cpi  Ni°.  (5.4. 19) 


In  the  following,  we  limit  our  discussion  to  the  situation  where 
cpo  » cpi  and 


-po 


» 


Ni 


•o 


1 


ElaEf 

KT 


(5.4.20) 


cpi  ' ' Nf  gla 
The  condition  cpo /cpi  »Ni°/Nf  is  obtained  when  the  acceptor  and 
donor  levels  are  very  close  to  each  other  and  the  Fermi  level  is  close 
to  both  Eid  and  Eia. 

Using  eqs.  (5.1.7),  (5.1.8),  (5.4.17),  (5.4.18)  and  (5.4.20),  we  obtain 


Mu  » M21,  M22.  M12, 

(5.4.21) 

A,i  » X2 . 

(5.4.22) 

Xi  = Mi  1 = Cpo  Nio'  = 7-  , 

ip 

(5.4.23) 

. (Mi  1M22-M12M21)  . 

A*=  (M11+M22)  ^ 

15.4.24} 

Using  eqs.  (5.4.14)  through  (5.4.16)  and  (5.4.23),  we  get 

S (CO)  = - — 1 +2°l  v for  %2«  0),  (5.4.25) 

pp  L w d (CO2  +>.i2) 
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ri+  r2  = (cpi  Nio°  +Cpo  Nio')P0  = MhPq, 


(5.4.26) 


S^tCQ)  Spp((Q)  4 cpo  Nio'  PQ  tp: 


V 


~r~  = 


P? LwdP  “(lWtp4)  f0t  ^2<<  “•  (5-4'27) 


O ‘ O * Q t.p 

Equation  (5.4.27)  can  be  approximated  by 

4 


Sw«°> 


V7"  LwdPQ  Cpo  Nio' 


for  X,2«  CO  « A,i. 


(5.4.28) 


Equations  (3.7),  (5.4.1)  to  (5.4.4)  and  (5.4.28)  can  be  used  to  extract 
trap  parameters  from  the  resistance  data  and  noise  data  if  we  use  the 
semi-empirical  mobility  formula  (2.19). 

The  theoretical  predictions  for  the  noise  and  resistance  can  be 
further  simplified  if  KT  log(gia  cpo/cpi)  > (Ela  - Ef)  > KT  log(gia)  and 
(Eid-Ef)  > KT  log(gid).  Then  eq.  (5.4.2)  can  be  approximated  by 

Ni 


Ni+  = 


(5.4.29) 


KT 


1 + gld  e 
Using  the  following  approximations, 

Ni'  < Ni,  Ni+,  Ni°, 

Ni+  = Na, 

Ni°  = Ni  - Na,  (5.4.30) 

and  eqs.  (2.23),  (5.4.29)  and  (5.4.30),  the  hole  concentration  at  low 
temperatures  is  given  by 


P = 


Sid  Na  Nv  e 


E - E,  , 
v Id 

KT 


(5.4.31) 


Ni-Na 

Using  eqs.  (5.1.1)  and  (5.4.31)  into  eq.  (5.4.28),  the  expression  for  the 
normalized  G-R  noise  becomes 
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V 2 ~ 


4 e 


E,  -E 
la  v 

KT 


for  Xi«  (0  « A,i.  (5.4.32) 


o L w d cpo  gla(Ni-Na  ) Ny 

Equations  (5.4.31)  and  (5.4.32)  illustrate  that  the  noise  contains  the 
activation  energy  of  the  acceptor  level,  whereas  the  carrier 
concentration  or  resistance  shows  the  activation  energy  associated 
with  the  the  donor  level. 


5. 2. 3. 2.  Deep  Defects  with  Energy  Levels  Far  Away  from  Each  Other  in 
Overcompensated  Material 

Instead  of  eq.  (5.4.20)  we  limit  our  discussion  in  this  section  to 
the  condition  cpo/cpi  « Ni°/Ni\  which  prevails  when  the  acceptor 
and  donor  levels  are  separated  by  more  than  a few  KT,  such  that  the 
Fermi  level  is  close  to  Eia  and  far  away  from  Eia.  We  find,  with  the  use 
of  eqs.  (5.1.1)  and  (5.1.2),  that  Ni°»Nf  . 

We  start  with  the  rate  equations  including  the  noise 

= - Mu  AP  - M12  ANi'  + F0  + Flt  (5.4.33) 

dANf 

— = - M2i  AP  - M22  ANi'  + F0,  (5.4.34) 

where 

Mn=  cpi  Nio°  +Cpo  Ni0'+epi  +cpi  Po  " Cpo  > M2i=  -epo+  Cpo  Nio', 

Mi2—  -epi+Cpo  P0+2epo  -2cpiP0,  M22  = Cpo  P + 2ePo*  (5.4.35) 
At  low  temperatures  PQ  is  small  in  overcompensated  devices, 
and  if  Nio°  is  large,  we  may  assume  PQ  « Nio°.  Using  eqs.  (5.1.8), 
(5.1.7)  and  similar  approximation  procedures  as  outlined  in  section 
5.2.3. 1,  we  get 

Sw^)  = 4 

V0^  LwdP0  cpi  Ni° 


for  X,2«  (0  « A-i.  (5.4.36) 
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This  is  of  the  same  form  as  eq.  (5.2.9).  Using  eq.  (5.4.30),  eq.  (5.4.36) 
becomes 


Sw(“) 


V0 


7—  = 


LwdP0  cpi  (Ni-Na) 


for  X,2«  (0  « A,i.  (5.4.37) 


The  hole  density  P can  be  approximated  by  eq.  (5.2.10).  The  acceptor- 
donor  model  reduces  to  the  single-donor  model  discussed  in  l.A.  The 
results  are  straightforward,  because  when  the  Fermi  level  is  close  to 
the  donor  level,  the  negatively  charged  defect  concentration  is 
negligible,  and  the  contribution  of  acceptor  levels  to  the  G-R  noise  can 
be  neglected.  In  this  case,  we  are  unable  to  measure  the  properties  of 
the  acceptor  level  using  G-R  noise  measurements. 


5. 2.3. 3.  Deep  Defects  with  Energy  Levels  Closely  Spaced  to  Each 
Other  in  Uncompensated  Material 

If  the  defect  concentration  is  smaller  than  the  shallow  impurity 

concentration  and  the  energy  levels  of  a defect  are  very  close,  then  we 

may  assume  P=  Na,  P » Ni,  Ni°,  Ni+,  Ni\  The  relaxation  matrix 

elements  (5.4.12)  become 

Mn=  Cpo-Cpi  Po+Cpo  Po»  M2i=  epi  + 2cpi  Po, 

Mj2=  -epi+CpoPo+2epo-2Cpi  Po,  M22=  Cpi  Po.  (5.4.38) 

Using  cpo  » cpi  as  assumed  in  section  5.2.3. 1,  the  relaxation  matrix 
elements  (5.4.38)  become 

Mn=  epo  + Cpo  P0,  M2i=  epi  + 2cpi  P0, 

Mi2=  Cpo  P+2ep0,  M22  = Cpi  Po-  (5.4.39) 

Using  eqs.  (5.4.17),  (5.4.18),  (5.4.20)  and  (5.4.39),  we  obtain 
Mu,  M12  » M22.  M21, 
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A<i  _ Mu  — 6po  + Cpo  Po  — f > 

Lp 

. __  (MiiM22_Mi2M2l) 

A2=  (Mii+M22) 

_ 6pi  Cpo  Pp  +2  ep0  epi  + Cpp  Cpi  Pp2  +3epp  Cpi  Pp 
Cpo  + Cpo  PG+  Cpi  PQ 

Since  Mu  is  of  the  same  order  as  Mi2,  -X2  is  of  the  same  order  as  M22, 
and  the  observed  noise  spectra  will  have  only  one  Lorentzian 
component  in  a certain  frequency  window  because  -X2  is  smaller  than 
A,i.  The  same  reasoning  as  used  in  section  3.1  applies  to  this  case.  The 
roll-off  time  constant  is  the  same  as  given  by  eq.  (5.2.22)  for  the 
single-acceptor  model. 

From  eq.  (5.4.20)  and  Table  5.1  we  get  ri»  r2.  Equation  (5.4.14) 
can  be  approximated  by 

4r 

S (C0)=  - , , °i  *1  £■>  for  -A,2«  (0,  (5.4.40) 

pp  Lwd  (0)24f) 

ri  = Cpo  Nio'  P0, 


Sv>>  _ 4 Cpp  Nio'  Pp  tp2 

VQ2  ~ Po2  "LwdPo2  (1+CO2 tp2) 

Using  r0  = go,  eq.  (5.4.41)  transforms  to 


for-^2«W-  (5.4.41) 


Sw(°)  _ Spp(Q)  4 Nio°  ePo  tp 
tpVo2  tpPo2  LwdPo2 


for  - \2«  CO  « X-i. 


(5.4.42) 


Note  that  eqs.  (5.4.42)  and  (5.2.17)  are  of  similar  form.  In  this  case. 


we  are  unable  to  measure  the  donor  level  of  defects. 


5. 2.3. 4.  Deep  Defects  with  Energy  Levels  Separated  by  More  than  a 
Few  KT  from  Each  Other  in  Uncompensated  Material 

+ - ° 

At  high  temperatures  we  may  assume  P » Ni  , Ni  , Ni  in 
devices  with  defect  concentrations  satisfying  Na  > Ni.  When  the  Fermi 
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level  is  close  to  the  acceptor  level  and  far  away  from  the  donor  level, 
we  may  assume  SFo  »Spi  which  implies  F0  » Fj . The  rate  equations 
including  the  noise  can  then  be  written  as 


dAP 

dt 

= -Mu  AP  - Mj2  ANi  + F0, 

(5.4.43) 

dANi" 

dt 

= - M21  AP  - M22  ANi  + F0, 

(5.4.44) 

with  the  relaxation  elements  given  by  eq.  (5.4.35). 

With  eqs.  (5.1.2),  (5.1.7),  (5.1.8),  (5.4.13)  and  Table  5.1,  SFo  »SFi 
results  in 


£0 

Cpl 


Ela  Ef 


ITT 


(5.4.45) 


The  carrier  fluctuations  will  come  mainly  from  the  acceptor  level, 
therefore  eq.  (5.4.9)  becomes 

ANf  = AP.  (5.4.46) 


With  an  additional  condition  cpo/cpi  »1,  eqs.  (5.4.43)  and  (5.4.44) 
have  the  same  time  constant  defined  by 


^ — Mu  — epo  + cpo  Po- 


(5.4.47) 


Consequently,  the  acceptor-donor  model  reduces  to  the  single- 
acceptor model  discussed  earlier. 

When  the  Fermi  level  is  close  to  the  donor  level  and  far  away 
from  the  acceptor  level,  we  may  assume  SFo  « SFi,  which  becomes 
with  eqs.  (5.1.2),  (5.1.7),  (5.1.8),  (5.4.13)  and  Table  5.1 


E,  - Er 
la  f 

Cpo  1 _ KT 

~~ — « - — e 

Cpl  gla 


(5.4.48) 


Equations  (5.4.10)  and  (5.4.11)  become 
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dAP 

dt 


= -Mu  AP  - M12  ANi+  + Fi, 


(5.4.49) 


dANi+ 

dt 


= -M2i  AP  - M22  ANi+ 


Fi, 


(5.4.50) 


with  relaxation  matrix  elements 

Mu  = Cpo'Cpi  Po+Cpo  Nio  "*'Cpi  Nio  +Cpo  Po. 

Mi2  = -Cpi+Cpo  Po+2epo-2cpi  Po, 

M2i  = epi  + 2cpi  P0. 

M22  = -cpi  Nio  + Cpi  Po.  (5.4.51) 

The  carrier  fluctuations  will  come  mainly  from  the  donor  level, 
therefore  eq.  (5.4.9)  becomes 


ANj+  = - AP.  (5.4.52) 

Since  we  discuss  the  condition  (Eia  - Eid)  larger  than  a few  KT,  we 
may  assume  epo/epi  «1,  i.e., 

E1  - E,  i 
la  Id 

Cpo  gld  r KT 

cpl  gla 


with  which  we  find  that  eqs.  (5.4.49)  and  (5.4.50)  have  the  same  time 
constant  defined  by 

^ = 6pi  + Cpi  Pq.  (5.4.53) 

Hence,  the  acceptor-donor  model  reduces  to  the  single-donor  model 
discussed  earlier. 


5.3.  Conclusions 

In  the  design  of  a device  for  investigating  defect  properties  from 
G-R  noise  spectra,  the  following  facts  should  be  considered: 

1.  Most  of  the  energy  levels  of  defects  can  be  detected  in 
uncompensated  devices.  Under  certain  conditions,  G-R  noise  spectra 
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of  double-level  defects  reduce  to  G-R  noise  spectra  of  single-level 
defects.  However,  in  uncompensated  devices,  the  noise  signal  is  rather 
small  compared  to  that  of  overcompensated  devices. 

2.  In  order  to  measure  a particular  energy  level  of  a defect  in 
devices  doped  with  an  intended  deep  defect  and  a shallow  dopant,  the 
doping  impurity  concentrations  should  satisfy  the  following  conditions: 
Na  > Ni  > 0.5Na  for  a second-donor  level,  and  Ni  > Na  for  a first-donor 
level.  We  are  able  to  measure  the  capture  coefficients  as  a function  of 
temperature  and  electric  field  if  the  values  of  the  device  resistance  are 
such  that  the  cut-off  frequency  of  the  circuit  is  larger  than  the  roll-off 
frequency  of  the  G-R  noise  spectra. 

3.  If  the  energy  levels  of  a double-level  defect  are  very  close  and 
there  is  a very  large  difference  in  the  magnitude  of  capture 
coefficients  associated  with  each  level,  the  noise  and  resistance  data  in 
a nearly  compensated  device  are  likely  to  show  different  activation 
energies  associated  with  the  different  levels  of  the  same  defect.  It  is 
only  possible  to  detect  the  energy  level  that  has  the  larger  capture 
coefficient  in  an  uncompensated  device. 

4.  For  a defect  with  an  acceptor  and  a donor  level,  or  double  or 
triple  acceptor  levels,  it  is  difficult  to  detect  the  acceptor  level  in  a p- 
type  device  doped  with  an  intended  defect  concentration  larger  than 
the  shallow  dopant  concentration,  because  the  Fermi  level  is  far  away 
from  the  deep  acceptor  energy  levels. 

5.  For  an  uncompensated  device,  defect  activation  energies  can 
be  extracted  from  the  roll-off  time  constants  of  G-R  noise  spectra 
measured  as  a function  of  temperature.  In  a p-type  device  doped  with 
an  intended  defect  concentration  larger  than  the  shallow  dopant 
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concentration,  defect  activation  energies  are  extracted  from  noise 
plateau  and/or  resistance  data. 

6.  For  a defect  with  multi-energy  levels,  which  are  separated  by 
more  than  a few  KT  from  each  other,  the  multi-level  G-R  noise  model 
reduces  to  the  single-level  G-R  noise  model  in  both  uncompensated 
and  overcompensated  devices. 

7.  For  an  uncompensated  device,  eqs.  (2.15)  and  (2.23)  can  be 
used  to  extract  the  defect  activation  energy.  We  can  not  determine 
whether  the  activation  energy  is  associated  with  an  acceptor  or  a 
donor  center.  This  can  only  be  deduced  from  the  magnitude  of  the 
carrier  capture  coefficient  associated  with  the  defect.  Using  a device 
doped  with  an  intended  deep  defect  concentration  larger  than  the 
shallow  dopant  concentration,  the  electrical  nature  of  a defect  can  be 
determined  since  a deep  acceptor  center  is  unlikely  to  be  detected  in 
a p-type  device,  whereas  a deep  donor  level  will  contribute  a strong 
noise  component. 

8.  The  methods  presented  in  this  chapter  can  be  extended  to 
design  a device  for  investigating  a defect  with  more  than  two  energy 
levels  in  the  band  gap. 


CHAPTER  VI 

CARRIER  CAPTURE,  EMISSION  AND  RECOMBINATION 

6.1.  Carrier  Capture  and  Emission 
In  a well-known  paper  [36],  Lax  proposed  a "cascade  capture" 
model:  the  electron  is  captured  by  one  of  the  highly  excited 
recombination  center  levels,  which  are  quasi-continuously  arranged, 
and  by  emitting  phonons  rolls  down  through  these  excited  levels.  The 
capture  cross  section  resulting  from  this  model  is  proportional  to  T '4. 
V.  N.  Abakumov  et  al.  [37]  argue  that  the  capture  is  dominated  by 
electrons  of  energy  m*vj2  rather  than  KT,  where  vi  is  the  velocity  of 
sound,  and  m*  is  the  electron  effective  mass.  In  the  latter  case,  the 
capture  cross  section  is  calculated  to  be  proportional  to  T '3. 

The  cascade  capture  in  charged  centers  relies  on  closely  spaced, 
Rydberg-like  bound  states.  Neutral  centers  however  do  not  have 
excited  states  [38],  and  therefore  the  energy  must  be  released  by 
other  mechanisms.  During  capture,  the  released  energy  can  be  (a) 
emitted  in  the  form  of  light  (radiative  transitions),  (b)  transferred  to 
another  electron  or  hole  (Auger  process),  or  (c)  transferred  to  phonon 
modes.  The  phonon  emission  processes  include  cascade  mechanisms 
and  multiphonon  emission  (nonradiative  transitions). 

Nonradiative  transitions  play  a role  in  many  processes,  examples 
include  the  broadening  of  sharp  optical  transitions,  and  various  device 
degradation  mechanisms.  Defects  are  likely  to  act  as  lifetime  killer 


104 


105 


centers  [39]  if  they  have  favorable  nonradiative  (i.e.,  vibrational) 
properties  such  as  large-amplitude  modes  promoting  transitions  and 
large-energy  modes  to  take  up  the  electronic  energy  to  facilitate 
nonradiative  carrier  capture. 

Emission  processes  can  be  considered  as  the  reverse  processes 
of  carrier  capture  in  equilibrium.  As  far  as  field-enhanced  emission 
processes  are  concerned,  there  are  three  mechanisms  used  to  account 
for  its  enhancement:  (1)  the  Poole-Frenkel  effect,  where  the  electron 
climbs  over  a barrier  lowered  by  the  presence  of  the  field,  (2)  pure 
tunneling,  and  (3)  phonon-assisted  tunneling,  where  the  electron 
absorbs  thermal  energy  from  the  lattice  and  then  tunnels  through  a 
field-lowered  barrier  at  a higher  energy  [40]. 

In  this  chapter,  we  discuss  nonradiative  transitions,  specifically, 
as  they  apply  to  carrier  transitions  via  neutral  centers. 

6.2.  Microscopic  Capture  Model 

6.2.1.  Introduction 

The  Monte  Carlo  technique  has  been  used  to  unravel  the 
nonlinear  transport  phenomena  of  carriers  in  semiconductors  for 
more  than  a decade.  The  advantages  of  using  the  Monte  Carlo  (MC) 
technique  in  device  simulations  are  that  it  provides  an  accurate 
microscopic  description  of  the  physical  processes,  that  it  does  not 
require  any  assumption  on  form  of  the  carrier  distribution  function, 
and  that  it  is  very  easy  to  include  band-structure  effects.  In  order  to 
implement  this  technique,  microscopic  models  for  the  physical 
processes  should  be  developed.  It  is  difficult  to  use  the  MC  technique 
to  investigate  phenomena  occurring  in  a semiconductor  device  with 
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deep  generation-recombination  (G-R)  centers  because  there  are  not 
many  well-accepted  physical  models  for  G-R  process  available  to  be 
implemented.  Although  a cascade  capture  G-R  mechanism  for  a 
shallow  impurity  has  been  developed  [37]  and  used  to  calculate  G-R 
noise  [15],  and  nonradiative  multiphonon  (NMP)  carrier  capture 
models  by  attractive  and  repulsive  centers  have  been  described  by 
R.  Passler  [41,42],  microscopic  theories  for  neutral  deep  levels  can 
not  be  found. 

Defects  create  localized  electronic  states  in  the  forbidden  band 
of  semiconductors.  Phenomenologically,  the  interaction  between  deep 
states  and  conduction  or  valence  bands  can  be  described  by 
Shockly-Hall-Read  (SHR)  statistics  [43].  In  SHR  statistics,  the  rate 
equations  involve  capture  and  emission  rates  of  carriers.  Carrier 
capture  efficiencies  of  traps  are  quantitatively  described  by  capture 
coefficients,  which  give  the  actual  size  of  the  crystal  volume  swept  out, 
in  unit  time,  due  to  capture  into  the  individual  trap  considered.  A 
theory  for  a realistic  capture  coefficient  C(T,E)  should  describe  both 
its  temperature  and  carrier  energy  dependence.  The 
phenomenological  capture  coefficient  C(T)  results  from  the 
microscopic  C(T,E)  averaged  over  the  carrier  distribution.  Preliminary 
knowledge  of  the  dependences  of  microscopic  carrier  capture 
coefficients  C(T,E)  on  E thus  represents  in  general  the  prerequisite 
for  actual  calculations  of  any  observed  capture  coefficient. 

The  theory  of  nonradiative  multiphonon  capture  of  free  carriers 
belongs  to  the  very  few  fields  of  solid  state  physics  where  a satisfactory 
solution  has  not  been  achieved  so  far.  The  multiphonon  problem 
involves  a complex  many-body  Hamiltonian  and,  in  order  to  induce 
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transitions,  some  perturbation  [44].  The  G-R  process  can  be  described 
only  by  (1)  approximating  the  unperturbed  Hamiltonian,  (2)  assuming 
some  sensible  wavefunctions  for  the  unperturbed  Hamiltonian,  and  (3) 
assuming  some  perturbation.  Once  the  basic  approach  is  selected, 
there  are  still  many  possible  variations  regarding  the  electron-phonon 
interaction,  the  vibration  modes,  and  the  electron  wavefunctions. 

The  microscopic  nonradiative  multiphonon  (NMP)  carrier 
capture  by  attractive  or  repulsive  centers  has  been  described  by 
R.  Passler  [41,42],  The  complete  dependence  of  microscopic  carrier 
capture  coefficients  C(T,E)  on  variables  T and  E is  derived  from  the 
semi-phenomenological  relation  C(T,E)  = So(E)  C°(T,E),  where  So(E) 
is  the  Sommerfeld  factor,  accounting  for  the  influence  of  non- 
vanishing center  charges  Z I e I on  a free-carrier  wavefunction,  and 
C°(T,E)  takes  into  account  the  energy-loss  mechanism.  The  condition 
of  energy  conservation  can  be  satisfied  when 

p HCO  = Ek  + Eto, 

where  Ek  is  the  kinetic  energy  of  the  carrier,  Eto  is  the  energy  depth 
of  the  traps  below  the  band  edge,  p is  the  number  of  phonons  emitted 
during  the  capture  process,  and  CO  is  the  vibrational  radian  frequency. 
In  Passler' s papers,  p is  not  necessary  an  integer. 

Generally  speaking,  capture  of  a carrier  by  a deep  center  occurs 
in  three  steps:  firstly,  the  free  wavefunction  is  modified  by  the  electric 
field  surrounding  the  capture  center  when  the  carrier  approaches  the 
deep  level  center,  i.e.,  coulomb  scattering:  secondly,  the  electron 
releases  its  energy  through  delocalized  phonon  emission.  If  the 
process  goes  on,  the  capture  becomes  a cascade  process.  The  carrier 
arrives  at  a bound  vibronic  state  by  either  tunnelling,  by  a thermally 
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activated  transition,  or  by  a direct  transition.  Thirdly,  the  carrier 
drops  to  the  ground  state  by  multi-phonon  emission.  The  capture 
cross  section  is  proportional  to  the  product  of  the  three  transition 
probabilities  associated  with  these  three  steps.  If  any  step  is  rate 
limiting,  then  that  process  will  be  sensed  by  an  experiment  designed 
to  probe  the  capture  cross  section. 

In  contrast  to  the  semi-phenomenological  approach  by  R.  Passler 
[41,42],  we  derive  a microscopic  nonradiative  multiphonon  carrier 
capture  coefficient  for  a neutral  deep  center  based  on  quantum 
mechanics  in  the  following  sections.  In  our  model,  the  carrier  releases 
its  energy  through  delocalized  optical  phonon  emission  and  directly 
jumps  from  a free  vibronic  state  to  a bound  vibronic  state.  We  assume 
that  this  step  is  the  rate-limiting  one.  Therefore  the  capture 
coefficient  can  be  calculated  by  evaluating  this  transition  probability.  In 
our  model,  the  energy  conservation  condition  with  non-integer  p as 
used  by  R.  Passler  [41,42]  is  avoided  by  a mechanism  in  which  the 
delocalized  phonon  frequency  is  adjustable.  Both  the  long-range  and 
short-range  electron-phonon  interactions  are  accounted  for  by 
including  localized  and  delocalized  phonons  in  the  representation  of 
vibronic  states.  We  focus  on  silicon  material,  but  the  basic  ideas  on 
which  this  theory  is  based  apply  to  other  materials  as  well. 

6.2.2.  Theory 

In  this  section,  we  derive  the  transition  probability  for  a free 
carrier  to  be  captured  by  a deep  neutral  center.  The  direct  influence 
of  the  externally  applied  electric  field  on  the  discrete  levels  may  be 
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appreciable  only  if  the  field  strength  E is  the  order  of  the  interior  field 
in  the  lattice,  and  is  therefore  neglected  in  the  following  discussion. 

6.2.2. 1.  Transition  Probability 

The  Hamiltonian  of  the  impurity-electron  system  can  be 
generally  expressed  as 


Htot  = He  + HL+  HeL’ 

(6.1) 

with 

hL  = Tn(Q>  + Ki2q2  . 

(6.2) 

where  Tn(Q)  is  the  vibrational  kinetic  energy,  and  Q is  the  lattice 
coordinate;  He,  Hl,  HeL  are  the  electronic  part,  the  vibrational  part 
and  the  electron-phonon  interaction  of  Htot,  respectively.  The 
electron-phonon  interaction  Hamiltonian  HeL  for  a free  state  is 
smaller  than  He  and  Hl,  and  is  therefore  neglected  in  the  discussion 
of  a free  state. 

The  problem  is  to  determine  the  probability  Wet  of  finding  the 
system  in  the  impurity  ground  state  as  t ->  oo,  if  it  was  in  the  free 
state  at  t ->  - oo.  The  exact  eigenstates  T satisfy 

( Htot  " Et  ) ^ = °-  (6.3) 

This  equation  can  not  be  solved  because  the  Hamiltonian  does  not 
allow  a separation  of  variables.  Hence  a static  approximation  is  used. 

In  the  static  approximation  [45],  the  electronic  wavefunction  *Pk 
is  assumed  to  be  independent  of  the  lattice  coordinate  and  taken  as 

He  <Pk=  Ee  <Pk  • (6.4) 

The  total  wavefunction  is  written  as  [45] 


no 


T= 


I <Pt,r)  X^.Q' , 


(6.5) 


where  Q is  the  configuration  coordinate,  and  ^kt  is  the  lattice 
wavefunction  which  is  derived  in  the  following. 

Taking  the  electron-phonon  interaction  to  be  of  the  form  HeL=UQ. 
where  U is  constant,  and  substituting  eq.  (6.5)  into  eq.  (6.3),  the 
following  expression  is  derived  for  the  lattice  wavefunctions 


(HL+Ee  + X <klU  Qlj>  ‘ Et>  Xkt(r’9)  - °.  (6.6) 

k 


i.e., 

( HL  + Eg  + u 9 - Et)  X Jr.Q)  = 0. 

Kl 

This  can  be  written  as 


(6.7) 


U ,2 


(T, 


n(Q) 


+ Kl(9+K?  +f  u2 

2 e 2Ki  l kt 


- Et>  XJr’Q)  = o- 


Equations  (6.4)  and  (6.7)  give  a vibronic  state  wavefunction 

^k.t=  1Pk(r-9o'  Xkt(r.Q) . 

The  total  vibronic  energy  is  given  by 


(6.8) 


(6.9) 


with 


Et  = ["2  + nt)h(0  - Eto, 


2 2 

Eto=  - Ee  + 2K7  = Eeg  + ^ , 


where  Eeg  is  the  electronic  ground  state  energy  of  a carrier. 
We  define 

u2 

2KT  = S H(0’ 


(6.10) 


(6.11) 


(6.11a) 


where  S is  the  Huang- Rhys  factor,  and  S hCG  is  the  relaxation  energy. 
For  a carrier  in  a free  state,  HeL  is  negligible,  and  therefore 
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1 fi  k 

Ek  = (2  + nk)R(0  + ♦ , (6.12) 

where  Ek  is  the  energy  of  the  carrier  in  the  free  state,  k is  the  carrier 

wavevector,  and  m*  is  its  effective  mass. 

The  wavefunctions  lk,nk>  and  lt,nt>  can  be  written  in  the 

static  approximation  as 
lk,nk>=lk>  lnk>, 

I t,nt>=  1 1>  lnt>.  (6.13) 

The  wavefunctions  lk,nk>  and  lt,nt>  are  orthogonal,  because  of  the 
orthogonality  relations 


J 

Q 

I 


Q 


9j  9kdr  = 

\t2  \tl  dr  = 


5 

j.k  , 

5 

t2,  tl  , 


(6.14) 


where  Q is  the  crystal  volume. 

In  general,  the  transition  rate  between  two  states  can  be 

computed  using  the  Golden  Rule 

pi2  = — I <2IH'I  1>|2  5(E2-Ei),  (6.15) 

E 

where  I 1>  and  I2>  represent  many-electron,  many-phonon  states, 

and  Ei,  E2  are  the  corresponding  total  energies. 

Equation  (6.15)  can  be  written  as 

P12  = ~ I <t,nt,n2'  I H'  I k,nk,nj->  I 2 5(E2-Ei). 

E 

where  lk,nk,nj->  and  It,nt,n2’>  refer  to  the  electronic,  delocalized 
and  localized  wavefunctions  of  the  states  I 1>  and  1 2>  respectively. 
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We  assume  that  the  optical  deformation  potential  is  the  perturbation 
Hamiltonian  H'  which  is  given  by  [14] 

v E 1/2  i^'r  + 

H'=  D X (— ) e (aq+aq  ),  (6.16) 

q 2pQ(0q  H H 

where  p is  the  density  of  the  crystal,  D is  the  optical  deformation 
potential  constant,  aq  is  the  phonon  creation  operator,  aq+  is  the 
phonon  annihilation  operator. 

If  however  acoustic  deformation  dominates  carrier  capture,  the 
perturbation  Hamiltonian  is  given  by  [14] 

H’=DaZ  q (~  — )1/2  e (aq+aq),  (6.16a) 

q A 2pQ(0q  M H 

where  p is  the  density  of  the  crystal,  Da  is  the  acoustic  deformation 
potential  constant,  aq  is  the  phonon  creation  operator,  aq+  is  the 
phonon  annihilation  operator. 

The  transition  rate  for  capture  of  a carrier  by  vibronic  states  through 
emitting  an  optical  phonon  is  given  by 

W kt=  “ avenk  I (nq+l)l<t,nt  1 V'  I k,nk>  I 2 5(Ek~E^— ^-), 

E n t E 

(6.17) 

where  EtOq  is  the  energy  of  a delocalized  optical  phonon,  and  avenk  is 
an  average  over  the  quantum  states  of  delocalized  phonons.  The 
vibrational  wavefunctions  I nt>,  and  I nk>  are  assumed  to  have  the 
same  vibration  frequency  CO,  the  unit  of  8 is  second. 


V=D  X 

q 


(; 


2pf2tOq 


1/2 


iq  • r 

e 


(6.18) 
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and 


1 


exp(- 


KCOc 

KT 


4-1 


(6.18a) 


Assuming  that  the  dynamics  of  delocalized  and  localized  phonons  are 
independent,  we  get 


<t,nt  I V I k,nk>=  <tl  V'  I k><  ntlnk>.  (6.19) 

The  probability  of  a transition  involving  p phonons  can  be 
determined  by  summing  over  the  transitions  between  an  initial  state 
with  quantum  number  n and  a final  state  with  quantum  number  n+p. 
Equation  (6.17)  can  be  written  as  a sum  of  the  transition  probabilities 
involving  p phonons 

Wkt  = JT  (Hq+ l)Kt  IVIk>|2XPp  §(Ek-Et-Kfflq^  (6  20) 
H P 27tfi 

with 

Pp=  (l-exp(-^^))  X exp(  -^^l<n+pln>  I2  . (6.21) 

n 

Let  the  position  of  the  energy  Et  be  so  smeared  out  [46]  that  we 
may  describe  it  only  by  the  probability  Y(Et)  dEt  of  finding  it  between 
Et  and  Et  + dEt,  and  with  g(ft)  = 27t  h Y(Et)  so  that 
J g(ft)  dft  = 1. 

With  the  energy  conservation  connection  and  eq.  (6.18),  the 
transition  probability  rate  becomes 

Wkt=^(n  + l)l<tlV'lk>|2XPpg(ft)  . 

H P 
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This  equation  can  be  further  approximated,  if  we  assume  that  Et  is 
peaked  around  Eto- 

The  transition  probability  rate  from  the  free  state  *Pk  to  the 
bound  state  *Pt  becomes 


Wkt^D2(^rvi) 


J 


lq  • r 


•[J  dr]  271  S(CDq)X  P p , 

Q P 


(6.23) 


where 


h2k2 

P = 2m*  + Et°  ‘ h wq>  (6.24) 

where  *Pt  is  the  bound  electronic  wavefunction,  *Pk  is  the  free  state 
wavefunction,  (Dq  is  the  frequency  of  a delocalized  optical  phonon 
mode.  We  define  S(C0q)  such  that  S(C0q)  dCOq  is  the  fraction  of  modes 
with  radian  frequency  between  (Oq  and  COq  + d(Oq. 

At  low  temperatures  (p  > n)  [47],  we  get 


P 


P = 


g*5  -S(2n+ 1) 

pTe 


(n+i)P 


2 

S n(n+ 1)  , 
(1+  p+T  K 


where 

1 

n = 

.hoy. 

CXp  ( " 1 


(6.25) 


(6.26) 


The  way  that  a carrier  is  trapped  in  a vibronic  ground  state  is 
illustrated  in  Fig.  6.1.  When  a free  carrier  loses  its  energy  htOq  by 


Electronic  +elastic  energy 
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Figure  6. 1 . The  configuration  coordinate  diagram- 
electronic  + elastic  energy  versus  Q 
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emitting  an  optical  phonon,  it  directly  transits  to  a vibronic  bound 
state,  and  subsequently  drops  to  the  ground  vibrational  state  by 
multiphonon  emission. 

6.2. 2. 2.  Density  of  States  of  Delocalized  Optical  Phonons 

The  dispersion  curves  of  phonons  in  silicon  in  the  [1,0,0] 
direction  can  be  described  by  the  following  branches  [48]:  the 
transverse  branch 

MWq2=4£+4Y(l-cos(7tqi)br2V2T V [£2+q2+(^-T|2)cos(7Cq-|)]  f (6.27) 
and  the  longitudinal  branch 

M(0q2=4£  + 8y  ( 1 -cos(7rqi ))  ^ 2y[2  ^\j  l+cos(7Cqi)  (6.28) 

where  the  - sign  refers  to  the  acoustic  phonon  branch  and  the  + sign 
indicates  the  optical  phonon  branch,  qi=  q/qmax  is  the  reduced 
phonon  wavevector,  qmax  = 7i/a,  a is  the  lattice  constant  of  silicon 
(5.43095  x 10’8cm),  and  M is  the  mass  of  host  lattice  (4.666  x 10’23g). 
We  fit  eqs.  (6.27)  and  (6.28)  to  the  experimental  results  [49],  and  get 
the  constants: 

£ = 0.5804  x 105  dynes/cm, 

Y=  0.4069  x 102  dynes/cm, 

T|=  0.4667  x 105  dynes/cm. 

The  fitted  results  are  shown  in  Fig.  6.2. 

The  density  of  states  for  optical  phonon  is  calculated  from  the 
following  formula 
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Figure  6.2.  Phonon  dispersion  curves  of 
silicon.  The  full  lines  are  the  fitted  results. 
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fCOu  „ 

! (a(fi)q)  )dCOq= 

Q rqxm 

rqym 

fqzm  dqxdqydqz. 

(6.29) 

^ (Oi 

87Z3J  „ 

qxm* 

■qym 

J -qZm 

where  (0U  and  (01  are 

the  upper 

and 

lower  limit  of  the 

radian 

frequency  for  each  dispersion  branch,  and  qjm  is  the  upper  limit  of  qj. 
For  simplicity,  we  approximate  the  density  of  states  for  phonons 


by 


H(mq)  = 2^(^)2-^-  . 
” 8jr  a dCOn~ 


(6.30) 


87t  “ uujqz 

From  eqs.  (6.27)  and  (6.29),  we  get  for  the  density  of  states  of 
the  transverse  optical  branch 


Et(C0q)= 


Q 

a 71 


~VM{4^+8Y(l-cos(7tqi))+2V2~V  [^2+q2+(^2-q2)cos(7tqT)]} 


(^-Tl2)sin(7tqi) 


(6.31) 


V2[^2+q2+(£“-rr)cos(7Cqi)] 


7T~ Z 


- 2Ysin(7iqi) 


The  density  of  states  of  the  longitudinal  optical  branch  is  given  by 


^ VM{4^+8Y(l-cos(7tqi))+2V2^Vl+cos(7rqi)} 

■=-l(Wq)=  -3 — — . (6.32) 

a K qsin(Ttqi)  „ . , 

",  = - 4y  sinftqi ) 

'V2[l+cos(7cqi)] 


6.2. 2. 3.  Evaluation  of  Electronic  Matrix  Elements 

A theory  for  the  magnitude  of  the  capture  cross  section  is 
difficult  to  develop  because  the  explicit  form  of  the  electron 
wavefunction  for  a deep  state  is  not  well  known.  The  initial  and  final 
electronic  wavefunctions  are  rather  complicated,  and  depend  on  what 
we  assume  for  the  potential  of  the  defect.  Since  the  defect  potential  is 
not  explicitly  known,  we  assume  the  potential  of  a neutral  deep  defect 
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to  be  a spherical  square  well  of  depth  Vo,  and  radius  b.  If  we  further 

assume  the  potential  is  a Lucovsky's  delta-function  [50],  then 
2 2 

Vo  b = 7t  H /(8m)  as  b ->  0.  The  electronic  wavefunctions  form  an 
orthogonal  set.  The  approximate,  normalized  ground  state  electronic 
wavefunction  is  given  by 

- a r 

<Pt(r)=V^  ^T~  • (6.33) 


with 


a= 


V2m  Fe~g 
li2 


(6.34) 


where  m is  the  electron  mass. 

For  simplicity,  the  free  electronic  wavefunction  is  approximated  by 


ik  • r 

e 


The  electronic  matrix  element  is  then  given  by 


(Pk(r)  = 


iq  r 

<Pte  <Pkdr  = 


!V 


27ta 

n 


l 


2 ,r 
\CC  +lk+ql 


(6.35) 


(6.36) 


Because  Ik+ql  =lk+q  cos(@)  I , where  0 is  the  angle  between  k and  q, 
the  average  electronic  matrix  element  is  given  by 


f 1 A 

2 — — 2 
\a  + lk+q  I J 


dcos(0)d(|) 


/ 2 2 2 \ 

(a  +k  +q  +2kq  ) 

2 2 2 

ya  +k  +q  -2kqy 


where  c|)  is  the  azimuthal  angle. 


(6.37) 


120 


6. 2. 2. 4.  Calculation  of  the  Capture  Coefficient 

With  eqs.  (6.23),  (6.25),  and  (6.36),  and  by  assuming  that  the 
transition  probability  mainly  comes  from  the  p-phonon  emission 
process,  the  transition  probability  via  optical  phonon  emission 
becomes 


Wkt(T,E,p)=  -^D2(._^Me2  . + 1}  2k  (Sl(C0q)+2  Et(COq)) 


j-2  2pQC0q 

SP  -S(2n+l)  _ p s^n(n+l) 

P~e 


(n+i)1  (1+  " 'p7i  " )•  (6.38) 

with  the  conditions  that  the  optical  phonon  wavevector  satisfies  eqs. 
(6.27)  and  (6.28),  and  the  carrier  energy  satisfies  eq.  (6.24). 

The  microscopic  capture  coefficient  is  given  by 

C(T,E,p)=  Q W ct(T,E,p).  (6.39) 


The  phenomenological  capture  coefficient  C(T)  results  from  the 
microscopic  C(T,E)  after  averaging  over  the  carrier  distribution.  The 
capture  cross  section  is  given  by 


G (T,E)= 


I 

P 


m*  Q Wct(T,E,p) 
hk 


The  average  capture  coefficient  reads 


(6.40) 


J n(E)C(T,E,p)  d E 

C(T)=  X 

P J n(E)  d E 


(6.41) 


where  n(E)  represents  the  carrier  distribution  function. 
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The  observed  capture  cross  section  is  C(T)  divided  by  the 
thermal  velocity  of  the  carriers.  In  the  following,  we  calculate  the  value 
of  the  electron  capture  coefficient  of  the  gold  acceptor  in  silicon. 
Equation  (6.41)  can  be  rewritten  as 

J D(E)f(E,T)C(T,E,p)  dE 

C(T)=  X . (6.42) 

P J D(E)f(E,T)dE 


where  D(E)  represents  the  density  of  states  in  the  conduction  band, 
and  f(E,T)  is  the  carrier  distribution  function. 

At  low  field  and  Ec-Ef  »KT,  f(E,T)  is  given  by 

h2  k2  V Ef 
2m*  KT  ‘ KT 

f(E,T)  = e e (6.43) 


With  the  following  expressions  for  D(E)  and  the  various  masses, 
V2MC  (mde)3/2  E1/2 


D(E)= 


h3  2 
n K 


V2Mc(mde)3/2E1/2k 
D(E)dE=  : — o dk, 


m*  h n 


1/3 


mde  =(mi*  mt*  mt*) 


nr 


=3  h* 


2 V1 

,*+  . 


mf  mt 


(6.44a) 

(6.44b) 


we  get  for  the  root  mean  square  velocity 

vrms  = ■ (6.45) 

Considering  that  the  transition  probability  is  dominated  by  the 
p-phonon  emission  process,  we  can  estimate  the  capture  coefficient 
from  eq.  (6.38), 
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_ , _ 1 , E , ^ 2„  , „ „ _ 2„ 

C(T,E,p)  = fi  — ' 
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} 


sp  _s(2n+l) 
* pfe 


2 

p s n(n+ 1)  % 
(n+l)  (1+  p+1  ’ 


(6.46) 


where  Di  is  the  longitudinal  optical  deformation  constant,  Dt  is  the 


transverse  optical  deformation  constant. 

If  we  take  the  radian  frequency  of  localized  optical  phonon  as 
(0=0.1218  x 10+15/s,  and  use  the  following  electronic  constants 


mi*=0.98  * 0.91095  x 10'27g,  mt*=0.19  * 0.91095  x 10'27g, 

a=5. 43095  x 10‘8cm,  Mc=l,  M=4.666  x 10'23g, 

T 3/2 

Nc=  2.8x  1019  (3^0)  /cm3,  p=2.328  g/cm3, 

Di=llx  108eV/cm,  Dt=2  x 108eV/cm 

for  the  gold  acceptor  level  in  silicon,  we  can  calculate  the  electron 
capture  cross  section  with  eqs.  (6.45)  and  (6.46  ) if  we  use  Eto=0.540 
eV,  S=2.743  and  p=6.  The  results  are  shown  in  Fig.  6.3.  The  theory 
gives  values  for  the  electron  capture  coefficient  which  are  of  the  same 
order  as  the  experimental  data.  The  slight  temperature  dependence  of 
our  model  may  be  due  to  some  of  the  approximations  we  used  and 
might  not  be  real. 
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Figure  6.3.  Electron  capture  cross  section  of  the  gold 
acceptor  level  versus  1000/T.  The  full  line  represents 
the  results  of  the  microscopic  model  outlined  in  the 
text. 


CHAPTER  VII 

DEFECT  PARAMETER  EXTRACTION  FROM  FIELD-EFFECT 

TRANSISTORS 


Junction-gate  transistors  are  particularly  attractive  for  use  in  the 
input  stages  of  low-frequency,  low-noise  amplifiers,  and  offer 
substantial  improvements  over  bipolar  amplifiers  when  the  signal 
source  impedance  is  large.  MOS  transistors,  on  the  other  hand,  tend 
to  generate  considerably  large  amounts  of  low-frequency  noise. 
However,  in  the  MHz  range  their  noise  is  comparable  to  the  junction- 
gate  device,  and  this,  together  with  a very  low  cross-modulation 
performance,  make  them  particularly  useful  for  low-noise  RF 
amplification. 

There  are  three  major  noise  sources  in  MOSFETs  and  JFETs: 
the  thermal  noise  in  the  channel,  the  1/f  noise,  and  the  generation- 
recombination  (G-R)  noise  due  to  random  emission  of  electrons  and 
holes  from  defect  centers  in  the  depletion  region  of  the 
semiconductor  and  the  G-R  noise  originating  from  the  random 
trapping  and  releasing  of  current  carriers  from  deep  centers  or 
shallow  impurities  in  the  channel. 

In  this  chapter,  we  analyze  the  posibilities  for  using  the  G-R 
noise  spectra  of  MOSFETs  and  JFETs  as  a tool  for  trap  parameter 
extraction.  We  discuss  MOSFETs  first,  and  then  describe  JFETs. 
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7.1.  Fluctuations  in  the  Charge  States  of  Shocklev-Hall-Read  Centers 
Located  in  the  Space  Charge  Region 


The  rate  equation  for  the  trapped  electron  concentration,  taking 
the  four  processes  indicated  in  Fig.  7. 1 into  account,  may  be  written  as 
dnt 

“dF  = Nt  [(cn  ftp  N - en  ft)  - (cp  ft  P - ep  ftp)],  (7.1) 

where  P,  N are  the  hole,  electron  concentration  in  the  valence  band 
and  conduction  band,  respectively;  nt,  Nt  are  the  occupied  and  total 
trap  concentration,  respectively;  cp,  ep  are  the  hole  capture  and 
emission  coefficient,  respectively;  cn,  en  are  the  electron  capture  and 
emission  coefficient,  respectively;  ft  is  the  fraction  of  occupied  traps. 
We  express  in  terms  of  Fermi-Dirac  statistics 


E 


1 +gte 


f ' n 

t rr~ 


(7.2) 


where  Et  is  the  trap  energy  level,  Ef  is  the  Fermi  level  and  gt  is  the 
degeneracy  factor. 

Using  the  principle  of  detailed  balance  and  eq.  (7.1),  eq.  (7.2) 
can  be  rewritten  as 


f.  =1-  f,  = 
tp  t 


cp  P +en 


cp  P +en  +cn  N +ep 
Using  the  principle  of  detailed  balance  and  eq.  (7.2),  we  get 


(7.3) 


eP  =cp  Nv£te 


E - E, 
v t 

"FT 


(7.4) 


Similarly, 
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Figure  7. 1 . The  electron  and  hole  generation-recombination 
processes:  (a)  electron  capture,  (b)  electron  emission, 

(c)  hole  capture,  (d)  hole  emission. 
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The  carrier  densities  can  be  expanded  around  their  steady-state  values 


2 

II 

o* 

+ AN, 

(7.6) 

P = Po 

+ AP, 

(7.7) 

nt  = nto 

+ Ant> 

(7.8) 

where  No,  Po,  nto  are  the  electron,  hole  and  occupied  trap 
concentration  in  steady-state,  respectively;  AN,  AP,  Ant  are  the 
corresponding  density  fluctuations. 

Using  eqs.  (7.1)  and  (7.6)  through  (7.8),  we  get  after  linearization 
dAnt 

dt  ='[(cn  N0+en)+(cp  P0+ep)]Ant  +cn(Nt-nto)AN-  cp  nt0  AP.  (7.9) 

The  noise  from  the  carrier  fluctuations  AN  and  AP  may  be  neglected, 
since  they  have  much  shorter  time  constants,  i.e.,  the  transit  time 
across  the  transition  layer,  than  Ant.  For  a nearly  depleted  region  in  a 
semiconductor,  the  solution  of  eq.  (7.9)  is  relatively  simple,  if  we 
concentrate  on  low-frequency  fluctuations  only.  Thus  we  get 

l(cn  N0  + en  )+(cp  P0  +ep  )].  (7.10a) 

When  the  carrier  densities  are  very  small,  it  becomes 

T~=  [en  +ep  ].  (7.10b) 

lp 

If  en  » ep  or  en  « ep  holds,  as  for  example  in  the  case  of  the 
acceptor  level  of  gold  in  silicon  where  en  is  20-50  times  greater  than 
ep  [8],  we  are  able  to  extract  the  activation  energy  of  the  trap  from  a 
measurement  of  the  lifetime  of  the  G-R  noise  spectrum  versus 
temperature.  Equations  (7.4)  and  (7.5)  are  often  assumed  to  be 
independent  of  electric  field,  and  the  activation  energies  and  capture 
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coefficients  are  therefore  extracted  from  these  equations.  In  reality 
however,  field-dependent  factors  should  be  taken  into  account  in  the 
interpretation  of  experimental  data. 


7.2.1.  Low-Field  Noise  Theory 

The  MOS  transistor  device  structure,  with  an  induced  n-channel 
structure  in  a p-type  substrate,  is  shown  in  Fig.  7.2.  According  to  the 
theory  of  low-frequency  G-R  noise  derived  by  Yau  and  Sah  [52],  the 
equivalent  gate  voltage  spectral  density  is  of  the  form 


7.2.  Generation-Recombination  Noise  in  Metal-Oxide- 
Semiconductor  Field-Effect  Transistors 


ft(l-ft)  1-^-  rdx  dV, 


-,2 


(7.11) 


with 


1 


tp  en  +ep 


(7.12) 


2 • 


(7.13) 


and 


i_  » 
2 


<>F  +V(y))  ] 

3. 
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Figure  7.2.  MOS  transistor  structure 
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Co=  ~r  . 


(7.14) 


o 


where  V=V(x)  is  the  voltage  along  the  channel,  W is  the  depletion 
region  thickness  of  the  MOSFET,  is  defined  by  = (Ef  - Ei)/q,  Ej 
is  the  intrinsic  Fermi  level  in  the  bulk;  x0  is  the  thickness  of  the 
silicon  dioxide  of  the  gate,  and  eQ  is  its  associated  dielectric  constant; 
Vgsa  is  the  gate-source  voltage  minus  the  threshold  voltage,  Vds  is  the 
drain-source  voltage;  x is  the  distance  from  the  silicon-dioxide 
interface,  and  Na  is  the  shallow  acceptor  concentration. 

Assuming  that  tp  and  ft  are  constants,  eq.  (7.11)  becomes 


Svgs(t0)-  Fs- 


(7.15) 


where  Fg  contains  the  geometrical  factors  and  the  d.c.  bias  voltages 
which  are  kept  constant  during  a particular  noise  measurement  versus 
temperature. 

Using  eq.  (7.13),  and  if  ep  « en,  eq.  (7.15)  may  be  simplified  to 


Sygs^  °>  . ££ 
Fg  tp  en 


(7.16) 


Combining  eqs.  (7.12)  and  (7.16),  we  get  ep  and  en.  Constructing 
Arrhenius  plots  and  using  eqs.  (7.3)  and  (7.4),  we  get  the  activation 
energies  of  the  trap  involved.  Recently,  Scholz  et  al.  [9]  used  this 
technique  and  DLTS  to  measure  trap  parameters  and  found  that  the 
results  from  both  techniques  are  in  agreement. 


7.2.2.  Field  Dependent  Noise  Theory 

According  to  the  arguments  of  Yau  and  Sah  [8],  the  major  noise 
contribution  comes  from  the  random  thermal  emission  from  gold 
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centers  located  near  the  Interface.  The  electric  field  can  be  varied  by 
varying  the  substrate  bias.  If  we  approximate  the  electric  field  by 
E = Emax  = 2Vbs/W(V=0),  where  Vbs  is  the  substrate-source  voltage, 
the  field  dependent  emission  coefficients  can  be  extracted. 


7.3.  Generation-Recombination  Noise  in  Junction-Gate  Field- 

Effect  Transistors 

7.3.1.  Low-Field  Junction  Generation-Recombination  Noise 

The  JFET  device  structure,  with  a n-channel  structure,  is  shown 
in  Fig.  7.3.  The  low-frequency  excess  noise  in  junction-gate  FET's  was 
first  analyzed  in  detail  by  Lauritzen  and  Sah  [53,54,55].  The  theoretical 
expression  for  G-R  noise  due  to  SHR  centers  in  the  gate-channel 
transition  region  of  an  abrupt  junction,  symmetrical  JFET  is  written  in 
the  form  [54,55] 


3 2. 
a 9 tp 


6 e 2 L Z(1-K02  tp2) 


ftd-ft)  F(us,ul), 


(7.17) 


where  £g  is  the  dielectric  constant  of  the  substrate. 
The  function  F(us,Ul)  is  given  by 


F(us,Ul)= 


|f3(us2-  ul2)-2(us3-  uL3)][ln(~)-2(us-  uL)+  ^tus2-  uL2)] 


(us-  uj5 


with 


and 


[Vgg+Vbgj1/2  =1 

V. 


U 


pd 


s. 


uI=l-[Vfis^ds+Vbfi]1/2  .1-Uri, 


pd 


ud. 


(7.18) 

(7.19) 

(7.20) 
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Fig.  7.3.  Schetch  of  an  n-channel  junction-gate  FET  showing 
coordinate  system  and  dimensions  used  in  the  analysis. 
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where  Vpd  is  the  pinch-off  potential,  and  Vbg  is  the  built-in  potential 
of  the  gate  junction. 

Above  pinch-off,  eq.  (7.20)  is  replaced  by 


with 


r aid'  1/2 


Id  /j,-  2 ri/  3 3* 

j — = 3(us  - uL  )-2(us  - u ) 

do  ^ 


(7.22) 


(7.23) 


Expressions  (7.12)  and  (7.13)  for  and  tp  respectively,  indicate 

that  defects  with  energies  away  from  the  center  of  the  gap  give 
negligible  noise.  K.  K.  Wang  et  al.  [56]  reevaluated  and  tp  by  not 

neglecting  the  free  carrier  density  density  at  the  edge  of  the  depletion 
region.  The  drain  current  noise  is  then  given  by 

„ „ [(Ud1/2-  Us1/2)-ln(“7 

Nd  a L"  (1+0)2  tp2) 

where  J0  depends  on  and  the  geometry  of  the  depletion  region.  Id 
is  the  d.c.  drain  current,  Nd  is  the  shallow  donor  concentration  in  the 


Sdi(<0)=  Wd  J°  [(Ud1/2-  Us1/2)-ln(E^T72i].  (7.24) 

Q1  Nd  a Lz  l+(02tn2  l- Us 


channel,  |in  is  the  electron  mobility,  and 
2 

tp~  3e„  ' 


(7.25) 


Hence,  a noise  measurement  of  tp  (T)  versus  T gives  the  defect 
activation  energy. 


7.3.2.  Generation-Recombination  Noise  in  the  Channel  of  Junction 
Field-Effect  Transistors 

The  expression  for  the  drain  current  G-R  noise  from  the 
channel  region  has  been  derived  by  different  authors  [4,54,57].  It  can 
be  written  as 
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Sdi(c°)  = 


4 Id 


V 


ds 


Ha  k dv 

(i+co2tp2)  * 


(7.26) 


with 


a = Nt  ft(l-ft)  / Nd. 

At  relatively  low  fields,  oc,  ft  and  tp  are  constants,  and  eq.  (7.26)  can  be 
written  as 

4 la  Vac  n a u„  t„ 

(7.27) 


4 Id  Vds  q OC  |in  tp 

Sdi(co^  = t2  , 

a L (1+CO2  tp2) 


Using  similar  techniques  as  presented  in  Chapters  II  through  V, 
trap  parameters  can  be  extracted  from  channel  G-R  noise. 

In  summary,  using  eqs.  (7.15)  and  (7.17),  we  can  extract,  from  a 
measurement  of  the  roll-off  time  constant  as  a function  of 
temperature,  the  trap  activation  energies  from  defect  centers  located 
in  junction  depletion  regions.  The  trap  parameter  extraction  from 
gate-junction  G-R  noise  is  complicated  by  the  geometrical  factors  of  a 
device  and  the  assumption  of  field-independent  capture  and  emission 
coefficients  required  to  take  the  frequency-dependent  factors 
tp / ( 1 -4- CO 2 tp2)  outside  the  integral.  G-R  noise  from  the  junction 
transition  region  of  a MOSFET  or  FET  is  generation  noise,  while  that 
from  the  channel  or  from  a bulk  p-type  or  n-type  semiconductor  is 
trapping  and  detrapping  noise. 


CHAPTER  VIII 

GENERATION-RECOMBINATION  NOISE  IN  POLYSILICON 

RESISTORS 

8.1  Introduction 

The  use  of  polycrystalline  silicon  has  become  increasingly 
widespread  in  solar  cells  and  VLSI  circuits,  including  polycrystalline 
MOS  devices,  static  RAM's,  PROM's,  ECL,  D/A  converters,  CCD 
detectors,  distributed  RC  filters,  and  as  redundancy  fusing  elements. 
The  polysilicon  film  is  composed  of  aggregates  of  single-crystalline 
silicon  grains  with  grain  boundaries  between  adjacent  grains.  The 
grain  boundaries  contain  trapping  states  that  are  capable  of  trapping 
mobile  carriers  and  contribute  to  the  creation  of  space  charge 
potentials,  and  may  also  serve  as  a sink  for  the  dopants,  which  cause 
them  to  become  inactive  [58].  Several  distribution  functions,  such  as 
an  ideal  5-distribution,  uniform  distribution,  Gaussian  distribution  [58], 
U-shaped  distribution  [59]  and  exponential  distribution  [60],  have 
been  used  to  describe  these  grain-boundary  trapping  states.  Since  G-R 
noise  is  generated  by  carrier  transitions  between  the  conduction  band 
and  trapping  states,  it  possibly  can  be  used  to  investigate  and  shine 
light  on  trapping  related  phenomena. 

Several  researchers  [61,62]  have  studied  the  noise  properties  of 
low  or  moderately  doped  n-type  polysilicon  films.  A large  amount  of 
pure  1/f  noise  was  observed  at  low  frequencies.  Number  fluctuations 
and  mobility  fluctuations  can  be  the  origins  of  1/f  noise 
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[17].  The  observed  1/f  noise  in  general  is  attributed  to  electron 
mobility  fluctuations  occurring  in  the  depletion  region  near  the  grain 
boundaries.  There  is  not  a general  theory  describing  number 
fluctuation  related  1/f  noise  in  polycrystalline  resistors.  Recently, 
M.  Y.  Luo  [61]  found  noise  spectra  of  the  form  C/(1+(G2  tp2)  and 
D/(l+(04  tp4)  on  LPCVD  phosphorus-doped  thin  films  with  a doping 
concentration  of  5.4  xlO  cm'  , where  D is  slightly  frequency 
dependent. 

In  this  chapter,  we  analyze  G-R  noise  models  associated  with  the 
defects  in  the  grain  boundary  regions.  We  explore  lightly  and 
moderately  doped  n-type  materials  for  the  sake  of  simplicity  and  hope 
that  it  will  stimulate  further  research  on  the  uncommon  D/(1+C04  tp4)- 
type  noise. 


8.2.  Theory 

The  polycrystalline  semiconductor  material  we  investigated  is 
assumed  to  satisfy  the  following  conditions:  (1)  the  grain-boundary 
space  charge  region  is  fully  depleted  and  its  width  is  small  compared 
to  the  grain  size,  (2)  its  physical  properties  are  the  same  in  adjacent 
grains,  (3)  the  boundaries  may  be  considered  as  surfaces  characterized 
by  an  effective  interface  state  density  nss. 

The  equilibrium  energy  band  diagram  of  the  investigated 
polysilicon  resistors  is  illustrated  in  Fig.  8.1.  The  symbol  Vb 
represents  the  potential  barrier,  w is  the  width  of  the  depletion 
region,  and  2d  is  the  width  of  the  grain-boundary.  We  use  the  simple, 
symmetrical  Schottky-barrier  band  model  [63],  and  assume  that 
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Figure  8.1.  Energy  band  diagram  of  a single  grain-boundary  region 
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the  voltage  drop  across  a grain  boundary  region  is  equally  distributed 
over  each  side  of  the  depletion  region  of  a grain  boundary.  To  model 
charge  transport  in  Schottky  diodes,  several  authors  [64]  have 
combined  the  diffusion  theory  and  the  thermionic  emission  theory  by 
considering  the  two  mechanisms  in  series.  The  current  flow  is 
described  by  the  diffusion  theory  in  the  depletion  region,  while  at  the 
metal-semiconductor  interface  the  net  current  is  described  by 


I = A q Vr  (N(d)-Nth(d)),  (8.1) 

where  N(d)  is  the  electron  density  at  x=d  when  the  current  is  flowing, 
and  Nth(d)  is  the  electron  density  at  x=d  when  the  current  is  equal  to 
0.  The  parameter  Vr  is  the  effective  recombination  velocity  at  the 
interface.  The  general  current  voltage  relationship  for  one  grain  can 
then  be  written  as  [64] 


VrVd 


ig  = 2 q Ag  e 


-qVb 

KT 


q vg 

sinhq^j^p. 


(8.2a) 


with 


VS=N^-  <8-2b> 

where  Va  is  the  applied  voltage,  Vg  is  the  voltage  drop  across  a grain, 
Ng  is  the  effective  doping  concentration  in  the  grain-bulk  region,  Nfg 
is  the  effective  number  of  grains,  and  Ag  is  the  cross  section  of  a grain. 
The  total  current  is  then  given  by 


I = 2 q A Ng(; 


VrVd 


Vr+Vd 
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-qVb 

KT 


sinhf 


q..Yg_ 

2 Nfg  KT 


(8.3) 


where  A is  the  cross  section  of  the  device. 
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In  polysilicon  resistors,  G-R  noise  sources  can  originate  from 
possibly  three  regions:  the  grain-bulk,  the  depletion  and  the  grain- 
boundary region.  The  G-R  noise  in  the  grain-bulk  regions  can  be 
calculated  using  the  bulk  G-R  noise  theories  given  in  Chapters  II  to  V 
and  some  small  signal  circuit  calculations  if  we  assume  that  the 
dynamic  resistance  of  the  grain-boundary  and  depletion  region  is 
larger  than  that  of  the  grain-bulk  region.  With  eq.  (2),  we  can  calculate 
the  dynamic  resistance  associated  with  a single  grain-boundary  region 
from  the  following  equation 


Rg 


(8.4) 


The  general  theory  for  G-R  noise  generated  by  traps  in  the  depletion 
region  has  been  described  elsewhere  [65].  The  equivalent  saturated 
diode  current  of  the  noise  in  the  depletion  region  can  be  written  as 
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dx, 


(8.5a) 


where  tp  is  the  trap  time  constant,  and  IeqL  and  IeqH  are  complicated 
functions  of  material  and  device  parameters. 

The  equivalent  output,  saturated  diode  current  of  the  noise  in  the 
depletion  regions  then  becomes  after  taking  circuit  considerations 
into  account 


A 2 
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(8.5b) 


Next,  we  focus  on  G-R  noise  from  the  trap  states  in  the  grain- 
boundary region.  Assuming  that  the  depletion  approximation  applies, 
Poisson's  equation  becomes 
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d2  V 
d2x 


q Ng 


d < x<  w+d  . 


(8.6) 


The  potential  barrier  height  (Vt>-Vgi)  is  the  difference  between  V(w+d) 
and  V(d),  that  is 


q w Nrt 


(8.7) 


where  Vgi  is  the  voltage  drop  in  the  depletion  region,  and  6S  is  the 
dielectric  constant  of  the  grain-bulk  material. 

Hence, 

(Vb-Vgi) 


w = 


Ng 


(8.8) 


For  simplicity,  we  limit  our  discussion  to  the  condition  Vb>Vgi.  From 
the  charge  conservation  relation,  we  get 


2 q w Ng  = q nss,  (8.9) 

where  nss  is  the  effective  grain-boundary,  charged  trap  density  per 
unit  area. 

The  current  equation  can  be  described  by 

I=AJ=  Aq[Dnf-p„N^i,  (8.10) 

where  |ln  is  the  electron  mobility,  and  Dn  is  the  electron  diffusion 
constant. 

The  electron  density  is  given  by 
, -q(ttn-yn  ) 

N=  Nc  e KT  (8.11) 

where  Tn  is  the  electron  potential,  and  On  is  the  electron  quasi-Fermi 
level. 
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In  the  following  discussion,  we  neglect  mobility  fluctuations. 
Expanding  the  carrier  densities  around  steady-state  in  first  order,  the 


current  equation  can  be  written  as 

dAN  dA'l'n  d'Fno 

AI=  A q[Dno  *T  A.  No  ~ +no  AN 

Using  eq.  (8. 12)  and  the  following  equations 
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Using  eq.  (8.14)  and  the  following  equations 
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we  get 
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The  derivative  of  the  potential  fluctuation  is  calculated  as  follows: 
dTn  g (w+d-x)  Ng 

to  = £s 

dA'Pp  q Aw  Ng 

Substituting  eq.  (8.19)  into  eq.  (8.18),  we  get 
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which  can  be  written  as 
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Using  eq.  (8.1),  the  current  fluctuation  can  be  written  as 
AI  = A qVr  AN(d). 

Thus  the  carrier  fluctuations  at  x=d  becomes 
ANId)= 

Substituting  eq.  (8.21b)  into  eq.  (8.20),  we  get 
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(8.22a) 
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(8.22b) 


Considering  only  the  low-frequency  noise  behavior  of  the 
device,  we  can  assume 

AN(w+d)  =0,  (8.23) 

then  eq.  (8.22b)  becomes 

q Aw  WefoNg 
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Fluctuations  in  the  charged  trap  density  of  the  grain-boundary  region 
cause  fluctuations  of  the  depletion  width  w, 

Anss  = 2 Aw  Ng.  (8.25) 

Thus 
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The  noise  spectral  density  becomes 

A q2  UN  W , 2 
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(8.26) 


Trap  carrier  fluctuations  are  calculated  as  follows.  Let  nss=2d  nt, 
where  nt  is  the  trap  density  of  a grain  boundary.  The  rate  equation  is 
given  by  [Chapter  VII] 
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The  spectral  density  of  G-R  noise  has  the  form  [17] 


S ((q)  = — — 

nt  2d  Ag  (1+CO2  tp2)  ’ 

where  g (=  en  nt ) is  the  electron  generation  rate. 

The  output  current  spectral  density  of  G-R  noise  becomes 

c fm,_  r cZefo,2  8 d g tp2 
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(8.29) 


(8.30) 


With  eq.  (8.20b),  and  assuming  that  the  Fermi  level  does  not  vary 
much  throughout  the  depletion  region,  we  get 


Wefo= 


w.+d  r -o0"! 

dx  = w 

-q°n(w+d) 

_ K T 

J Le  J 

Le  J 

= w e 


-q(Vb  -Vgl  +Vn) 

KT 


(8.31) 


where  Vn  = Ln(Ng/Nc). 

Equation  (8.30)  describes  the  G-R  noise  from  one  depletion 
region.  Assuming  that  the  fluctuations  of  barrier  heights  are 
independent.  After  some  circuit  manipulation  we  get  for  the  total 
device  1/f  noise 
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In  the  preceding  discussion,  we  have  assumed  that  the  trap 
states  have  a 8- distribution  function.  If  however  the  energy  of  the  trap 
states  is  distributed  over  a specific  energy  range  such  that  [17] 
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(8.33) 


where  (3,  Ti,  T0  are  functions  of  material  properties,  the  output  current 
spectral  density  of  G-R  noise  becomes 


Equation  (8.34a)  describes  the  1/f  noise  from  one  depletion 
region  of  a grain,  we  get  after  some  circuit  manipulation 


Equation  (8.34b)  describes  1/f  noise  produced  by  G-R  processes  in  the 
grain-boundary  region.  Since  the  observed  1/f  noise  is  proportional  to 
the  square  of  the  device  current,  the  fluctuations  of  the  charged  trap 
density  in  the  grain-boundary  can  not  be  used  to  account  for  the 
published  1/f  noise  in  polysilicon  resistors. 
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CHAPTER  IX 

CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

9.1.  Conclusions 

We  have  studied  deep  trap  parameter  extraction  methods  in 
semiconductors  using  G-R  noise  measurements.  In  the  course  of  our 
experimental  investigations,  we  found  that  proper  device  design  is 
critical  for  application  of  these  methods.  Deep  trap  parameters  can  be 
measured  on  devices  with  an  intended  deep  defect  concentration 
larger  or  smaller  than  the  shallow  dopant  concentration. 

Gold  in  silicon  has  been  widely  studied  in  the  past,  its  behavior 
is  not  always  predictable  however.  In  this  thesis  noise  measurements 
are  described  with  which  we  look  at  the  role  of  this  impurity  in 
silicon.  We  carefully  set  up  three  experiments  to  fully  exploit  noise 
measurement  techniques.  The  first  experiment  deals  with  low-field 
noise  measurements,  which  were  done  on  p+-p‘-p+  silicon  devices 

I c o 

doped  with  a boron  concentration  of  10  cm'  and  an  initially 
unknown  gold  concentration.  The  measured  trap  energy,  hole  capture 
cross  section  and  degeneracy  factor,  which  were  extracted  from  the 
experimental  noise  and  resistance  data,  are  in  good  agreement  with 
the  generally-accepted  data  for  the  gold  donor  level,  indicating  that 
this  method  is  accurate  for  investigating  deep-defect  properties. 

Since  defects  can  exist  in  several  charged  states  in  a device, 
other  states  may  affect  the  measured  data  collected  for  a particular 

energy  level.  Our  second  experiment  was  done  on  p+-p"-p+  silicon 
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devices  doped  with  a boron  concentration  of  2.7  xlO  cm'  . The  plot 
of  the  device  resistance  versus  temperature  indicates  that  gold 
produces  an  activation  energy  of  0.267  eV,  and  the  plot  of  the  voltage 
noise  versus  temperature  indicates  that  gold  produces  an  activation 
energy  of  0.336  eV.  These  results  are  explained  in  terms  of  double 
donor  centers  interacting  with  the  valence  band.  We  conclude  that 
when  the  energy  levels  of  a defect  are  within  a few  KT  from  each 
other,  the  collected  data  for  a particular  energy  level  will  be  affected 
by  the  behavior  of  a nearby  energy  level. 

We  analyzed  the  noise  of  several  multi-level  trap  centers,  and 
derived  the  conditions  for  which  the  noise  of  a multi-level  trap  center 
reduces  to  the  noise  of  a single-level  trap  center.  For  a defect  with 
multi-energy  levels,  which  are  separated  by  more  than  a few  KT  from 
each  other,  the  multi-level  G-R  noise  model  reduces  to  the  single-level 
G-R  noise  model  in  both  uncompensated  and  overcompensated 
devices.  Intuitively,  this  is  consistent  with  the  fact  that  very  often  a 
single  level  G-R  noise  model  is  used  to  extract  defect  parameters.  We 
also  devised  a method  by  which  G-R  noise  measurements  can  be  used 
to  discriminate  between  donor-like  and  acceptor-like  trap  behavior. 
Using  a device  doped  with  an  intended  deep  defect  concentration 
larger  than  the  shallow  dopant  concentration,  we  can  determine 
whether  a defect  is  an  acceptor  or  a donor  center.  A deep  acceptor 
center  is  unlikely  to  be  detected  in  a p-type  device  because  the  Fermi 
level  is  far  away  from  the  deep  acceptor  energy  levels.  When  the  deep 
defect  is  donor-like,  the  Fermi  level  is  close  to  the  trap  energy  level 
and  the  noise  becomes  a strong  function  of  temperature  and  defect 
parameters. 
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Since  emission  and  capture  coefficients  may  be  a function  of 
electric  field,  and  the  free  carrier  concentration  in  a strongly  doped 
device  is  determined  by  the  balance  of  emission  and  capture  processes 
in  traps,  a field-induced  decrease  in  the  emission  coefficient  or  a 
field-enhanced  capture  coefficient,  or  both  can  reduce  the  free  carrier 
concentration  in  a semiconductor  device  when  subject  to  electric 
field.  The  current-voltage  (I-V)  characteristics  of  the  devices  used  in 
the  first  experiment  were  measured  and  they  showed  sublinear 
behavior  at  moderately  high  fields.  With  the  use  of  low  frequency 
voltage  noise  data  measured  as  a function  of  electric  field,  we  found 
that  the  sublinear  I-V  behavior  was  caused  by  a decrease  in  the  hole 
concentration  with  increasing  electric  field.  This  decrease  is  caused 
by  a field-induced  decrease  in  the  hole  emission  coefficient,  which  is 
larger  than  the  associated  field-induced  decrease  in  the  hole  capture 
coefficient. 

In  addition  to  the  above  described  macroscopic  G-R  noise 
mechanisms,  a microscopic  model  for  calculating  the  carrier  capture 
coefficient  for  deep  traps  is  derived.  In  our  model,  the  carrier  releases 
its  energy  through  delocalized  optical  phonon  emission  and  directly 
jumps  from  a free  vibronic  state  to  a bound  vibronic  state.  The 
calculated  capture  coefficient  shows  a slight  temperature  dependence 
and  agrees  with  the  experimental  data  measured  between  160  and 
200K. 

Theories  for  G-R  noise  generated  in  field-effect  devices  are 
discussed.  Trap  parameters  can  be  extracted  from  G-R  noise 
generated  in  junction  devices.  Its  interpretation  however  is  more 
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involved  than  is  the  case  for  bulk  devices.  Finally,  we  analyzed  G-R 
noise  models  associated  with  defects  in  the  grain  boundary  regions  of 
n-type  silicon  resistors. 

In  summary,  this  thesis  provides  comprehensive  information  on 
how  to  use  noise  data  from  G-R  processes  to  investigate  trap 
parameters. 


9.2.  Suggestions  for  Future  Research 

Firstly,  the  capture  coefficients  in  our  experiment  are  extracted 
using  optimization  techniques.  From  eq.  (2.9)  we  find  that  the  roll-off 
time  constant  provides  information  about  the  temperature- 
dependence  of  the  capture  coefficient.  One  can  reduce  the  doping 
concentration  of  shallow  acceptors  and  deep  impurities  and  increase 
the  device  cross  section  such  that  the  roll-off  time  constants  can  be 
measured.  Then  the  device  can  be  used  to  extract  directly  the 
information  about  the  temperature  and  field  dependence  of  capture 
coefficients. 

Secondly,  with  further  modifications  of  our  microscopic  capture 
model  developed  in  Chapter  VI  and  coupling  this  model  with  the 
Monte  Carlo  simulation  technique,  high  field  studies  of  G-R  processes 
and  carrier  transport  become  possible. 

Thirdly,  knowledge  of  the  field  dependence  of  capture  and 
emission  coefficients  is  very  useful  for  the  interpretation  of  capture 
and  emission  mechanisms  and  the  possible  deformation  of  bound 
carrier  potentials.  This  knowledge  can  be  obtained  by  measuring  the 
I-V  and  noise  characteristics  of  a symmetric  device  with  doping 
profiles  and  a geometry  that  eliminate  a bias  polarity  dependence  of 
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the  I-V  characteristics,  so  that  the  I-V  behavior  is  determined  only  by 
the  field-dependent  G-R  processes  and  the  electric  field  distribution. 

Finally,  we  developed  a simple  G-R  noise  model  for  polysilicon  in 
Chapter  VIII  based  on  the  symmetrical  Schottky-barrier  band  model. 
This  model  needs  further  modifications.  It  should  include  the 
asymmetrical  behavior  of  the  crystallites,  of  the  voltage  distribution, 
and  of  the  band  structure  of  grain  boundaries.  We  also  may  dope  an 
intended  defect  in  the  grain  boundary  such  that  the  Fermi  level  is 
nearby  the  defect  energy  level  and  is  therefore  known.  This  will  make 
it  easier  to  characterize  grain-boundary  defects.  Once  these  steps  have 
been  taken,  it  may  be  easier  to  find  a solution  for  the  D/(1-K04  tp4)-type 
noise  problem. 
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